题意:
n*m的格子 每个格子有个数字 你可以任选起点 每一步向下或者向右 每一步的花费是曼哈顿距离-1 如果一步的两端的格子数字一样那么你会得到这个数字的花费 问 最多选k次起点 在遍历完所有格子前提下最多剩下多少花费
思路:
题目可以表示为用最多k条路径去覆盖n*m 思路就转向了二分图
X集合表示每个点的出 Y集合表示入 如果f格子能走到g格子 那么建边gx->fy流量为1费用为花费减得到
那么如何满足至多k个起点呢?
加入一个虚点 从这个点向所有Y集合连边流量1费用0
那么这个点可以认为是超级起点 从他开始走向真正起点的入
最后从源点S向所有X点连边流量1费用0 向虚点连边流量k费用0 从所有Y点向汇点T连边流量1费用0
这时最小费用流所得的费用的相反数即为剩下的花费 流即表示经过的点数
代码:
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
#define M 200
#define inf ((1U<<31)-1)
int t,cas,n,m,p,ansflow,ans,tot,S,T;
int head[M],pre[M],vis[M],dis[M],maz[M][M];
struct edge
{
int u,v,w,c,next;
}ed[M*M];
queue<int> qu;
void add(int U,int V,int W,int C)
{
ed[tot].u=U;
ed[tot].v=V;
ed[tot].w=W;
ed[tot].c=C;
ed[tot].next=head[U];
head[U]=tot++;
}
int spfa()
{
int i,u,v;
while(!qu.empty()) qu.pop();
for(i=0;i<M;i++)
{
vis[i]=0;
dis[i]=inf;
pre[i]=-1;
}
qu.push(S);
vis[S]=1;
dis[S]=0;
while(!qu.empty())
{
u=qu.front();
qu.pop();
vis[u]=0;
for(i=head[u];~i;i=ed[i].next)
{
if(!ed[i].w) continue;
v=ed[i].v;
if(dis[v]>dis[u]+ed[i].c)
{
dis[v]=dis[u]+ed[i].c;
pre[v]=i;
if(!vis[v])
{
vis[v]=1;
qu.push(v);
}
}
}
}
return dis[T]!=inf;
}
int mcmf()
{
int i,res=0,flow=0,tmp;
while(spfa())
{
tmp=inf;
for(i=pre[T];~i;i=pre[ed[i].u])
{
if(tmp>ed[i].w) tmp=ed[i].w;
}
for(i=pre[T];~i;i=pre[ed[i].u])
{
ed[i].w-=tmp;
ed[i^1].w+=tmp;
res+=ed[i].c;
}
flow+=tmp;
}
ansflow=flow;
return res;
}
int main()
{
int i,j,k,cost;
char input[110];
scanf("%d",&t);
while(t--)
{
printf("Case %d : ",++cas);
scanf("%d%d%d",&n,&m,&p);
for(i=1;i<=n;i++)
{
scanf("%s",input+1);
for(j=1;j<=m;j++) maz[i][j]=input[j]-'0';
}
S=0;
T=n*m*2+2;
memset(head,-1,sizeof(head));
tot=0;
add(S,n*m*2+1,p,0);
add(n*m*2+1,S,0,0);
for(i=1;i<=n;i++)
{
for(j=1;j<=m;j++)
{
add(S,(i-1)*m+j,1,0);
add((i-1)*m+j,S,0,0);
add(n*m+(i-1)*m+j,T,1,0);
add(T,n*m+(i-1)*m+j,0,0);
add(n*m*2+1,n*m+(i-1)*m+j,1,0);
add(n*m+(i-1)*m+j,n*m*2+1,0,0);
for(k=j+1;k<=m;k++)
{
cost=k-j-1;
if(maz[i][j]==maz[i][k]) cost-=maz[i][j];
add((i-1)*m+j,n*m+(i-1)*m+k,1,cost);
add(n*m+(i-1)*m+k,(i-1)*m+j,0,-cost);
}
for(k=i+1;k<=n;k++)
{
cost=k-i-1;
if(maz[i][j]==maz[k][j]) cost-=maz[i][j];
add((i-1)*m+j,n*m+(k-1)*m+j,1,cost);
add(n*m+(k-1)*m+j,(i-1)*m+j,0,-cost);
}
}
}
ans=mcmf();
if(n*m==ansflow) printf("%d\n",-ans);
else puts("-1");
}
return 0;
}