码迷,mamicode.com
首页 > 其他好文 > 详细

HDU 3315 My Brute(费用流)

时间:2015-07-14 13:14:03      阅读:114      评论:0      收藏:0      [点我收藏+]

标签:

职务地址: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(费用流)

标签:

原文地址:http://www.cnblogs.com/lcchuguo/p/4645129.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!