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

Description
给出一个数字串，现在要选出一个价值-代价最大的子序列，子序列价值和代价计算方法如下：
子序列的价值: sum(w[id(i)][id(j)]) (i !=j) id(i)为这个子序列中第i个字符在原串中的位置
子序列的代价: 设x为数字i出现的次数，则代价为a[i]*(x-1)+b[i] (x>0) 0 (x=0)
Input
第一行为一整数T表示用例组数，每组用例首先输入数字串的长度n，然后是一个长度为n的数字串，之后10行每行两个整数ai和bi，最后一个n*n矩阵表示w
(T<=20,n<=100,0<=ai<=bi<=1000,0<=w[i][j]<=50)
Output
对于每组用例，输出最大化的价值-代价值
Sample Input
1
3
135
1 2
1 2
1 2
1 2
1 2
1 2
1 2
1 2
1 2
1 2
0 0 3
1 0 0
4 0 0
Sample Output
Case #1: 3
Solution
考虑三种点
1.Pij(1<=i< j<=n)，表示子序列中有第i个字符和第j个字符，这个点的点权是w[i][[j]+w[j][i]
2.i(1<=i<=n)，表示子序列中有i字符x，这个点的点权是-ax
3.x(0<=x<=9)，表示子序列中有x这个数字，这个点的点权是-(bx-ax)
选Pij就要选i点和j点，选i点就要选s[i]这个数字，选j点就要选s[j]这个数字
问题变成从上面三种点中选择若干点使得所选点满足上述关系且总点权最大，问题转化成求最大权闭合子图问题
Code

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
typedef long long ll;
#define maxn 11111 
#define maxm 111111
#define INF 0x3f3f3f3f
int head[maxn],cur[maxn],d[maxn],st[maxm],s,e,no;//s为源点,e为汇点,n为点数,no为边数 
struct point
{
    int u,v,flow,next;
    point(){};
    point(int x,int y,int z,int w):u(x),v(y),next(z),flow(w){};
}p[maxm];
void add(int x,int y,int z)//从x到y建容量为z的边 
{
    p[no]=point(x,y,head[x],z);//前向弧,标号为偶 
    head[x]=no++;
    p[no]=point(y,x,head[y],0);//后向弧,标号为奇 
    head[y]=no++;
}
void init()//初始化 
{
    memset(head,-1,sizeof(head));
    no=0;
}
bool bfs()
{
    int i,x,y;
    queue<int>q;
    memset(d,-1,sizeof(d));
    d[s]=0; 
    q.push(s);
    while(!q.empty())
    {
        x=q.front();    
        q.pop();
        for(i=head[x];i!=-1;i=p[i].next)
        {
            if(p[i].flow&& d[y=p[i].v]<0)
            {
                d[y]=d[x]+1;
                if(y==e)    
                    return true;
                q.push(y);
            }
        }
    }
    return false;
}
int dinic()//最大流 
{
    int i,loc,top,x=s,nowflow,maxflow=0;
    while(bfs())
    {
        memcpy(cur,head,sizeof(head));
        top=0;
        while(true)
        {
            if(x==e)
            {
                nowflow=INF;
                for(i=0;i<top;i++)
                {
                    if(nowflow>p[st[i]].flow)
                    {
                        nowflow=p[st[i]].flow;
                        loc=i;
                    }
                }
                for(i=0;i<top;i++)
                {
                    p[st[i]].flow-=nowflow;
                    p[st[i]^1].flow+=nowflow;
                }
                maxflow+=nowflow;
                top=loc;    
                x=p[st[top]].u;
            }
            for(i=cur[x];i!=-1;i=p[i].next)
                if(p[i].flow&&d[p[i].v]==d[x]+1) 
                    break;
            cur[x]=i;
            if(i!=-1)
            {
                st[top++]=i;
                x=p[i].v;
            }
            else 
            {
                if(!top)    
                    break;
                d[x]=-1;
                x=p[st[--top]].u;
            }
        }
    }
    return maxflow;
}
int T,n,a[11],b[11],w[111][111],Case=1;
char c[111];
int main()
{
    scanf("%d",&T);
    while(T--)
    {
        init();
        scanf("%d%s",&n,c+1);
        for(int i=0;i<10;i++)scanf("%d%d",&a[i],&b[i]);
        int sum=0;
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)
            {
                scanf("%d",&w[i][j]);
                sum+=w[i][j];
            }
        s=0,e=n*n+n+11;
        for(int i=1;i<=n;i++)
            for(int j=i+1;j<=n;j++)
            {
                int t=i*n+j;
                add(s,t,w[i][j]+w[j][i]);
                add(t,n*n+i,INF);
                add(t,n*n+j,INF);
            }
        for(int i=1;i<=n;i++)
        {
            add(n*n+i,e,a[c[i]-'0']);
            add(n*n+i,n*n+n+c[i]-'0'+1,INF);
        }
        for(int i=0;i<10;i++)
            add(n*n+n+i+1,e,b[i]-a[i]);
        int ans=dinic();
        printf("Case #%d: %d\n",Case++,sum-ans);
    }
}