HDU 3395 Special Fish | 费用流

*本题的输入有点坑，直接看输入样例，应该是没有其他字符的。如果我不是在测数据的时候多按了下回车，也不知道自己栽在输入上。后来根据别人的代码改了下输入方式就AC了（之前我是逐个输入的，%c）。

题意：

一些奇怪的鱼，都认为自己雄性，它认为别的鱼是雌性就会Attack。

每条鱼最多可以被Attack一次以及Attack别的鱼一次。注意，是最多。因此每条鱼都可以不进行Attack。

每条鱼都有一个权值，Attack后可以得到一个值，是由这两条鱼的权值异或所得，记为sum

求出sum的最大值。

思路：

主要是题意理解的问题，一开始我也理解错了。

建图：

一条鱼可以攻击和被攻击最多一次。拆点

每条鱼都可以不攻击别的鱼，每条鱼再连一条边到汇点，容量1，费用0

AC代码：

<span style="font-family:Microsoft YaHei;">#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <queue>
using namespace std;

#define debug cout<<"??"<<endl;
const int MAXN = 205;
const int MAXEDGE = 1e5 + 5;
const int INF = 0x3f3f3f3f;
int n, a[105];
struct CEdge
{
        int from, to, cap, flow, cost, next;
}edge[MAXEDGE];
struct CMCMF
{
        int s, t, pp;
        int head[MAXN], a[MAXN], d[MAXN], p[MAXN];
        bool inq[MAXN];
        CMCMF(int ss, int tt)
        {
                s = ss, t = tt;
                pp = 0;
                memset(head, -1, sizeof(head));
        }
        inline void addEdge(int u, int v, int cap, int cost)
        {
                edge[pp] = (CEdge){u, v, cap, 0, cost, head[u]};
                head[u] = pp++;
                edge[pp] = (CEdge){v, u, 0, 0, -cost, head[v]};
                head[v] = pp++;
        }
        bool bellmanFord(int &flow, int &cost)
        {
                memset(d, INF, sizeof(d));
                memset(inq, false, sizeof(inq));
                queue <int> q;
                q.push(s), inq[s] = true;
                a[s] = INF, d[s] = 0, p[s] = -1;
                while(!q.empty())
                {
                        int u = q.front();
                        q.pop(), inq[u] = false;
                        int next = head[u];
                        while(next != -1)
                        {
                                CEdge &e = edge[next];
                                if(e.cap > e.flow && d[e.to] > d[u] + e.cost)
                                {
                                        d[e.to] = d[u] + e.cost;
                                        a[e.to] = min(a[u], e.cap - e.flow);
                                        p[e.to] = next;
                                        if(!inq[e.to])  inq[e.to] = true, q.push(e.to);
                                }
                                next = e.next;
                        }
                }
                if(d[t] == INF) return false;
                cost += d[t] * a[t];
                flow += a[t];
                int u = t;
                while(u != s)
                {
                        edge[p[u]].flow += a[t];
                        edge[p[u]^1].flow -= a[t];
                        u = edge[p[u]].from;
                }
                return true;
        }
};
int main()
{
        while(scanf("%d", &n), n)
        {
                int s = 0, t = 2*n+1;
                CMCMF mcmf(s, t);
                for(int i = 1;i <= n; i++)
                {
                        scanf("%d",&a[i]);
                }
                getchar();
                //s, t
                for(int i = 1;i <= n; i++)
                {
                        mcmf.addEdge(s, i, 1, 0);
                        mcmf.addEdge(i+n, t, 1, 0);
                        mcmf.addEdge(i, t, 1, 0);
                }
                //others

                int tmp;
                for(int i = 1;i <= n; i++)
                {
                        for(int j = 1;j <= n; j++)
                        {
                                scanf("%1d", &tmp);
                                if(tmp)
                                        mcmf.addEdge(i, j+n, 1, -(a[i]^a[j]));
                        }
                        getchar();
                }
                int cost = 0, flow = 0;
                while(mcmf.bellmanFord(flow, cost));
                cout<<-cost<<endl;
        }
        return 0;
}</span>