
3 4 5 6 6 8 10 12 14 16 7 7 6 7 3 5 3 4 5 6 6 8 10 12 14 16 5 5 5 5 5 5 0
7 33.333% Oh, I lose my dear seaco!
输出最大得分和比例。
假设分数<=0则输出Oh, I lose my dear seaco!
#include<stdio.h>
#include<string.h>
#include<queue>
using namespace std;
const int MAXN = 1010;
const int MAXM = 100100;
const int INF = 1<<29;
struct EDG{
int to,next,cap;
int cost,flag;
}edg[MAXM];
int head[MAXN],eid;
int pre[MAXN], cost[MAXN] ; //点0~(n-1)
void init(){
eid=0;
memset(head,-1,sizeof(head));
}
void addEdg(int u,int v,int cap,int cst,int flag){
edg[eid].to=v; edg[eid].next=head[u]; edg[eid].cost = cst;
edg[eid].cap=cap; edg[eid].flag=flag; head[u]=eid++;
edg[eid].to=u; edg[eid].next=head[v]; edg[eid].cost = -cst;
edg[eid].cap=0; edg[eid].flag=-flag; head[v]=eid++;
}
bool inq[MAXN];
int q[MAXN],flg[MAXN];
bool spfa(int sNode,int eNode,int n){
int l=0 , r=0;
for(int i=0; i<n; i++){
inq[i]=false; cost[i]= -INF; flg[i]=-INF;
}
cost[sNode]=0; flg[sNode]=0; inq[sNode]=1; pre[sNode]=-1;
q[r++]=sNode;
while(l!=r){
int u=q[l++];
if(l==MAXN)l=0;
inq[u]=0;
for(int i=head[u]; i!=-1; i=edg[i].next){
int v=edg[i].to;
if(edg[i].cap<=0)continue;
if( cost[v]<cost[u]+edg[i].cost || cost[v]==cost[u]+edg[i].cost&&flg[v]<flg[u]+edg[i].flag){ //在满足可增流的情况下,最小花费
cost[v] = cost[u]+edg[i].cost;
flg[v] = flg[u]+edg[i].flag;
pre[v]=i; //记录路径上的边
if(!inq[v]){
if(r==MAXN)r=0;
q[r++]=v;
inq[v]=1;
}
}
}
}
return cost[eNode]!=-INF; //推断有没有增广路
}
//反回的是最大流,最小花费为minCost
int minCost_maxFlow(int sNode,int eNode ,int& minCost,int n){
int ans=0;
while(spfa(sNode,eNode,n)){
ans+=flg[eNode];
minCost+= cost[eNode];
for(int i=pre[eNode]; i!=-1; i=pre[edg[i^1].to]){
edg[i].cap-=1; edg[i^1].cap+=1;
}
}
return ans;
}
int main(){
//输入,初始化init()
int n,valu[MAXN],H[MAXN],P[MAXN],A[MAXN],B[MAXN];
while(scanf("%d",&n)>0&&n){
for(int i=1; i<=n; i++)
scanf("%d",&valu[i]);
for(int i=1; i<=n; i++)
scanf("%d",&H[i]);
for(int i=1; i<=n; i++)
scanf("%d",&P[i]);
for(int i=1; i<=n; i++)
scanf("%d",&A[i]);
for(int i=1; i<=n; i++)
scanf("%d",&B[i]);
init();
int s=0 , t= 2*n+1;
for(int i=1; i<=n; i++)
{
addEdg(s , i , 1 , 0 , 0);
addEdg(i+n, t , 1 , 0 , 0);
}
for(int i=1; i<=n; i++)
for(int j=1; j<=n; j++)
if(P[j]/A[i]+(P[j]%A[i]!=0?1:0)<=H[i]/B[j]+(H[i]%B[j]!=0?
1:0)){
if(i==j)
addEdg(i,j+n,1,valu[i],1);
else
addEdg(i,j+n,1,valu[i],0);
}
else{
if(i==j)
addEdg(i,j+n,1,-valu[i],1);
else
addEdg(i,j+n,1,-valu[i],0);
}
int maxcost=0;
int ans=minCost_maxFlow(s , t , maxcost, t+1);
if(maxcost<=0)
printf("Oh, I lose my dear seaco!\n");
else
printf("%d %.3f%%\n",maxcost,100.0*ans/n);
}
}
原文:http://www.cnblogs.com/mthoutai/p/6748809.html