Description:

CZ市为了欢迎全国各地的同学，特地举办了一场盛大的美食节。作为一个喜欢尝鲜的美食客，小M自然不愿意错过这场盛宴。他很快就尝遍了美食节所有的美食。然而，尝鲜的欲望是难以满足的。尽管所有的菜品都很可口，厨师做菜的速度也很快，小M仍然觉得自己桌上没有已经摆在别人餐桌上的美食是一件无法忍受的事情。于是小M开始研究起了做菜顺序的问题，即安排一个做菜的顺序使得同学们的等待时间最短。小M发现，美食节共有n种不同的菜品。每次点餐，每个同学可以选择其中的一个菜品。总共有m个厨师来制作这些菜品。当所有的同学点餐结束后，菜品的制作任务就会分配给每个厨师。然后每个厨师就会同时开始做菜。厨师们会按照要求的顺序进行制作，并且每次只能制作一人份。此外，小M还发现了另一件有意思的事情: 虽然这m个厨师都会制作全部的n种菜品，但对于同一菜品，不同厨师的制作时间未必相同。他将菜品用1, 2, ..., n依次编号，厨师用1, 2, ..., m依次编号，将第j个厨师制作第i种菜品的时间记为 ti,j 。小M认为：每个同学的等待时间为所有厨师开始做菜起，到自己那份菜品完成为止的时间总长度。换句话说，如果一个同学点的菜是某个厨师做的第k道菜，则他的等待时间就是这个厨师制作前k道菜的时间之和。而总等待时间为所有同学的等待时间之和。现在，小M找到了所有同学的点菜信息: 有 pi 个同学点了第i种菜品（i=1, 2, ..., n）。他想知道的是最小的总等待时间是多少。

Hint:

就是如何安排的问题,和修车那题一样,枚举所有二元组

不过由于这题数据加强了,一来不能加太多边,spfa会挂

考虑每次必定是从小到大找到一条容量为1的增广路

我们就从小到大动态加边就行了,这样不会超过$ p_sum $次

#include <map>
#include <set>
#include <stack>
#include <cmath>
#include <queue>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#define ls p<<1 
#define rs p<<1|1
using namespace std;
typedef long long ll;
const int mxn=1e5+5;
int S,T,n,m,cnt,ans,cost,sum,hd[mxn],tm[mxn],ck[mxn],p[805],a[805][805];
int pre[mxn],f[mxn],dis[mxn],vis[mxn],bl[mxn];

inline int read() {
    char c=getchar(); int x=0,f=1;
    while(c>'9'||c<'0') {if(c=='-') f=-1;c=getchar();}
    while(c<='9'&&c>='0') {x=(x<<3)+(x<<1)+(c&15);c=getchar();}
    return x*f;
}
inline void chkmax(int &x,int y) {if(x<y) x=y;}
inline void chkmin(int &x,int y) {if(x>y) x=y;}

struct ed {
    int to,nxt,w,val;
}t[mxn<<1];

inline void add(int u,int v,int w,int z) {
    t[cnt]=(ed) {v,hd[u],w,z}; hd[u]=cnt++;
    t[cnt]=(ed) {u,hd[v],0,-z}; hd[v]=cnt++;
}

int spfa() {
    memset(vis,0,sizeof(vis));
    memset(f,0x3f,sizeof(f));
    memset(dis,0x3f,sizeof(dis));
    memset(pre,0,sizeof(pre));
    queue<int > q; q.push(S); dis[S]=0;
    while(!q.empty()) {
        int u=q.front(); q.pop(); vis[u]=0;
        for(int i=hd[u];i!=-1;i=t[i].nxt) {
            int v=t[i].to;
            if(t[i].w>0&&dis[v]>dis[u]+t[i].val) {
                dis[v]=dis[u]+t[i].val; pre[v]=u;
                bl[v]=i; f[v]=min(f[u],t[i].w);
                if(!vis[v]) q.push(v),vis[v]=1;
            }
        }
    }
    return pre[T];
}

void Dinic() {
    while(spfa()) {
        ans+=f[T];
        cost+=f[T]*dis[T]; 
        int tp=t[bl[T]^1].to;
        add(tp+1,T,1,0);
        for(int i=1;i<=n;++i) add(i,tp+1,1,a[i][ck[tp+1]]*tm[tp+1]);
        for(int i=T;i!=S;i=pre[i]) {
            t[bl[i]].w-=f[T];
            t[bl[i]^1].w+=f[T];
        }
    }
}

int main()
{
    n=read(),m=read(); S=0; memset(hd,-1,sizeof(hd));
    for(int i=1;i<=n;++i) p[i]=read(),sum+=p[i]; T=sum*(m+1)+n+n+1;
    for(int i=1;i<=n;++i) add(S,i,p[i],0);  
    for(int i=1;i<=n;++i) 
        for(int j=1;j<=m;++j) 
            a[i][j]=read(),add(i,j*sum+1,1,a[i][j]);
    for(int i=1;i<=m;++i) add(i*sum+1,T,1,0);
    for(int i=1;i<=m;++i) 
        for(int j=1;j<=sum;++j) {
            int id=i*sum+j;
            tm[id]=j,ck[id]=i;
        }
    Dinic(); printf("%d",cost); 
    return 0;
}