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>