HDU 5772 String problem（最大权闭合子图）

题意：

      有一个长度为n(n<=100)的串，字符集为'0'-'9'。求分数最大。

      分数定义： Score=Value - Cost

         Value=0;
        for(int i=1;i<=length(substr);++i){
           for(int j=1;j<=length(substr);++j){
             if(i!=j)
                  Value+=w[id[i]][id[j]];
              }
        }

对于字符x：

cost[x]=ax∗(kx−1)+bx,kx≠0

（kx为x在子序列出现的次数，ax，bx输入会给出）

Cost所有cost[x]相加

解析：   一眼不会0。0   先搬官方题解。

首先将点分为3类

第一类：Pij 表示第i个点和第j个点组合的点，那么Pij的权值等于w[i][j]+w[j][i]（表示得分）

第二类：原串中的n个点每个点拆出一个点，第i个点权值为 –a[s[i]] （表示要花费）

第三类：对于10种字符拆出10个点，每个点的权值为  -(b[x]-a[x])

 

那么我们可以得到一个关系图 ，对于第一类中的点Pij，如果想要选择Pij，你就必须要选中第二类中的点i和j，对于第二类中的点如果你想选中第i个点，其对应的字符s[i]，那么就必须选中第三类中s[i] 对应的点，因为每个种类的点第一次选中时花费是b[s[i]]，而第二类中花费都是a[s[i]],一定要补上b[s[i]]-a[s[i]]，而且只需要补上一次。

得到上面的关系图后然后就是普通的最大权闭合子图问题,直接求解即可。

先了解下最大权闭合子图：

所以再把最大权闭合子图转换下，用最大流求解：

对于点权为正的节点，从S连一条容量为点权的边到该点，对于点权为负的边，从该点连一条容量为点权绝对值的边到T。原图中的边保留，容量为inf。

答案 = 所有正点权之和 - 最大流（最小割）。

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<string>
#include<vector>
#include<cctype>
#include<set>
#include<map>
#include<queue>
#include<stack>
#include<iomanip>
#include<sstream>
#include<limits>
#define ll long long
#define inf 0x3f3f3f3f
using namespace std;
const int maxn = 2e5+10;
const ll MOD = 1000000007;
const double EPS = 1e-10;
const double Pi = acos(-1.0);
#define INF 0x7fffffff

struct Edge{
    int to , cap,rev;  // rev反向边
};
vector<Edge>G[maxn];
int level[maxn],iter[maxn];  // 顶点到源点的标号   当前弧
void add_edge(int from,int to,int cap) //cap容量
{
    G[from].push_back((Edge){to,cap,G[to].size()});
    G[to].push_back((Edge){from,0,G[from].size()-1 });
}
void bfs(int s) // 分层
{
    memset(level,-1,sizeof(level));
    queue<int>que;
    level[s] = 0;
    que.push(s);
    while(!que.empty())
    {
        int v = que.front();que.pop();
        for(int i = 0; i < G[v].size(); i++)
        {
            Edge &e = G[v][i];
            if (e.cap>0 && level[e.to] < 0)
            {
                level[e.to] = level[v]+1;
                que.push(e.to);
            }
        }
    }
}
int dfs(int v, int t, int f) //寻增广路
{
    if (v==t) return f;
    for(int &i = iter[v]; i < G[v].size(); i++)
    {
        Edge &e = G[v][i];
        if (e.cap > 0 && level[v] < level[e.to])
        {
            int d = dfs(e.to , t, min(f,e.cap));
            if (d > 0)
            {
                e.cap -= d;
                G[e.to][e.rev].cap += d;
                return d;
            }
        }
    }
    return 0;
}
int max_flow(int s, int t) // s到t最大流
{
    int flow = 0;
    for(;;)
    {
        bfs(s);
        if (level[t] < 0) return flow;
        memset(iter,0,sizeof(iter));
        int f;
        while((f = dfs(s,t,INF)) > 0) { flow += f;}
    }
}
int s[200];
char str[200];
int a[15],b[15],w[200][200];
int main(){
#ifdef LOCAL
	freopen("C:\\Users\\lanjiaming\\Desktop\\acm\\in.txt","r",stdin);
	//freopen("output.txt","w",stdout);
#endif
//ios_base::sync_with_stdio(0);
    int T,kase = 0;
    scanf("%d",&T);
    while(T--)
    {
        int n;
        scanf("%d",&n);
        scanf(" %s",str+1);
        for(int i = 1; i <= n; i++) s[i] =str[i]- '0';
        for(int i = 0; i < 10; i++) scanf("%d%d",&a[i],&b[i]);
        int ans = 0;
         for(int i = 1; i <= n;i++)
            for(int j = 1; j <= n;j++)
            {
                scanf("%d",&w[i][j]);
                ans += w[i][j];
            }
        int S = 0,T = n*(n-1)/2+n+11,cnt = 0;
        for(int i= 0; i <= T; i++) G[i].clear();
        for(int i = 1; i <= n; i++)         //第一类点
            for(int j = i+1; j <= n; j++)
            {
                cnt++;
                add_edge(S,cnt,w[i][j]+w[j][i]);
                add_edge(cnt,n*(n-1)/2+i,inf);
                add_edge(cnt,n*(n-1)/2+j,inf);
            }
        for(int i = 1; i <= n; i++)      //第二类点
        {
            add_edge(n*(n-1)/2+i,n*(n-1)/2+n+s[i]+1,inf);
            add_edge(n*(n-1)/2+i,T,a[s[i]]);
        }
        for(int i = 0; i < 10; i++)    //第三类点
            add_edge(n*(n-1)/2+n+i+1,T,b[i]-a[i]);

        ans -= max_flow(S,T);
        printf("Case #%d: %d\n",++kase,ans);
    }
    return 0;
}