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;
}