2879: [Noi2012]美食节

题目链接

题目大意：n个厨师，做m道菜有不同时间，最小化总等待时间

题解：和1070一样……但是这题数据范围太大了会TLE
神奇优化：动态加边，只保留当前有用的边
某位厨师在做倒数第k道菜后
再把他做倒数第k+1道菜的边加上

我的收获：动态加边黑科技……

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
using namespace std;
const int M=105050;
#define INF 0x3f3f3f3f

int n,m,t,st,ed,ans,tot,nod;
int w[M];
bool vis[M];
int head[M],d[M],pre[M],T[65][105];

queue<int> q;

struct edge{int fro,to,nex,c,val;}e[3000005];

void add(int i,int j,int ca,int co){e[t].fro=i,e[t].to=j,e[t].nex=head[i],e[t].c=ca,e[t].val=co;head[i]=t++;}
void insert(int i,int j,int w,int z){add(i,j,w,z),add(j,i,0,-z);}

bool spfa()
{
    for(int i=0;i<=ed;i++) vis[i]=0,d[i]=INF;
    d[st]=0;q.push(st);
    while(!q.empty())
    {
        int u=q.front();q.pop();vis[u]=false;
        for(int i=head[u];i!=-1;i=e[i].nex){
            int v=e[i].to;
            if(e[i].c&&d[v]>d[u]+e[i].val){
                d[v]=d[u]+e[i].val;
                pre[v]=i;
                if(!vis[v]){
                    vis[v]=true;
                    q.push(v);
                }
            }
        }
    }
    return d[ed]!=INF;
}
void flow()
{
    int mx=INF;
    for(int u=ed;u!=st;u=e[pre[u]].fro)
        mx=min(mx,e[pre[u]].c);
    for(int u=ed;u!=st;u=e[pre[u]].fro){
        e[pre[u]].c-=mx;e[pre[u]^1].c+=mx;
        ans+=mx*e[pre[u]].val;
    }
    int x=e[pre[ed]].fro;  
    int a=(x-n-1)/tot+1,b=(x-n)%tot+1;  
    if(b<=tot){  
        for(int i=1;i<=n;i++){  
            insert(i,n+(a-1)*tot+b,1,b*T[i][a]);  
        }  
    }  
    insert(n+(a-1)*tot+b,ed,1,0);  
}

void work()
{
    while(spfa()) flow();
    cout<<ans<<endl;
}

void build()
{
    for(int i=1;i<=n;i++){  
        for(int j=1;j<=m;j++)
            insert(i,n+(j-1)*tot+1,1,T[i][j]);  
    }  
    for(int i=1;i<=m;i++) insert(n+(i-1)*tot+1,ed,1,0);  
}

void init()
{
    t=0;memset(head,-1,sizeof(head));
    cin>>n>>m;
    for(int i=1;i<=n;i++){
        cin>>w[i],tot+=w[i];
        insert(0,i,w[i],0);
    }
    st=0,ed=m*tot+n+1;
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            scanf("%d",&T[i][j]);
    build();
}

int main()
{
    init();
    work();
    return 0;
}