[poj3614] Sunscreen

题意：

给你m头奶牛，n瓶防晒霜，奶牛们要晒太阳，每头奶牛所能承受的阳光为([l_i,r_i])，每瓶防晒霜的防晒值为(s_i)，可以给(c_i)头奶牛用，求最多有多少个奶牛能防晒。

题解：

最大流；

最开始是打算用堆去做的，第一眼看上去就像是二分图匹配，再看到(c_i)就很像是网络流了；

源点向防晒霜连(c_i)的边，防晒霜向奶牛连1的边，奶牛向汇点连1的边。

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<queue>
#define ll long long
#define N 10010
#define M 1000010
using namespace std;

int s,t,n,m,e_num=-1,inf=1<<30;
int nxt[M<<1],to[M<<1],w[M<<1],h[N],lev[N],cur[N],l[N],r[N],spf[N],cc[N];

queue<int> q;

int gi() {
  int x=0,o=1; char ch=getchar();
  while(ch!='-' && (ch<'0' || ch>'9')) ch=getchar();
  if(ch=='-') o=-1,ch=getchar();
  while(ch>='0' && ch<='9') x=x*10+ch-'0',ch=getchar();
  return o*x;
}

void add(int x, int y, int z) {
  nxt[++e_num]=h[x],to[e_num]=y,w[e_num]=z,h[x]=e_num;
}

bool bfs() {
  queue<int> q;
  for(int i=0; i<=t; i++) lev[i]=0;
  lev[s]=1,q.push(s);
  while(!q.empty()) {
    int u=q.front(); q.pop();
    for(int i=h[u]; i!=-1; i=nxt[i]) {
      int v=to[i];
      if(!lev[v] && w[i]>0) {
	lev[v]=lev[u]+1;
	q.push(v);
	if(v==t) return true;
      }
    }
  }
  return false;
}

int dfs(int u, int f) {
  if(u==t) return f;
  int tag=0,c;
  for(int &i=cur[u]; i!=-1; i=nxt[i]) {
    int v=to[i];
    if(lev[v]==lev[u]+1 && w[i]>0) {
      c=dfs(v,min(f-tag,w[i]));
      w[i]-=c,w[i^1]+=c,tag+=c;
      if(tag==f) return tag;
    }
  }
  return tag;
}

int Dinic() {
  int ans=0;
  while(bfs()) {
    for(int i=0; i<=t; i++) cur[i]=h[i];
    ans+=dfs(s,inf);
  }
  return ans;
}

int main() {
  memset(h,-1,sizeof(h));
  m=gi(),n=gi(),s=0,t=n+m+1;
  for(int i=1; i<=m; i++) l[i]=gi(),r[i]=gi();
  for(int i=1; i<=n; i++) spf[i]=gi(),cc[i]=gi();
  for(int i=1; i<=n; i++)
    for(int j=1; j<=m; j++) {
      if(spf[i]>=l[j] && spf[i]<=r[j]) add(i,j+n,1),add(j+n,i,0);
    }
  for(int i=1; i<=n; i++) add(s,i,cc[i]),add(i,s,0);
  for(int i=1; i<=m; i++) add(i+n,t,1),add(t,i+n,0);
  printf("%d", Dinic());
  return 0;
}