[SDOI2006] 仓库管理员的烦恼 - 二分图最大权匹配

最小化代价,即最大化“本土”货物的数量

于是就是个二分图最大权匹配裸题

#include <bits/stdc++.h>
using namespace std;
#define reset(x) memset(x,0,sizeof x)
// Input: g[v][u] (v in II, u in I)
// Method: solve(n1,n2)
// Output: ans, mat[u] (u in I)
namespace km {
const int inf=1e+9;
const int MX=405;
int n,m;
int py[MX],vy[MX],pre[MX];
int slk[MX],g[MX][MX],kx[MX],ky[MX],ans;
int mat[MX];
void clear() {
n=m=0;
reset(py); reset(vy); reset(pre);
reset(slk); reset(g); reset(kx); reset(ky);
}
void KM(){
int i,j,k,x,p=0;
int d,t;
for(i=1;i<=n;i++)
for(j=1;j<=n;j++)
kx[i]=max(kx[i],g[i][j]);
for(i=1;i<=n;i++){
memset(vy,0,sizeof(int)*(n+1));
for(j=0;j<=n;j++) slk[j]=inf;
memset(pre,0,sizeof(int)*(n+1));
for(py[k=0]=i;py[k];k=p){
d=inf;vy[k]=1;x=py[k];
for(j=1;j<=n;j++)if(!vy[j]){
if((t=kx[x]+ky[j]-g[x][j])<slk[j])slk[j]=t,pre[j]=k;
if(slk[j]<d)d=slk[j],p=j;
}
for(j=0;j<=n;j++)
if(vy[j])kx[py[j]]-=d,ky[j]+=d;
else slk[j]-=d;
}
for(;k;k=pre[k])py[k]=py[pre[k]];
}
} void solve(int n1,int n2){
n=max(n1,n2);
KM();
ans=0;
for(int i=1;i<=n;i++)ans+=kx[i]+ky[i];
for(int i=1;i<=n1;i++)mat[i]=(g[py[i]][i]?py[i]:0);
}
} int main() {
int n,tot=0;
cin>>n;
for(int i=1;i<=n;i++) {
for(int j=1;j<=n;j++) {
cin>>km::g[i][j];
tot+=km::g[i][j];
}
}
km::solve(n,n);
cout<<tot-km::ans;
}
上一篇:通过jdbc获取数据库中的表结构 主键 各个表字段类型及应用生成实体类


下一篇:Python Web学习笔记之Python多线程和多进程、协程入门