hdu 5772 String problem（2016多校第四场1009）最大权闭合图

每个w[i][j]+w[j][i]存在依赖于字母i,j的存在，i,j的存在依赖于0~9对应数字的存在，明显的最大权闭合图。

#include<stdio.h>
#include<iostream>
#include<algorithm>
#include<vector>
#include<string.h>
#include<queue>
#define maxn 10005
#define inf 100000000
using namespace std;
struct pi
{
    int to;
    int cost;
    int rev;
}pp;
vector<pi >g[maxn];
queue<int>q;
int s,t;
int line[maxn],leve[maxn];
int min(int a,int b)
{
    int p;
    p=a;
    if(b<a)
        p=b;
    return p;
}
void add(int a,int b,int cos)
{
    pp.to=b;
    pp.cost=cos;
    pp.rev=(int)g[b].size();
    g[a].push_back(pp);
    pp.to=a;
    pp.cost=0;
    pp.rev=(int)g[a].size()-1;
    g[b].push_back(pp);
    return ;
}
void bfs(void)
{
    q.push(s);
    int p,f,i;
    while(!q.empty())
    {
        p=q.front();
        q.pop();
        f=(int)g[p].size();
        for(i=0;i<f;i++)
        {
            pi &e=g[p][i];
            if(e.cost>0&&line[e.to]<0)
            {
                line[e.to]=line[p]+1;
                q.push(e.to);
            }
        }
    }
    return ;
}
int dfs(int v,int f)
{
    int p;
    if(v==t)
        return f;
    for(int &i=leve[v];i<g[v].size();i++)
    {
        pi &e=g[v][i];
        if(line[v]<line[e.to]&&e.cost>0)
        {
            p=dfs(e.to,min(f,e.cost));
            if(p>0)
            {
                e.cost-=p;
                g[e.to][e.rev].cost+=p;
                return p;
            }
        }
    }
    return 0;
}
int  dinic()
{
    int f;
    int flow=0;
    while(1)
    {
        memset(line,-1,sizeof(line));
        memset(leve,0,sizeof(leve));
        line[s]=0;
        bfs();
        if(line[t]<0)
            return flow;
        f=dfs(s,inf);
        if(f==0)
            continue;
        flow+=f;
        while((f=dfs(s,inf))>0)
        {
            flow+=f;
        }
    }
}
int a[10];
int b[10];
int w[105][105];
char c[105];
int main()
{
    int t1,n,N=0;
    cin>>t1;
    while (t1--) {
        scanf("%d",&n);
        scanf("%s",c);
        for(int i=n;i>=1;i--) c[i]=c[i-1];
        for(int i=0;i<10;i++){
            scanf("%d%d",&a[i],&b[i]);
        }
        for(int i=1;i<=n;i++){
            for(int j=1;j<=n;j++){
                scanf("%d",&w[i][j]);
            }
        }
        s=0;
        t=10+n+(n*(n-1))/2+1;
        int all=0;
        for(int i=0;i<=t;i++) g[i].clear();
        for(int i=1;i<=10;i++){
            add(i,t,b[i-1]-a[i-1]);
        }
        for(int i=1;i<=n;i++){
            add(10+i,c[i]-'0'+1, 1000000000);
            add(10+i,t,a[c[i]-'0']);
        }
        int tot=10+n;
        for(int i=1;i<=n;i++){
            for(int j=i+1;j<=n;j++){
                tot++;
                add(tot,10+i,1000000000);
                add(tot,10+j,1000000000);
                add(s,tot,w[i][j]+w[j][i]);
                all+=w[i][j]+w[j][i];
            }
        }
        all-=dinic();
        printf("Case #%d: %d\n",++N,all);
    }
}