首页 > 代码库 > hdu4780 费用流 (机器任务工作不中断问题)
hdu4780 费用流 (机器任务工作不中断问题)
机器,任务 ,每个任务有有时间,不可中断。
题意:m个机器,n个糖果要加工,给出每个糖果的工作时间(s,t),以及糖果之间、机器预备时间以及费用,求最小费用。
这题开始受原来可以时间中断那题影响,开始用时间建图,巨麻烦,后来学习了,才觉悟时间只是计算费用的,没有帮毛钱关系, s-->机器-》糖果-》t;
因为要每个糖果都加工一次,糖果拆点,必经过(-inf)。有一个没过,则无解。
网上的添加一组新点法也可以。我的使用释放压力+无穷碧流法也好。
跪点:注意!用inf=0x3f3f3f3f,最后加上n*inf,爆int! 这种以后要注意!顺变 e[][3]也要longlong !
本来这样我就A了!但是竟然因为一个变量的命名搞混了WA一天!!!SB了! 自己搞变量不要搞这么相似啊 !!!!!!!!!!
#include<iostream> #include<cstdio> #include<queue> #include<cstring> #include<algorithm> using namespace std; const long long inf=0x3f3f3f3f; const int maxn=505,maxe=80000; int n,m,k; int ss,tt; int head[maxn];long long e[maxe][4];int nume=0; void inline adde(int i, int j, int c,long long w) { e[nume][0]=j;e[nume][1]=head[i];head[i]=nume; e[nume][2]=c;e[nume++][3]=w; e[nume][0]=i;e[nume][1]=head[j];head[j]=nume; e[nume][2]=0;e[nume++][3]=-w; } int inq[maxn];long long d[maxn];int pre[maxn];int prv[maxn]; bool spfa(long long &sum) { for(int i=0;i<=tt;i++) { inq[i]=0;d[i]=inf; } queue<int>q; q.push(ss); inq[ss]=1;d[ss]=0; while(!q.empty()) { int cur=q.front(); q.pop(); inq[cur]=0; for(int j=head[cur];j!=-1;j=e[j][1]) { int v=e[j][0]; if(d[v]>d[cur]+e[j][3]&&e[j][2]>0) { d[v]=d[cur]+e[j][3]; pre[v]=j; prv[v]=cur; if(!inq[v]) { q.push(v); inq[v]=1; } } } } if(d[tt]==inf)return 0; int minf=inf; int cur=tt; while(cur!=ss) { minf=minf<e[pre[cur]][2]?minf:e[pre[cur]][2]; cur=prv[cur]; } cur=tt; while(cur!=ss) { e[pre[cur]][2]-=minf; e[pre[cur]^1][2]+=minf; cur=prv[cur]; } sum+=d[tt]*minf; return 1; } long long mincost() { long long sum=0; while(spfa(sum)); return sum; } int mac[105][2];int cmtime[105][105];int cmony[105][105]; int cgtime[105][105];int cgmony[105][105]; void build() { for(int i=1;i<=m;i++) { adde(ss,i,1,0); adde(i,tt,1,0); for(int j=1;j<=n;j++) { if(cmtime[j][i]<mac[j][1]) { if(cmtime[j][i]<=mac[j][0]) adde(i,j+m,1,cmony[j][i]); else adde(i,j+m,1,cmony[j][i]+k*(cmtime[j][i]-mac[j][0])); } } } for(int i=1;i<=n;i++) { adde(i+m,i+m+n,1,-inf); adde(i+m+n,tt,1,0); for(int j=1;j<=n;j++) { if(i!=j&&mac[j][1]>mac[i][1]+cgtime[i][j]) { if(mac[i][1]+cgtime[i][j]<=mac[j][0]) adde(i+m+n,j+m,1,cgmony[i][j]); else adde(i+n+m,j+m,1,cgmony[i][j]+k*(mac[i][1]+cgtime[i][j]-mac[j][0])); } } } } void init() { nume=0; memset(head,-1,sizeof(head)); ss=0,tt=m+n+n+1; } int main() { while(~scanf("%d%d%d",&n,&m,&k)&&(n||m||k)) { init(); for(int i=1;i<=n;i++) scanf("%d%d",&mac[i][0],&mac[i][1]); for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) { scanf("%d",&cmtime[i][j]); } for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) { scanf("%d",&cmony[i][j]); } for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) { scanf("%d",&cgtime[i][j]); } for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) { scanf("%d",&cgmony[i][j]); } build(); long long ans=mincost()+n*inf; if(ans>=inf) printf("-1\n"); else printf("%I64d\n",ans); } }
hdu4780 费用流 (机器任务工作不中断问题)
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。