Codeforces 884F Anti-Palindromizes [费用流]

Descripiton:
定义反回文串为每个位置和对称的位置字符都不相同。现在给一个偶数长度的串，可以随便排列字符，给出每位一个价值，定义总价值为改变后的字符串与原先字符串对应相等位置的价值的和，问最大价值。

---

Solution:
巧妙的费用流模型。
建图有5层：源，字符，\frac_{n}{2}对位置，n个位置，汇。
源向字符连边，容量为出现次数，费用为0。
字符向对应位置对连边，具体讨论见代码。
每对位置向对应位置连边。
位置向汇连边，容量为2。
这样成功满足了反回文串的条件。

---

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
#include <cmath>
using namespace std;
const int N = 2005, inf = 0x3f3f3f3f;
struct edge {
    int nxt, to, f, w;
} e[2000005];
int n, m, source, sink, cnt = 1;
int h[N], d[N], Prev[N], pree[N], vis[N], b[N];
char s[N];
void link(int u, int v, int f, int c) {
    e[++cnt].nxt = h[u];
    h[u] = cnt;
    e[cnt].to = v;
    e[cnt].w = c;
    e[cnt].f = f;
}
void insert(int u, int v, int f, int c) {
    link(u, v, f, c);
    link(v, u, 0, -c);
}
bool spfa() {
    queue<int> q;
    memset(d, 0x3f3f, sizeof(d));
    d[source] = 0;
    q.push(source);
    while(!q.empty()) {
        int u = q.front();
        q.pop();
        vis[u] = 0;
        for(int i = h[u]; i; i = e[i].nxt) {
            if(d[e[i].to] > d[u] + e[i].w && e[i].f) {
                d[e[i].to] = d[u] + e[i].w;
                Prev[e[i].to] = u;
                pree[e[i].to] = i;
                if(!vis[e[i].to]) {
                    vis[e[i].to] = 1;
                    q.push(e[i].to);
                }
            }
        }
    }
    return d[sink] != inf;
}
int mcf() {
    int ret = 0;
    while(spfa()) {
        int now = sink, delta = inf;
        while(now != source) {
            delta = min(delta, e[pree[now]].f);
            now = Prev[now];
        }
        now = sink;
        while(now != source) {
            e[pree[now]].f -= delta;
            e[pree[now] ^ 1].f += delta;
            now = Prev[now];
        }
        ret += delta * d[sink];
    }
    return ret;
}
int main() {
    scanf("%d%s", &n, s + 1);
    m = n / 2;
    sink = 27 + m + n;
    for(int i = 1; i <= n; ++i) {
        insert(source, s[i] - 'a' + 1, 1, 0);
    }
    for(int i = 1; i <= n; ++i) {
        scanf("%d", &b[i]);
    }
    for(int j = 1; j <= 26; ++j) {
        for(int i = 1; i <= m; ++i) {
            if(s[i] - 'a' + 1 == j && s[n - i + 1] - 'a' + 1 == j) {
                insert(j, i + 26, 1, -max(b[i], b[n - i + 1]));
            } else if(s[i] - 'a' + 1 == j) {
                insert(j, i + 26, 1, -b[i]);
            } else if(s[n - i + 1] - 'a' + 1 == j) {
                insert(j, i + 26, 1, -b[n - i + 1]);
            } else {
                insert(j, i + 26, 1, 0);
            }
        }
    }       
    for(int i = 1; i <= m; ++i) {   
        insert(i + 26, i + 26 + m, 1, 0);
        insert(i + 26, n - i + 27 + m, 1, 0);
    }
    for(int i = 1; i <= n; ++i) {
        insert(i + m + 26, sink, 2, 0);
    }
    printf("%d\n", -mcf());
    return 0;
}