jzoj 3514. 最小比例(ratio)

Description
图中共有N个点的完全图，每条边都有权值，每个点也有权值。要求选出M个点和M-1条边，构成一棵树，使得：
即所有边的权值与所有点的权值之和的比率最小。
给定N和M，以及N个点的权值，和所有的边权，要求M个点的最小比率生成树。
Input
第一行包含两个整数N和M(2<=N<=15，2<=M<=N)，表示点数和生成树的点数。
接下来一行N个整数，表示N个点的边权。
最后N行，每行N列，表示完全图中的边权。所有点权和边权都在[1,100]之间。
Output
输出最小比率生成树的M个点。当答案出现多种时，要求输出的第一个点的编号尽量小，第一个相同，则第二个点的编号尽量小，依次类推，中间用空格分开。编号从1开始。
Sample Input
输入1：
3 2
30 20 10
0 6 2
6 0 3
2 3 0
输入2：
2 2
1 1
0 2
2 0
Sample Output
输出1：
1 3
输出2：
1 2
Data Constraint
对于30%数据，N<=5。

//written by zzy

题目大意：

求最小(总边权/总点权)的生成树

题解：

考虑暴力枚举点的情况，跑最小生成树，更新答案。

#include<bits/stdc++.h>
#define N 1005

using namespace std;
int n,m,k,i,j,u,p,q,g,v;
int a[N],fa[N];
double ans=9999999997;

struct E {
	int x,y,z;
}l[N*N];

bool cmp(E a,E b) {
	return (a.z<=b.z);
}

int get(int x)
{
	if (x==fa[x]) return x;
	return fa[x]=get(fa[x]);
}

void merge(int x,int y)
{
	fa[get(x)]=get(y);
}

int mintree()
{
	int s=0;
	for (int i=1;i<=u;i++)
	 if (fa[l[i].x]!=0&&fa[l[i].y]!=0&&get(l[i].x)!=get(l[i].y))
	  {
	  	s+=l[i].z; merge(l[i].x,l[i].y);
	  }
	return s;
}

int main()
{
	freopen("ratio.in","r",stdin); 
	freopen("ratio.out","w",stdout); 
	scanf("%d%d",&n,&m);
	for (i=1;i<=n;i++)
	 scanf("%d",&a[i]);
	for (i=1;i<=n;i++)
	 for (j=1;j<=n;j++)
	  if (i!=j)
	  {
	  	u++;
	  	l[u].x=i; l[u].y=j; scanf("%d",&l[u].z);
	  } else scanf("%d",&v);
	sort(l+1,l+u+1,cmp);
	for (i=1;i<=(1<<n)-1;i++)
	 if (__builtin_popcount(i)==m)
	 {
	 	p=0; q=0;
	 	memset(fa,0,sizeof(fa));
	 	for (j=1;j<=n;j++)
	 	 if (i&1<<j-1) fa[j]=j,q+=a[j];
	 	p=mintree();
	 	if ((double)p/q<ans) {
	 		ans=(double)p/q; g=i;
	 	}
	 }
	for (j=1;j<=n;j++)
	 if (g&1<<j-1)
	  printf("%d ",j);
	return 0;
}