My Brute HDU - 3315（KM || 费用流）

题意：

有S1到Sn这n个勇士要和X1到Xn这n个勇士决斗,初始时,Si的决斗对象是Xi. 如果Si赢了Xi,那么你将获得Vi分,否则你将获得-Vi分. Si和Xi对决时,Si有初始生命Hi,初始攻击Ai, Xi有初始生命Pi,初始攻击Bi. 且Si先出手,然后Xi失去Ai生命,之后如果Xi没死,那么Xi出手,Si失去Bi生命. 直到有一方的生命值<=0时,决斗结束.

现在要你重新安排S和X的决斗顺序,使得你能获得的分最多.如果有多个最优解,你要选取那个维持初始决斗顺序最多的解.

解析：

　　仔细看题是不是攻击力和血量除了确定连接关系外并无用处又因为边有权值且每个只能用一次那么自然而然想到 KM 或费用流

　　这里是用费用流做的

　　费用流是求最小这里求最大所以取反即可

　　那怎么求出没有改变顺序的勇士

　　想一下求最小割集的做法

　　如果 i == j 建立费用为 -w * (n - 1) - 1的边 - 1 表示没有改变顺序

　　否则建立费用为 w * (n - 1) 的边

　　容量都为1

　　那么最后求出的最小费用取反之后是不是就是最大的分数设为value

　　那么value / （n + 1）是不是就是实际的分数 value % （n + 1）是不是就是没有改变顺序的勇士仔细想一想

　　还有%一定要 %% 不然超时。。。emm。。。

#include <iostream>

#include <cstdio>

#include <sstream>

#include <cstring>

#include <map>

#include <cctype>

#include <set>

#include <vector>

#include <stack>

#include <queue>

#include <algorithm>

#include <cmath>

#include <bitset>

#define rap(i, a, n) for(int i=a; i<=n; i++)

#define rep(i, a, n) for(int i=a; i<n; i++)

#define lap(i, a, n) for(int i=n; i>=a; i--)

#define lep(i, a, n) for(int i=n; i>a; i--)

#define rd(a) scanf("%d", &a)

#define rlld(a) scanf("%lld", &a)

#define rc(a) scanf("%c", &a)

#define rs(a) scanf("%s", a)

#define pd(a) printf("%d\n", a);

#define plld(a) printf("%lld\n", a);

#define pc(a) printf("%c\n", a);

#define ps(a) printf("%s\n", a);

#define MOD 2018

#define LL long long

#define ULL unsigned long long

#define Pair pair<int, int>

#define mem(a, b) memset(a, b, sizeof(a))

#define _  ios_base::sync_with_stdio(0),cin.tie(0)

//freopen("1.txt", "r", stdin);

using namespace std;

const int maxn =  + , INF = 0x7fffffff, LL_INF = 0x7fffffffffffffff;

int n, m, s, t, value, flow, cnt;

int V[maxn], H[maxn], P[maxn], A[maxn], B[maxn];

int head[maxn], d[maxn], p[maxn], f[maxn], vis[maxn];

struct node

{

    int u, v, w, c, next;

}Node[];

void add_(int u, int v, int w, int c)

{

    Node[cnt].u = u;

    Node[cnt].v = v;

    Node[cnt].w = w;

    Node[cnt].c = c;

    Node[cnt].next = head[u];

    head[u] = cnt++;

}

void add(int u, int v, int w, int c)

{

    add_(u, v, w, c);

    add_(v, u, -w, );

}

int spfa()

{

    queue<int> Q;

    mem(vis, );

    mem(p, -);

    for(int i = ; i < maxn; i++) d[i] = INF;

    d[s] = ;

    Q.push(s);

    vis[s] = ;

    p[s] = ; f[s] = INF;

    while(!Q.empty())

    {

        int u = Q.front(); Q.pop();

        vis[u] = ;

        for(int i = head[u]; i != -; i = Node[i].next)

        {

            node e = Node[i];

            if(d[e.v] > d[u] + e.w && e.c > )

            {

                d[e.v] = d[u] + e.w;

                p[e.v] = i;

                f[e.v] = min(f[u], e.c);

                if(!vis[e.v])

                {

                    vis[e.v] = ;

                    Q.push(e.v);

                }

            }

        }

    }

    if(p[t] == -) return ;

    flow += f[t]; value += f[t] * d[t];

    for(int i = t; i != s; i = Node[p[i]].u)

    {

        Node[p[i]].c -= f[t];

        Node[p[i] ^ ].c += f[t];

    }

    return ;

}

void max_flow()

{

    flow = value = ;

    while(spfa());

    if(-value <= )

        printf("Oh, I lose my dear seaco!\n");

    else

        printf("%d %.3f%%\n", -value / (n + ), -value % (n + ) / (double) n * );

}

int main()

{

    while(scanf("%d", &n) != EOF && n)

    {

        mem(head, -);

        cnt = ;

        s = , t = n *  + ;

        for(int i = ; i <= n; i++)

        {

            rd(V[i]);

            add(s, i, , );

            add(i + n, t, , );

        }

        for(int i = ; i <= n; i++)

            rd(H[i]);

        for(int i = ; i <= n; i++)

            rd(P[i]);

        for(int i = ; i <= n; i++)

            rd(A[i]);

        for(int i = ; i <= n; i++)

            rd(B[i]);

        for(int i = ; i <= n; i++)

        {

            for(int j = ; j <= n; j++)

            {

                if(ceil(P[j] / (double) A[i]) <= ceil(H[i] / (double) B[j]))

                {

                    add(i, n + j, -V[i] * (n + ) - ((i == j) ?  : ), );

                }

                else

                {

                    add(i, n + j, V[i] * (n + ) - ((i == j) ?  : ), );

                }

            }

        }

        max_flow();

    }

    return ;

}

巴特西

My Brute HDU - 3315（KM || 费用流）

最新文章

热门文章