记录编号 412444 评测结果 AAAAAAAAAA
题目名称 [NOI 2012]美食节 最终得分 100
用户昵称 GravatarCooook 是否通过 通过
代码语言 C++ 运行时间 9.914 s
提交时间 2017-06-09 11:11:43 内存使用 207.30 MiB
显示代码纯文本
#include <cstdio>
#include <cstring>
#include <queue>
#include <iostream>
using namespace std;
const int MAXN = 1000005,MAXN_edge = 10000005;
const int inf = 0x3f3f3f3f;
int n,m,w[105][105],num[45],s,t;
int first[MAXN],e=2,dis[MAXN],tot,pre[MAXN];
int flow[MAXN];
int ans = 0;
bool vis[MAXN];


struct edge{
    int u,v,next,w,c;
}a[MAXN_edge];

inline void add(int u,int v,int w,int c){
    a[e].u=u;
    a[e].v=v;
    a[e].w=w;
    a[e].c=c;
    a[e].next=first[u];
    first[u]=e++;
}

inline void push(int u,int v,int w,int c){
    add(u,v,w,c);add(v,u,0,-c);
}

inline bool bfs(){
    memset(flow,0x3f,sizeof(flow));
    memset(dis,0x3f,sizeof(dis));
    memset(pre,0,sizeof(pre));
    dis[s]=0;
    vis[s]=1;
    queue<int>q;
    q.push(s);
    while(!q.empty()){
        int k=q.front();
        q.pop();vis[k]=0;
        for(int i=first[k];i;i=a[i].next){
            int v=a[i].v;
            if(dis[v]>dis[k]+a[i].c && a[i].w){
                dis[v]=dis[k]+a[i].c;
                pre[v]=i;
                flow[v]=min(flow[k],a[i].w);
                if(!vis[v])vis[v]=1,q.push(v);
            }
        }
    }
    return dis[t]!=0x3f3f3f3f;
}

inline void change(){
    ans+=dis[t]*flow[t];
    for(int i=pre[t];i;i=pre[a[i].u]){
    	a[i].w-=flow[t];
    	a[i^1].w+=flow[t];
	}
    int p = a[pre[t]].u/tot,b=a[pre[t]].u%tot;
    for(int i=1;i<=n;i++)
        push(i,(1+p)*tot+b,1,w[b][i]*(p+1));
    push((1+p)*tot+b,t,1,0);
}


int main(){

    freopen("noi12_delicacy.in","r",stdin);

    freopen("noi12_delicacy.out","w",stdout);

    scanf("%d%d",&n,&m);
    

    for(int i=1;i<=n;i++)scanf("%d",&num[i]),tot+=num[i];
    
    t=tot*tot+tot+1;

    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            scanf("%d",&w[j][i]);

    for(int i=1;i<=n;i++) push(s,i,num[i],0);
    
    for(int i=1;i<=m;i++)push(tot+i,t,1,0);

    for(int i=1;i<=n;i++)
        for(int k=1;k<=m;k++){
            push(i,tot+k,1,w[k][i]);
		}
        
	while(bfs())change();
    
    printf("%d\n",ans);
}