bzoj1570: [JSOI2008]Blue Mary的旅行(二分+网络流)

1570: [JSOI2008]Blue Mary的旅行

题目：传送门

---

题解：

　　 get到拆点新姿势，还是做题太少了...ORZ

　　 因为每天就只能有一个航班，那就不能直接连了，所以要拆点（然后就被卡住了）

　　 二分的话肯定是二分天数啦，这个就不说了

　　 拆点的方法：把每种航班拆成天数+1（太强辣%%%Star_Feel），因为第一个点的编号是1，而到达的天数为0

　　 那么这时候如果说x机场可以到y机场，那就把x机场的第i天连到y机场的第i+1天就好（代表花费了一天），流量就为机票数。

　　 特别神奇...然后一通乱check就没了...ORZ%%%

　　 注意答案求的是第几天，那么check的时候参数是总天数，那么就是mid+1

　　

---

代码：

 #include<cstdio>
 #include<cstring>
 #include<cstdlib>
 #include<cmath>
 #include<algorithm>
 #define inf 999999999
 using namespace std;
 struct node
 {
     int x,y,c,next,other;
 }a[];int len,last[];
 void ins(int x,int y,int c)
 {
     int k1,k2;
     k1=++len;
     a[len].x=x;a[len].y=y;a[len].c=c;
     a[len].next=last[x];last[x]=len;

     k2=++len;
     a[len].x=y;a[len].y=x;a[len].c=;
     a[len].next=last[y];last[y]=len;

     a[k1].other=k2;
     a[k2].other=k1;
 }
 int n,m,T,st,ed,head,tail;
 int list[],h[];
 bool bt_h()
 {
     memset(h,,sizeof(h));h[st]=;
     list[]=st;head=;tail=;
     while(head!=tail)
     {
         int x=list[head];
         for(int k=last[x];k;k=a[k].next)
         {
             int y=a[k].y;
             if(h[y]== && a[k].c>)
             {
                 h[y]=h[x]+;
                 list[tail++]=y;
             }
         }
         head++;
     }
     if(h[ed]>)return true;
     return false;
 }
 int find_flow(int x,int flow)
 {
     if(x==ed)return flow;
     int s=,t;
     for(int k=last[x];k;k=a[k].next)
     {
         int y=a[k].y;
         if(h[y]==h[x]+ && a[k].c> && s<flow)
         {
             s+=t=find_flow(y,min(a[k].c,flow-s));
             a[k].c-=t;a[a[k].other].c+=t;
         }
     }
     if(s==)h[x]=;
     return s;
 }
 struct edge
 {
     int x,y,c;
 }e[];
 bool check(int k)
 {
     len=;memset(last,,sizeof(last));
     for(int i=;i<=m;i++)
     {
         int x=e[i].x,y=e[i].y,c=e[i].c;
         for(int j=;j<k;j++)ins((x-)*k+j,(y-)*k+j+,c);
     }
     st=,ed=n*k+;
     for(int i=;i<=k;i++)ins(st,i,T),ins((n-)*k+i,ed,inf);
     int ans=;while(bt_h())ans+=find_flow(st,T);
     if(ans>=T)return true;
     return false;
 }
 int main()
 {
     scanf("%d%d%d",&n,&m,&T);for(int i=;i<=m;i++)scanf("%d%d%d",&e[i].x,&e[i].y,&e[i].c);
     int l=,r=n*T,ans=;
     while(l<=r)
     {
         int mid=(l+r)>>;
         if(check(mid+))ans=mid,r=mid-;
         else l=mid+;
     }
     printf("%d\n",ans);
     return ;
 }