hdu2883

这道题目比前面一道求满流的题目多了一点不一样，我理解就是离散化，因为si，ei的范围实在太大，但n的范围却不大，则

把所有si，ei列出来，排序，得到最多2*n-1个区间，则每个区间就是各节点，到汇点的容量限制就是m*区间长度。

#include<cstdio>#include<cstring>#include<algorithm>using namespace std;const int INF=1<<30;int si[210],ei[210],ni[210],ti[210],head[2100],h[2100],num[4200];int a[410],n,m,sum;int cnt,ans,c;struct edge{       int x,y,flow,nxt,op;}em[1000010];void add(int x,int y,int c){     em[++cnt].x=x;     em[cnt].y=y;     em[cnt].flow=c;     em[cnt].nxt=head[x];     head[x]=cnt;     em[cnt].op=cnt+1;     cnt++;     em[cnt].x=y;     em[cnt].y=x;     em[cnt].flow=0;     em[cnt].nxt=head[y];     head[y]=cnt;     em[cnt].op=cnt-1;}int dfs(int x,int flow){    if(x==n+c)    return flow;    int temp=flow,pos=n+c+1;    for(int j=head[x];j;j=em[j].nxt)    {            int y=em[j].y;            int w=em[j].flow;            if(h[x]==h[y]+1&&w>0)            {                                 int temp_flow=dfs(y,min(temp,w));                                 temp-=temp_flow;                                 em[j].flow-=temp_flow;                                 em[em[j].op].flow+=temp_flow;                                 if(!temp||h[0]==n+c+1)                                 return flow-temp;            }            if(w>0&&h[y]<pos)            pos=h[y];    }    if(temp==flow)    {                  num[h[x]]--;                  if(num[h[x]]==0)                  h[0]=n+c+1;                  else                  {                      h[x]=pos+1;                      num[h[x]]++;                  }    }    return flow-temp;}void isap(){     memset(h,0,sizeof(h));     memset(num,0,sizeof(num));     while(h[0]<n+c+1)     {                    ans+=dfs(0,INF);     } //    printf("%d %d\n",sum,ans);     if(ans==sum)     printf("Yes\n");     else     printf("No\n");}int main(){    while(2==scanf("%d%d",&n,&m))    {                                 sum=0;                                 int tot=0;                                 cnt=0;                                 memset(head,0,sizeof(head));                                 for(int i=1;i<=n;i++)                                 {                                         scanf("%d%d%d%d",si+i,ni+i,ei+i,ti+i);                                                                                  a[++tot]=si[i];                                         a[++tot]=ei[i];                                         add(0,i,ni[i]*ti[i]);                                         sum+=ni[i]*ti[i];                                 }                                 sort(a+1,a+tot+1);                                 c=0;                                 for(int i=1;i<=tot;i++)                                 if(a[c]!=a[i])                                 a[++c]=a[i];                               /*  for(int i=1;i<=c;i++)                                 printf("%d\n",a[i]);*/                                 for(int i=1;i<=c-1;i++)                                 {                                         int num=(a[i+1]-a[i])*m;                                         add(n+i,n+c,num);                                 }                                 for(int i=1;i<=c-1;i++)                                 {                                         for(int j=1;j<=n;j++)                                         if(si[j]<=a[i]&a[i+1]<=ei[j])                                         add(j,i+n,INF);                                 }                                  ans=0;                                 isap();    }        return 0;}