传送门
题目同上
史上速度最快ac
#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<queue>
using namespace std;
const int maxn=2005;//点的最大上线
const int maxx=44005;//边的最大上线
const int INF=1<<30;
int to[maxx],nex[maxx],cap[maxx],flow[maxx],cost[maxx];
int head[maxn];
int pre[maxn],dis[maxn];
int eid;
void init()//每个样例都需要使用一次
{
eid=0;
memset(head,-1,sizeof(head));
}
void addEdge(int u,int v,int ca,int co)//建边
{
int i;
for(i=head[u];~i;i=nex[i])
if(to[i]==v)
break;
if(~i)
{
if(cost[i]>co)
cost[i]=co,cost[i^1]=-co;
return;
}
to[eid]=v;
cap[eid]=ca;
cost[eid]=co;
flow[eid]=0;
nex[eid]=head[u];
head[u]=eid++;
to[eid]=u;
cap[eid]=0;
cost[eid]=-co;
flow[eid]=0;
nex[eid]=head[v];
head[v]=eid++;
}
bool inQ[maxn];
bool spfa(int s,int e,int n)//最短路算法代替广搜
{
queue<int> que;
memset(inQ,false,sizeof(inQ));
for(int i=0;i<=n;i++)
dis[i]=INF;
dis[s]=0,inQ[s]=true;
pre[s]=-1;
que.push(s);
while(!que.empty())
{
int u=que.front();
que.pop();
inQ[u]=false;
for(int i=head[u];~i;i=nex[i])
{
int v=to[i];
if(cap[i]-flow[i]>0&&dis[v]>dis[u]+cost[i])
{
dis[v]=dis[u]+cost[i];
pre[v]=i;//这是为了修改边的容量
if(!inQ[v])
que.push(v),inQ[v]=true;
}
}
}
return dis[e]!=INF;
}
int minCostMaxFlow(int s,int e,int &minCost,int n)//建完边使用
{
int ans=0;
while(spfa(s,e,n))
{
int mint=INF;
for(int i=pre[e];~i;i=pre[to[i^1]])
{
if(mint>cap[i]-flow[i])
mint=cap[i]-flow[i];
}
ans+=mint;
for(int i=pre[e];~i;i=pre[to[i^1]])
{
flow[i]+=mint;flow[i^1]-=mint;
minCost+=mint*cost[i];
}
}
return ans;
}
int main()
{
int t;
scanf("%d",&t);
int cal=1;
int n,m;
int x,y,w;
while(t--)
{
scanf("%d%d",&n,&m);
init();
while(m--)
{
scanf("%d%d%d",&x,&y,&w);
addEdge(x,n+y,1,w);
addEdge(y,n+x,1,w);
}
for(int i=1;i<=n;i++)
addEdge(0,i,1,0),addEdge(n+i,n+n+1,1,0);
int minCost=0;
int ans=minCostMaxFlow(0,n+n+1,minCost,n+n+2);
if(ans<n)
printf("Case %d: NO\n",cal++);
else
printf("Case %d: %d\n",cal++,minCost);
}
return 0;
}