题面
题解
这是一道很裸很裸的最大权值闭合子图。然而被暴力蒙蔽了双眼的我真的没看出来,我还是太菜了!
一看cost的定义,以及分段的方式,非常符合网络流的特性。
构图:
- 对于 wi,j + wj,i 建一个点,
- 对于每个位置建一个点,代价为 −a
- 对于每个字符建一个点,代价为 −b + a
- 一类点连向二类点,二类点连向三类点
代码
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define MAXN 450000
#define MAXM 120
int End[MAXN],Next[MAXN],Last[MAXN],Len[MAXN],e=1;
void add(int x,int y,int L)
{
Len[++e]=L;
End[e]=y;
Next[e]=Last[x];
Last[x]=e;
return ;
}
void _r(int& x)
{
char c=getchar();
while(c<'0'||c>'9')
{
c=getchar();
}
for(x=0;c>='0'&&c<='9';c=getchar())
{
x=(x<<1)+(x<<3)+c-'0';
}
return ;
}
int A[MAXM],B[MAXM],w[MAXM][MAXM],tot,S,E,n;
char s[MAXM];
int dis[MAXN],cnt[MAXN],id[MAXM][MAXM],pp[MAXM],qq[20];
const int inf=0x3c3c3c3c;
int dfs(int u,int flow)
{
int delta=0,tmp,v;
if(u==E)
{
return flow;
}
for(int t=Last[u];t;t=Next[t])
{
v=End[t];
if(Len[t]>0&&dis[u]==dis[v]+1)
{
tmp=dfs(v,min(flow,Len[t]));
delta+=tmp;
Len[t]-=tmp;
Len[t^1]+=tmp;
if(delta==flow||dis[S]>=tot)
{
return delta;
}
}
}
if(dis[S]>=tot)
{
return delta;
}
cnt[dis[u]]--;
if(cnt[dis[u]]==0)
{
dis[S]=tot;
}
cnt[++dis[u]]++;
return delta;
}
void work()
{
memset(dis,0,sizeof(dis));
memset(cnt,0,sizeof(cnt));
e=1;
memset(Last,0,sizeof(Last));
_r(n);
tot=0;
scanf("%s",s+1);
for(int i=0;i<10;i++)
{
_r(A[i]);
_r(B[i]);
}
S=++tot;
E=++tot;
int ans=0;
for(int i=0;i<10;i++)
{
qq[i]=++tot;
add(tot,E,B[i]-A[i]);
add(E,tot,0);
}
for(int i=1;i<=n;i++)
{
pp[i]=++tot;
add(tot,E,A[s[i]-'0']);
add(E,tot,0);
add(tot,qq[s[i]-'0'],inf);
add(qq[s[i]-'0'],tot,0);
}
for(int i=1;i<=n;i++)
{
for(int j=1;j<=n;j++)
{
_r(w[i][j]);
if(i!=j)
{
ans+=w[i][j];
}
}
}
for(int i=1;i<=n;i++)
{
for(int j=i+1;j<=n;j++)
{
id[i][j]=++tot;
add(S,tot,w[i][j]+w[j][i]);
add(tot,S,0);
add(tot,pp[i],inf);
add(pp[i],tot,0);
add(tot,pp[j],inf);
add(pp[j],tot,0);
}
}
while(dis[S]<tot)
{
ans-=dfs(S,inf);
}
printf("%d\n",ans);
return ;
}
int main()
{
int T;
scanf("%d",&T);
for(int i=1;i<=T;i++)
{
work();
}
return 0;
}