BZOJ 1570: [JSOI2008]Blue Mary的旅行(迭代加深最大流)

版权声明：本文为博主原创文章，未经博主允许必须转载。 https://blog.csdn.net/qq_35950004/article/details/88416476

题目
好裸啊。
看到数据范围秒懂：每天拆点就行了。

AC Code:

#include<bits/stdc++.h>
#define maxn 10000
#define maxm maxn * 105
#define inf 0x3f3f3f3f
using namespace std;

int n,m,t;
int mark[305][105],x[maxm],y[maxm],w[maxm];
int S,T,tot,dis[maxn];

int info[maxn],Prev[maxm],to[maxm],cap[maxm],cnt_e=1;
inline void Node(int u,int v,int c){ Prev[++cnt_e]=info[u],info[u]=cnt_e,to[cnt_e]=v,cap[cnt_e]=c; }
inline void Line(int u,int v,int c){ Node(u,v,c),Node(v,u,0); }

int aug(int now,int Max)
{
    if(now == T) return Max;
    int inc , st = Max;
    for(int i=info[now];i;i=Prev[i])
        if(cap[i] && dis[to[i]] == dis[now] - 1)
        {
            inc = aug(to[i] , min(cap[i] , st));
            if(inc) st-= inc, cap[i]-=inc , cap[i^1] += inc;
            else dis[to[i]] = -1;
            if(!st) break;
        }
    return Max - st;
}

bool BFS()
{
    static queue<int>q;
    memset(dis,-1,sizeof dis);
    q.push(T),dis[T] = 0;
    for(int now;!q.empty();)
    {
        now = q.front() , q.pop();
        for(int i=info[now];i;i=Prev[i])
            if(cap[i^1] && dis[to[i]] == -1)
            {
                dis[to[i]] = dis[now] + 1;
                q.push(to[i]);
            }
    }
    return dis[S] != -1;
}

int main()
{
    scanf("%d%d%d",&n,&m,&t);
    for(int i=1;i<=m;i++)
        scanf("%d%d%d",&x[i],&y[i],&w[i]);
    S = ++ tot , T = ++ tot;
    int ans = 0 , stm = 0;
    Line(S,3,t);
    for(;;ans++)
    {
        for(int i=1;i<=n;i++) mark[ans][i] = ++tot;
        Line(mark[ans][n],T,t);
        if(ans)
        {
            for(int i=1;i<=m;i++)
                Line(mark[ans-1][x[i]],mark[ans][y[i]],w[i]);
            for(int i=1;i<=n;i++)
                Line(mark[ans-1][i],mark[ans][i],inf);
        }
        for(;BFS();)
            stm += aug(S,inf);
        if(stm == t) break;
    }
    printf("%d\n",ans);
}