【CODEVS】1034 家园

【算法】网络流-最大流(dinic)

【题解】

飞船有可承载人数限制，地球为源点，月球为汇点，人像水流一样从以飞船上限为容量的边流向汇点。

人在各站点都面临着上船与否的选择，难以用DP解决最优策略，于是这样的取舍问题可以使用网络流。

其实主要是看数据范围。

构图思路：

由于我们要计算时间，所以使变量T从1开始循环直到运送完毕。

对于每个时间相应加边，从而可以不用在意周期问题。

时间T，对于每个点增加一个T*22+i的分点，并从(T-1)分点向T分点连一条容量为inf的边。

对于飞行路径从(T-1)*22+u向T*22+v连一条容量为p[i]的边。

构图完毕。

最后推一下极限情况下最长的时间，超过了输出0即可。

tot=1，不是0，不是2！

要加当前弧优化，速度会快很多(101ms→2ms)！

#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
const int inf=0x3f3f3f3f;
struct edge{int from,v,flow;}e[];
int S=,T,tot=,first[],q[],a[][],cur[],d[],p[],n,m,k;
void insert(int u,int v,int w)
{
    tot++;e[tot].v=v;e[tot].flow=w;e[tot].from=first[u];first[u]=tot;
    tot++;e[tot].v=u;e[tot].flow=;e[tot].from=first[v];first[v]=tot;
}
bool bfs()
{
    memset(d,-,sizeof(d));
    int head=,tail=;q[]=;d[]=;
    while(head!=tail)
     {
         int x=q[head++];if(head>=)head=;
         for(int i=first[x];i;i=e[i].from)
          if(d[e[i].v]==-&&e[i].flow>)
           {
               q[tail++]=e[i].v;if(tail>=)tail=;
               d[e[i].v]=d[x]+;
           }
     }
    if(d[T*+n+]==-)return ;
    return ;
}
int dfs(int x,int a)
{
    if(x==T*+n+||a==)return a;
    int flow=,f;
    for(int& i=cur[x];i;i=e[i].from)
     if(d[e[i].v]==d[x]+&&(f=dfs(e[i].v,min(a,e[i].flow)))>)
      {
          e[i].flow-=f;
          e[i^].flow+=f;
          flow+=f;
          a-=f;
          if(a==)break;
      }
    return flow;
}
int main()
{
    scanf("%d%d%d",&n,&m,&k);
    for(int i=;i<=m;i++)
     {
         scanf("%d%d",&p[i],&a[i][]);
         for(int j=;j<=a[i][];j++)
          {
              scanf("%d",&a[i][j]);
              if(a[i][j]==-)a[i][j]=n+;
          }
     }
    for(T=;T<=;T++)
     {
         for(int i=;i<=n+;i++)
          {
              insert((T-)*+i,T*+i,inf);
          }
         for(int i=;i<=m;i++)
          {
              int u=(T-)*+a[i][(T-)%a[i][]+];
              int v=T*+a[i][T%a[i][]+];
              insert(u,v,p[i]);
          }
         while(bfs())
          {
              for(int i=;i<=T*+n+;i++)cur[i]=first[i];
              k-=dfs(,inf);
          }
         if(k<=){printf("%d",T);return ;}
     }
    printf("");
    return ;
}