hdu 3395 Special Fish 最大费用流

28 篇文章 0 订阅

题意:n条鱼可以相互攻击,得到的值为这两条鱼的异或值。每条鱼只能被攻击一次。求能得到的亦或者的最大值。

很容易想到最大费用流。拆点,对每个 i ,连边(s, i, 1, 0), (i + n, t, 1, 0)。当 ij = 1 时, 对 i —> j + n 连边(i, j + n, 1, -(val[i] ^ val[j]))。

由于题目的特殊性,当最大费用出现时, 并不一定是最大流, 而最大费用流的前提是最大流,所以再建立一些额外的边,保证流量不破坏费用的前提下保证满流。比较方便的是直接建边(i, t, 1, 0)。

#include <bits/stdc++.h>
#include <map>
#include <set>
#include <queue>
#include <stack>
#include <cmath>
#include <time.h>
#include <vector>
#include <cstdio>
#include <string>
#include <iomanip>
///cout << fixed << setprecision(13) << (double) x << endl;
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <algorithm>

using namespace std;

#define lson l, mid, rt << 1
#define rson mid + 1, r, rt << 1 | 1
#define ls rt << 1
#define rs rt << 1 | 1
#define pi acos(-1.0)
#define eps 1e-8
#define Mp(a, b) make_pair(a, b)
#define asd puts("asdasdasdasdasdf");
typedef long long ll;
typedef pair <int, int> pl;
//typedef __int64 LL;
const int inf = 0x3f3f3f3f;
const int N = 440;
const int st = N-9, ed = N-8;

struct node{
    int v, w, c, nxt;   //c--->cap  w--->cost(dis)
}e[N*N*2];

int cnt, n, m;
int head[N];
int dis[N]; //费用看成距离
int cap[N];     //流量
int cur[N];     //当前弧
int vis[N];
queue <int> q;
char mat[N][N];
int flow, cost;
vector <pl> H, M;

void init()
{
    cnt = 0;
    memset( head, -1, sizeof( head ) );
    memset( mat, inf, sizeof( mat ) );
    H.clear(), M.clear();
}

void add( int u, int v, int c, int w )
{
    e[cnt].v = v;
    e[cnt].c = c;
    e[cnt].w = w;
    e[cnt].nxt = head[u];
    head[u] = cnt++;

    e[cnt].v = u;
    e[cnt].c = 0;
    e[cnt].w = -w;
    e[cnt].nxt = head[v];
    head[v] = cnt++;
}

bool spfa()
{
    memset( dis, inf, sizeof( dis ) );
    memset( vis, 0, sizeof( vis ) );
    while( !q.empty() ) q.pop();
    cap[st] = inf;
    cur[st] = -1;
    dis[st] = 0;
    q.push( st );
    while( !q.empty() ) {
        int u = q.front();
        q.pop();
        vis[u] = 0;
        for( int i = head[u]; ~i; i = e[i].nxt ) {
            int v = e[i].v, c = e[i].c, w = e[i].w;
            if( c && dis[v] > dis[u] + w ) {
                dis[v] = dis[u] + w;
                cap[v] = min( c, cap[u] );
                cur[v] = i;
                if( !vis[v] ) {
                    vis[v] = 1;
                    q.push(v);
                }
            }
        }
    }
    if( dis[ed] == inf )
        return 0;
    if( dis[ed] >= 0 )  //费用开始减少时停止
        return 0;
    flow += cap[ed];
    cost += cap[ed] * dis[ed];
    for( int i = cur[ed]; ~i; i = cur[e[i^1].v] ) {
        e[i].c -= cap[ed];
        e[i^1].c += cap[ed];
    }
    return 1;
}

int MCMF()
{
    flow = cost = 0;
    while( spfa() );
    return cost;
}

int main()
{
    while( ~scanf("%d", &n ) && n ) {
        init();
        int val[N];
        for( int i = 1; i <= n; ++i ) {
            scanf("%d", &val[i]);
        }
        char ch[132];
        for( int i = 1; i <= n; ++i ) {
            scanf("%s", ch+1);
            for( int j = 1; j <= n; ++j ) {
                if( ch[j] == '1' )
                    add( i, j+n, 1, -(val[i]^val[j]));
            }
        }
        for( int i = 1; i <= n; ++i ) {
            add( st, i, 1, 0 );
            add( i+n, ed, 1, 0 );
            //由于题目的特殊性,当最大费用出现时, 并不一定是最大流,
            //而最大费用流的前提是最大流,所以再建立一些额外的边,
            //保证流量不破坏费用的前提下保证满流。比较方便的是直接建边(i, t, 1, 0)。
            //add( i, ed, 1, 0 );
        }
        int ans = MCMF();
        printf("%d\n", -ans);
    }
    return 0;
}
</span>

见:http://www.cnblogs.com/ac-luna/p/3757446.html

  • 1
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值