【POJ3164】Command Network 最小树形图 朱刘算法

板子题不解释。。。

所谓最小树形图，直白一点的叫法就是“有向图最小生成树”，就是给定一个起点，在一张带权有向图中找一个边集的子集，使得从起点可以到达所有点（有向图哦），并且边集权值和最小。

主要思想就是为每个点找一条最小入边，如果不成环就是答案，如果成环就把环缩成一个点，并把｛（【指向属于该环的点】的边）的权值｝减去【（该点的最小入边）的权值】（你萌一定没有看懂），重复该过程直到不成环。

证明的话就类似如果要选一条指向该环的边，那么这个环就会断开那么相当于有一条边不选blablabla。。。。。

\

#include<cstdlib>
#include<cstdio>
#include<iostream>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<queue>
#include<vector>
using namespace std;
#define MAXN 10005
#define MAXM 10005
#define oo 1999999999
double _read()
{
	int x=0; char ch=getchar(); bool flag=false;
	while(ch<'0' || ch>'9'){if(ch=='-')flag=true; ch=getchar();}
	while(ch>='0'&&ch<='9'){x=x*10+ch-'0'; ch=getchar();} if(flag)x=-x; return (double)x;
}
int _readint()
{
	int x=0; char ch=getchar(); bool flag=false;
	while(ch<'0' || ch>'9'){if(ch=='-')flag=true; ch=getchar();}
	while(ch>='0'&&ch<='9'){x=x*10+ch-'0'; ch=getchar();} if(flag)x=-x; return x;
}
int N,M,f[MAXN],belong[MAXN],v[MAXN];
double in[MAXN];
struct P
{
	double x,y;
}A[MAXN];
struct EE
{
	int from,to; double d;
}E[MAXM];
void Init()
{
	int x,y; 
	for(int i=1;i<=N;i++){scanf("%lf%lf",&A[i].x,&A[i].y);}
	for(int i=1;i<=M;i++)
	{
		scanf("%d%d",&x,&y);
		E[i]=(EE){x,y,sqrt((A[x].x-A[y].x)*(A[x].x-A[y].x)+(A[x].y-A[y].y)*(A[x].y-A[y].y))};
	}
	return ;
}
double work(int &root)
{
	double ans=0; int cnt=0,j;
	while(1)
	{
		for(int i=1;i<=N;belong[i]=f[i]=-1,in[i]=oo,v[i]=-1,i++);
		for(int i=1;i<=M;i++)if(E[i].from!=E[i].to)
		{
			j=E[i].to;
			if(f[j]==-1 || in[j]>E[i].d)
			{
				f[j]=E[i].from; in[j]=E[i].d;
			}
		}
		for(int i=1;i<=N;i++)if(f[i]==-1 && i!=root )return -1;
		cnt=0; in[root]=0;
		for(int i=1;i<=N;i++)
		{
			ans+=in[i];
			j=i;
			while(j!=root && v[j]!=i && belong[j]==-1)
			{
				v[j]=i; j=f[j];
			}
			if(belong[j]==-1 && j!=root )
			{
				cnt++; belong[j]=cnt;
				for(int w=f[j];w!=j;w=f[w])
				{
					belong[w]=cnt;
				}
			}
		}
		if(cnt==0)break;
		for(int i=1;i<=N;i++)if(belong[i]==-1){cnt++;belong[i]=cnt;}
		for(int i=1;i<=M;i++)
		{
			if(belong[E[i].from]!=belong[E[i].to])E[i].d-=in[E[i].to];
			E[i].from=belong[E[i].from];
			E[i].to=belong[E[i].to];
		}
		N=cnt; root=belong[root];
	}
	return ans;
}
int main()
{
//	freopen("in.txt","r",stdin);
	while(scanf("%d%d",&N,&M)!=EOF)
	{
		Init();
		int rt=1; 
		double ans=work(rt);
		if(ans==-1)printf("poor snoopy\n");
		else printf("%.2lf\n",ans);
	}
	return 0;
}