HDU 3315 My Brute(费用流)
职务地址: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;
}
版权声明:本文博客原创文章。博客,未经同意,不得转载。
HDU 3315 My Brute(费用流)的更多相关文章
- My Brute HDU - 3315(KM || 费用流)
题意: 有S1到Sn这n个勇士要和X1到Xn这n个勇士决斗,初始时,Si的决斗对象是Xi. 如果Si赢了Xi,那么你将获得Vi分,否则你将获得-Vi分. Si和Xi对决时,Si有初始生命Hi,初始攻击 ...
- HDU 5988 Coding Contest(费用流+浮点数)
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=5988 题目大意: 给定n个点,m条有向边,每个点是一个吃饭的地方,每个人一盒饭.每个点有S个人,有B盒 ...
- POJ 2135 Farm Tour && HDU 2686 Matrix && HDU 3376 Matrix Again 费用流求来回最短路
累了就要写题解,近期总是被虐到没脾气. 来回最短路问题貌似也能够用DP来搞.只是拿费用流还是非常方便的. 能够转化成求满流为2 的最小花费.一般做法为拆点,对于 i 拆为2*i 和 2*i+1.然后连 ...
- hdu 1533 KM或费用流
以前用KM写过,现在再用费用流写. #include <iostream> #include <cstdio> #include <cstring> #includ ...
- HDU 2686 Matrix(最大费用流)
Matrix Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) Total Sub ...
- HDU 3315 My Brute(二分图最佳匹配+尽量保持原先匹配)
http://acm.hdu.edu.cn/showproblem.php?pid=3315 题意: 有S1到Sn这n个勇士要和X1到Xn这n个勇士决斗,初始时,Si的决斗对象是Xi. 如果Si赢了X ...
- Cyclic Tour HDUOJ 费用流
Cyclic Tour Time Limit: 1000/1000 MS (Java/Others) Memory Limit: 32768/65535 K (Java/Others)Total ...
- hdu 2686 Matrix 最小费用最大流
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=2686 Yifenfei very like play a number game in the n*n ...
- hdu3315-My Brute(费用流 or KM算法)
题目:My Brute Seaco是一个漂亮的妹子,喜欢玩一款名叫My Brute的游戏.情人节快到了,starvae和xingxing都想邀请妹子过节,但是妹子只能陪一个啊,于是两个人决定打一架,用 ...
随机推荐
- 研究一下TForm.WMPaint过程(也得研究WM_ERASEBKGND)——TForm虽然继承自TWinControl,但是自行模仿了TCustomControl的全部行为,一共三种自绘的覆盖方法,比TCustomControl还多一种
先擦除背景: procedure TCustomForm.WMEraseBkgnd(var Message: TWMEraseBkgnd); begin if not IsIconic(Handle) ...
- js获取当前日期,网页头部用
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN" "http://www.w3.org/ ...
- ZigBee研究之旅(一)
*********************************************************************** 以下有引用webee公司的文档的内容,版权属于webee ...
- alv 显示 汇总数据
alv中有 字段设置 gs_fieldcat-do_sum = ‘1’.字段 用于控制汇总
- EFI/GPT探索(为何win7分区时创建100M隐藏分区)
EFI/GPT探索(为何win7分区时创建100M隐藏分区) 转自 http://blog.tomatoit.net/article.asp?id=348 EFI/GPT是新一代的固件/启动管理技术, ...
- java.lang.NoSuchMethodError: org.springframework.beans.factory.annotation.InjectionMetadata.<init>(L
关于错误: java.lang.NoSuchMethodError: org.springframework.beans.factory.annotation.InjectionMetadata.&l ...
- Qt Windows下开机自动启动自己的程序
源地址:http://blog.csdn.net/chrisfxs/article/details/13279491 版权声明:本文为博主原创文章,未经博主允许不得转载. void Widget::R ...
- [Python] heapq简介
[Python] heapq简介 « Lonely Coder [Python] heapq简介 judezhan 发布于 2012 年 8 月 8 日 暂无评论 发表评论 假设你需要维护一个列表,这 ...
- Android实现位图剪切
我们不能总是依赖于BitmapFactory 以下告诉大家怎么从Bitmaqp中截取某一部分创建新的Bitmap 系统会有一个默认png图片:icon.png 可是这个图片中最外层会有白色的 比較讨 ...
- tomcat各版本和jsp、jstl、servlet的依赖关系(转)
Servlet / JSP / Tomcat Version Servlet/ JSP Tomcat 2.5/2.1 6.0.18 2.4/2.0 5.5.27 2.3/1.2 4.1.3 ...