HDU 3315 My Brute(费用流)
来源:互联网 发布:java gbk转utf8 乱码 编辑:程序博客网 时间:2024/06/05 05:17
题目地址:HDU 3315
这个题的思路完全是自己想出来的,自我感觉挺巧妙的。。。(大牛勿喷。。。)对大胆建图又多了一份信心。
具体思路是构造一个二分图,Si连源点,Xi连汇点,流量都是1,费用0.然后当Si可以赢Xj的时候,就对这两人连一条边,费用值为-Vi*1000,如果i==j的话,费用值就再减1,因为题目要求尽量不改变原先的顺序,所以说应该尽量让序号相同的对打。而费用值减1的话,会优先考虑序号相同的,而且让费用扩大了1000倍,此时也不会改变主要的分数因素大小。同理,输的话,费用值为Vi*1000,如果i==j的话,费用值同样减1。
最后算出的最大费用cost/1000就是正确的费用值,cost%1000就是改变顺序了的数目。然后做相应判断与计算即可。
代码如下:
#include <iostream>#include <cstdio>#include <string>#include <cstring>#include <stdlib.h>#include <math.h>#include <ctype.h>#include <queue>#include <map>#include <set>#include <algorithm>using namespace std;const int INF=0x3f3f3f3f;int head[1000], source, sink, cnt, fei[10000], q[10000], flow, cost;int d[1000], vis[1000], cur[1000], f[1000];struct node{ int u, v, cap, cost, next;}edge[1000000];void add(int u, int v, int cap, int cost){ edge[cnt].v=v; edge[cnt].cap=cap; edge[cnt].cost=cost; edge[cnt].next=head[u]; head[u]=cnt++; edge[cnt].v=u; edge[cnt].cap=0; edge[cnt].cost=-cost; edge[cnt].next=head[v]; head[v]=cnt++;}int spfa(){ memset(d,INF,sizeof(d)); memset(vis,0,sizeof(vis)); queue<int>q; q.push(source); d[source]=0; f[source]=INF; cur[source]=-1; 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].v; if(d[v]>d[u]+edge[i].cost&&edge[i].cap) { d[v]=d[u]+edge[i].cost; f[v]=min(f[u],edge[i].cap); cur[v]=i; if(!vis[v]) { vis[v]=1; q.push(v); } } } } if(d[sink]==INF) return 0; flow+=f[sink]; cost-=f[sink]*d[sink]; for(int i=cur[sink];i!=-1;i=cur[edge[i^1].v]) { edge[i].cap-=f[sink]; edge[i^1].cap+=f[sink]; } return 1;}void mcmf(int n){ flow=cost=0; while(spfa()) ; if(cost/1000<=0) { printf("Oh, I lose my dear seaco!\n"); //printf("%d\n",cost); } else { int x=cost/1000, y=cost%1000; printf("%d %.3lf%%\n",x,y*100.0/n); }}int pan(int x1, int x2, int y1, int y2){ while(1) { x2-=y1; if(x2<=0) return 1; x1-=y2; if(x1<=0) return 0; }}int main(){ int n, i, V[100], H[100], P[100], A[100], B[100], j; while(scanf("%d",&n)!=EOF&&n) { for(i=1;i<=n;i++) scanf("%d",&V[i]); for(i=1;i<=n;i++) scanf("%d",&H[i]); for(i=1;i<=n;i++) scanf("%d",&P[i]); for(i=1;i<=n;i++) scanf("%d",&A[i]); for(i=1;i<=n;i++) scanf("%d",&B[i]); memset(head,-1,sizeof(head)); cnt=0; source=0; sink=2*n+1; for(i=1;i<=n;i++) { add(source,i,1,0); add(i+n,sink,1,0); } for(i=1;i<=n;i++) { for(j=1;j<=n;j++) { if(pan(H[i],P[j],A[i],B[j])) { if(i==j) add(i,j+n,1,-V[i]*1000-1); else add(i,j+n,1,-V[i]*1000); } else { if(i==j) add(i,j+n,1,V[i]*1000-1); else add(i,j+n,1,V[i]*1000); } } } mcmf(n); } return 0;}
1 0
- HDU 3315 My Brute(费用流)
- HDU 3315 My Brute(费用流)
- HDU 3315 My Brute 费用流
- HDU 3315 My Brute(费用流)
- HDU 3316 My Brute(二维费用流)经典
- HDU 3315 My Brute 费用流+求错位率(好题)
- hdu 3315 My Brute 费用流,费用最小且代价最小
- hdu 3315 My Brute 网络流
- HDU 2853 Assignment & HDU 3315 My Brute
- HDU 3315 My Brute(KM最大匹配)
- HDOJ 3315 - My Brute 维护两个最值..构图最大费用最大流 or KM模板
- HDU 3315 My Brute KM算法
- HDU 3315 My Brute(KM算法)
- My Brute (hdu 3315 二分图最大权匹配KM算法)
- HDU 3315 My Brute(二分图最优匹配)
- 【网络流】hdu3315 My Brute
- My Brute
- HDU 3315 My Brute(二分图最优匹配:优先用原匹配边)
- django中的时区设置TIME_ZONE,USE_TZ
- 黑马程序员----Sizeof与Strlen的区别与联系
- 阿里巴巴2014前端线上笔试题
- 男生对女生最具诱惑的10句话
- How to map android keys to your usb/bluetooth keyboard
- HDU 3315 My Brute(费用流)
- 使用IDEA创建Scala程序
- 剑指Offer18树的子结构
- Java 不同类型属性在字节码文件class中
- EventDispatcher Mechanism
- 杀人事件测你的潜意识??
- impdp 导入物化视图(mv)时报错 :ORA-31685: Object type MATERIALIZED_VIEW
- linux下如何恢复rm命令删除的文件
- How to Enable Multi-Touch