题意: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