博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
bzoj 1070 费用流
阅读量:4541 次
发布时间:2019-06-08

本文共 2068 字,大约阅读时间需要 6 分钟。

//可以网络流,但是要怎么分配每辆车让谁维修以及维修顺序呢。可以考虑每辆车维修时间对总结果的贡献,把每个修车人拆成n个点共n*m个点,//n辆车连向这n*m个点,流量1,费用k*修车时间,其中k(1=
<=n)表示这辆车是这个技术人员修的倒数第k俩车,他后面k-1辆被同一个人维修的车//都要等待他修好,因此他的贡献值就是k*维修时间。然后源点连向n辆车,n*m个点连向汇点。#include
#include
#include
#include
using namespace std;const int INF=0x7fffffff;const int MAXN=1000;const int MAXM=1000000;int N,a[70][12],tot,head[MAXN+10],pre[MAXN],dis[MAXN];bool vis[MAXN+10];struct Edge{ int to,next,cap,flow,cost;}edge[MAXM+10];void init(int x){ N=x; tot=0; memset(head,-1,sizeof(head));}void add(int x,int y,int w,int c){ edge[tot].to=y; edge[tot].cap=w; edge[tot].cost=c; edge[tot].flow=0; edge[tot].next=head[x]; head[x]=tot++; edge[tot].to=x; edge[tot].cap=0; edge[tot].cost=-c; edge[tot].flow=0; edge[tot].next=head[y]; head[y]=tot++;}bool Spfa(int s,int t){ queue
q; for(int i=0;i<=N;i++){ vis[i]=0; dis[i]=INF; pre[i]=-1; } vis[s]=1; dis[s]=0; q.push(s); while(!q.empty()){ int u=q.front();q.pop(); vis[u]=0; for(int i=head[u];i!=-1;i=edge[i].next){ int v=edge[i].to; if(edge[i].cap>edge[i].flow&&dis[v]>dis[u]+edge[i].cost){ dis[v]=dis[u]+edge[i].cost; pre[v]=i; if(!vis[v]){ vis[v]=1; q.push(v); } } } } if(pre[t]==-1) return 0; return 1;}int MaxCostFlow(int s,int t){ int Flow=0,Cost=0; while(Spfa(s,t)){ int Min=INF; for(int i=pre[t];i!=-1;i=pre[edge[i^1].to]) Min=min(Min,edge[i].cap-edge[i].flow); for(int i=pre[t];i!=-1;i=pre[edge[i^1].to]){ edge[i].flow+=Min; edge[i^1].flow-=Min; Cost+=edge[i].cost*Min; } Flow+=Min; } return Cost;}int main(){ int n,m; while(scanf("%d%d",&m,&n)==2){ init(n+n*m+1); for(int i=1;i<=n;i++){ add(0,i,1,0); for(int j=1;j<=m;j++) scanf("%d",&a[i][j]); } for(int i=1;i<=n;i++){ for(int j=1;j<=m;j++){ for(int k=1;k<=n;k++){ add(i,j*n+k,1,k*a[i][j]); } } } for(int i=n+1;i<=n+n*m;i++) add(i,n+n*m+1,1,0); printf("%.2lf\n",1.0*MaxCostFlow(0,n+n*m+1)/(1.0*n)); } return 0;}

 

转载于:https://www.cnblogs.com/--ZHIYUAN/p/7469704.html

你可能感兴趣的文章
如何将一段文本编译成C#内存程序的过程
查看>>
PAT——1070. 结绳
查看>>
【23.33%】【codeforces 664C】International Olympiad
查看>>
java-网络编程-使用URLDecoder和URLEncoder
查看>>
最短路之dijkstra算法
查看>>
SHDP--Working With HBase (二)之HBase JDBC驱动Phoenix与SpringJDBCTemplate的集成
查看>>
Lua语法基础(一)
查看>>
.Net Core2.*学习手册
查看>>
实验一、命令解释程序的编写实验
查看>>
2018年11月14日 学习字符串用法2
查看>>
2019年5月26日 re模块2
查看>>
Python学习笔记(一)——初学Python
查看>>
顺序表应用8:最大子段和之动态规划法(SDUT 3665)
查看>>
Python内置函数(52)——range
查看>>
正则表达式
查看>>
c# 执行 sql service 的存储过程
查看>>
《软件构架实践》读后感03
查看>>
jQuery入门(4)jQuery中的Ajax应用
查看>>
Java 发送http GET/POST请求
查看>>
索引之详解一
查看>>