【最大流】BAPC2014 A Avoiding the Apocalypse (Codeforces GYM 100526)

题目链接：

　　http://codeforces.com/gym/100526

　　http://acm.hunnu.edu.cn/online/?action=problem&type=show&id=11664&courseid=0

题目大意：

　　总共N个点，现在有P个人，在T时间内要从起点S走到C个医疗站中的任意一个，M条X到Y的有向边，每条边每个单位时间可以通过的人数是pi,走完这条边耗时ti。

　　人可以停留在任意一个点，求最多能有多少人到达医疗站。

　　(1<=C<=N<=1000,0<=M<=1000,1<=P,T<=100,1<=pi,ti<=100)

题目思路：

　　【最大流】

　　SPFA居然T了。。

　　将每个点按照时间拆点{X节点j时间的状态为点X(j)}，把X到Y的边pi,ti拆成T-ti条边，每条边从X(j)到Y(j+ti)容量为pi，j为时间从0~T-ti，表示在j时间有pi个人可以从X到Y，到达Y的时间为j+ti。

　　因为人可以停留在节点上所以X(i)到X(i+1)连一条容量P的边。

　　每个医疗站的最终时间Ci(T)到超级汇T_T连一条容量为P的边（到达医疗站的人就无限停留直到时间为T到超级汇T_T）

　　超级源S_S到起点S(0)连一条容量为P的边表示初始起点S有P个人。

　　最终统计到达超级汇的人数即可。跑一遍最大流。

 //
 //by coolxxx
 //#include<bits/stdc++.h>
 #include<iostream>
 #include<algorithm>
 #include<string>
 #include<iomanip>
 #include<map>
 #include<memory.h>
 #include<time.h>
 #include<stdio.h>
 #include<stdlib.h>
 #include<string.h>
 //#include<stdbool.h>
 #include<math.h>
 #define min(a,b) ((a)<(b)?(a):(b))
 #define max(a,b) ((a)>(b)?(a):(b))
 #define abs(a) ((a)>0?(a):(-(a)))
 #define lowbit(a) (a&(-a))
 #define sqr(a) ((a)*(a))
 #define swap(a,b) ((a)^=(b),(b)^=(a),(a)^=(b))
 #define mem(a,b) memset(a,b,sizeof(a))
 #define eps (1e-8)
 #define J 10
 #define mod 1000000007
 #define MAX 0x7f7f7f7f
 #define PI 3.14159265358979323
 #define N 1004
 #define M 104
 using namespace std;
 typedef long long LL;
 int cas,cass;
 int n,m,lll,ans;
 int s,p,t,S,T,nn;
 int c[N];
 int last[N*M+M],vd[N*M+M],d[N*M+M];
 struct xxx
 {
     int next,to,q;
 }a[N*M*M];
 void add(int x,int y,int f)
 {
     a[++lll].next=last[x];
     a[lll].to=y;
     a[lll].q=f;
     last[x]=lll;
 }
 int sap(int u,int f)
 {
     int i,v,tt,asp=,mix=nn-;
     if(u==T)return f;
     for(i=last[u];i;i=a[i].next)
     {
         v=a[i].to;
         if(a[i].q>)
         {
             if(d[u]==d[v]+)
             {
                 tt=sap(v,min(f-asp,a[i].q));
                 asp+=tt;
                 a[i].q-=tt;
                 a[i^].q+=tt;
                 if(asp==f || d[S]==nn)
                     return asp;
             }
             mix=min(mix,d[v]);
         }
     }
     if(asp!=)return asp;
     if(!--vd[d[u]])d[S]=nn;
     else vd[d[u]=mix+]++;
     return asp;
 }
 int main()
 {
     #ifndef ONLINE_JUDGE
 //    freopen("1.txt","r",stdin);
 //    freopen("2.txt","w",stdout);
     #endif
     int i,j,k;
     int x,y,f,dd;
     for(scanf("%d",&cas);cas;cas--)
 //    for(scanf("%d",&cas),cass=1;cass<=cas;cass++)
 //    while(~scanf("%s",s+1))
 //    while(~scanf("%d",&n))
     {
         ans=;lll=;mem(vd,);mem(d,);mem(last,);
         scanf("%d%d%d%d%d",&n,&s,&p,&t,&cass);
         S=n*t+t+;T=S+;
         nn=T;vd[]=nn-t+;
         add(S,t,p);
         add(t,S,);
         add(t,s*t+,p);
         add(s*t+,t,);
         for(i=;i<=n;i++)
         {
             for(j=;j<t;j++)
             {
                 add(i*t+j,i*t+j+,p);
                 add(i*t+j+,i*t+j,);
             }
         }
         for(i=;i<=cass;i++)
         {
             scanf("%d",&c[i]);
             add(c[i]*t+t,T,p);
             add(T,c[i]*t+t,);
         }
         scanf("%d",&m);
         for(i=;i<=m;i++)
         {
             scanf("%d%d%d%d",&x,&y,&f,&dd);
             for(j=;j+dd<=t;j++)
             {
                 add(x*t+j,y*t+j+dd,f);
                 add(y*t+j+dd,x*t+j,);
             }
             if(x==s)
             {
                 add(t,y*t+dd,f);
                 add(y*t+dd,t,);
             }
         }
         while(d[S]<nn)
         {
             f=sap(S,MAX);
             ans+=f;
         }
         printf("%d\n",ans);
     }
     return ;
 }
 /*
 //
 
 //
 */