BZOJ2879 [Noi2012]美食节
AC通道:http://www.lydsy.com/JudgeOnline/problem.php?id=2879
这题codevs上也有,不过数据不同:http://codevs.cn/problem/1935/
如果你觉得你的程序跑得很快,那么你可以交codevs,如果你觉得你的程序可能是错的,那么请你交BZOJ
BZOJ数据貌似厉害一些?不过BZOJ时限也放得开些。
这题就是上题的升级版,在修车的基础上,数据很大,所以需要通过动态加边来使其变快。
因为上题也分析到了,只有厨师走完短的路,才会走长的,那么如果没走短的路,我长的路留在这里只是增加spfa的复杂度。
然后动态加边就可以了。
这个题刚开始的思路是错的,就是我在动态加边的同时,可以动态加点,但这样需要一个数组记下这是哪个厨师,也要记下这是厨师的第几个菜。
但是发现如果我的费用流需要退流的话,我记录第几个菜的数组需要-1,但是我做不到,这会导致这个点被重复算几次,然后就比答案小了。
所以必须确定下位置,这样的话,你就不会重复的添加点了,最多重复添加相同意义的边,而因为流量限制,这条边没有实际意义。
不过这样的话,速度就会变慢...
好慢啊好慢啊...
![](https://images.cnblogs.com/OutliningIndicators/ContractedBlock.gif)
#include<queue> #include<cstdio> #include<cstring> #include<algorithm> using namespace std; const int maxn=50; const int maxm=110; const int maxp=100010; const int INF=0x3f3f3f3f; struct Node{ int data,next,low,cost; }node[maxp*maxn*2]; #define now node[point].data #define www node[point].low #define ccc node[point].cost #define then node[point].next int n,m,cnt,ans; int s,t,Idex,tot; int head[maxp]; int c[maxp],ind[maxp],ink[maxp]; int dis[maxp],pre[maxp]; int a[maxm][maxm]; bool inq[maxp]; deque<int >q; void add(int u,int v,int w,int c){ node[cnt].data=v;node[cnt].next=head[u];node[cnt].low=w;node[cnt].cost=c;head[u]=cnt++; node[cnt].data=u;node[cnt].next=head[v];node[cnt].low=0;node[cnt].cost=-c;head[v]=cnt++; } bool SPFA(){ int x,Low=INF,ta,tb; memset(dis,0x3f,sizeof(dis)); q.push_back(s);dis[s]=0;pre[s]=pre[t]=-1; while(!q.empty()){ x=q.front();q.pop_front();inq[x]=false; for(int point=head[x];point!=-1;point=then) if(www && dis[now]>dis[x]+ccc){ dis[now]=dis[x]+ccc;pre[now]=point^1; if(!inq[now]){ inq[now]=true; if(q.empty() || dis[q.front()]>dis[now]) q.push_front(now); else q.push_back(now); } } } if(pre[t]<0) return false; for(int point=pre[t];point!=-1;point=pre[now]){ if(now==s){ x=node[point^1].data; ta=(x-1)/tot+1;tb=x%tot+1; } Low=min(Low,min(Low,node[point^1].low)); } for(int point=pre[t];point!=-1;point=pre[now]) node[point^1].low-=Low,www+=Low; for(int i=1;i<=m;i++) add((ta-1)*tot+tb,n*tot+i,1,tb*a[ta][i]); ans+=Low*dis[t]; return true; } int main(){ #ifndef ONLINE_JUDGE freopen("2879.in","r",stdin); freopen("2879.out","w",stdout); #endif scanf("%d%d",&m,&n); for(int i=1;i<=m;i++) scanf("%d",&c[i]),tot+=c[i]; for(int i=1;i<=m;i++) for(int j=1;j<=n;j++) scanf("%d",&a[j][i]); t=n*tot+m+1; for(int i=s;i<=t;i++) head[i]=-1; for(int i=1;i<=n*tot;i++) add(s,i,1,0); for(int i=1;i<=m;i++) add(n*tot+i,t,c[i],0); for(int i=1;i<=n;i++) for(int k=1;k<=m;k++) add((i-1)*tot+1,n*tot+k,1,a[i][k]); while(SPFA()); printf("%d",ans); return 0; }