2016 Multi-University Training Contest 4 T9

http://acm.hdu.edu.cn/showproblem.php?pid=5772

最大权闭合子图。

得到价值w[i][j]的条件是选了i，j这两个位置的字符。选择位置的i字符花费为

第一次选s[i]：a[s[i]]       不是第一次选s[i]：b[s[i]]

所以对于选i位置字符前提为选了花费为b[s[i]]-a[s[i]]的字符i。

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

 #include<cstdio>
 #include<cstring>
 #define M 120010
 #define N 6010
 #define inf 1000000000
 int len[M <<],e[M <<],nex[M <<],other[M <<],head[N],last[N],d[N],num[N];
 int te,sum,_,j,num1,num2,w[][],a[],b[],m,n,i,ans,tot,ss,tt,ee,u,v,c;
 char s[];
 int init()
 {
     memset(head,,sizeof(head));
     memset(num,,sizeof(num));
     memset(d,,sizeof(d));
     ans=ee=;
 }
 int min(int a,int b)
 {
     return a<b?a:b;
 }
 void add(int u,int v,int c)
 {
     //printf("%d %d %d\n",u,v,c);
     other[++ee]=ee+;
     e[ee]=v;nex[ee]=head[u];head[u]=ee;len[ee]=c;
     other[++ee]=ee-;
     e[ee]=u;nex[ee]=head[v];head[v]=ee;len[ee]=;
 }
 int dfs(int x,int flow)
 {
     int rec,j,p;
     if (x==tt) return flow;
     rec=;j=last[x];
     while (j!=)
     {
         if (len[j]> && d[x]==d[e[j]]+)
         {
             last[x]=j;
             p=dfs(e[j],min(len[j],flow-rec));
             len[j]-=p;len[other[j]]+=p;
             rec+=p;
             if (rec==flow) return rec;
         }
     j=nex[j];
     }
     if (d[ss]>tot) return rec;
      if (--num[d[x]]==) d[ss]=tot;
      last[x]=head[x];
      num[++d[x]]++;
      return rec;
 }
 int main()
 {
     scanf("%d",&_);
     while (_--)
     {
         sum=;init();
         scanf("%d",&n);
         scanf("%s",s+);
         for (i=;i<=;i++)
             scanf("%d%d",&a[i],&b[i]);
         for (i=;i<=n;i++)
             for (j=;j<=n;j++)
                 scanf("%d",&w[i][j]);
         num1=;ss=;
         for (i=;i<=n;i++)
             for (j=i+;j<=n;j++)
                 if (w[i][j]+w[j][i]!=)
                 {
                     add(ss,++num1,w[i][j]+w[j][i]);
                     sum+=w[i][j]+w[j][i];
                 }
         num2=;tt=num1+n+;
         for (i=;i<=n;i++)
             for (j=i+;j<=n;j++)
                    if (w[i][j]+w[j][i]!=)
                   {
                     add(++num2,num1+i,inf);
                     add(num2,num1+j,inf);
                 }
         for (i=;i<=n;i++)
         {
             add(num1+i,num1+n+s[i]-''+,inf);
             add(num1+i,tt,a[s[i]-'']);
         }
         for (i=;i<=;i++)
             add(num1+n+i+,tt,b[i]-a[i]);
         tot=num[]=tt;
         for (i=ss;i<=tt;i++)
             last[i]=head[i];
         while (d[ss]<tot)
             ans+=dfs(ss,);
         printf("Case #%d: %d\n",++te,sum-ans);
     }
     return ;

 }

Run ID	Submit Time	Judge Status	Pro.ID	Exe.Time	Exe.Memory	Code Len.	Language	Author
17819971	2016-07-29 15:40:45	Accepted	5772	15MS	2052K	2416 B	G++	lbz007