这个算法是我最近一直在学习的,刷了几天的题,现在也对这个算法有一些理解了【滑稽】
最小生成树可以用prime算法和kruskal两种算法解决,然而我只学会了kruskal算法QAQ,那么我就来浅谈一下最小生成树kruskal算法的定义与应用。
这个算法的思路分为四步:
***1.搜集图中所有的边
2.把搜集到的所有的边按照长度从小到大进行排序
3.从小到大选普通边,总是期望小的普通边能成为最小生成树的边。但是有个条件:新选进来的边不能和之前选中的边“闹矛盾”(就是不能和之前选中的边构成环),如何判断新选的边会不会和之前的边构成环,这个事情交给并查集来做。
4.直到选够了N-1条边为止*。
这样就可以做大部分的最小生成树的题目了(部分需要prime算法的题除外),只需要改一下建边的部分和加一些优化就可以了。
接下来我们看一些题吧~~~
1.pku1251
模板题,按照算法思路做即可,不过要先理解并查集的概念哦~
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<ctime>
#include<set>
#include<map>
#include<iostream>
#include<algorithm>
using namespace std;
struct node
{
int x,y,d;
}a[110000]; int len;
int cmp(const void *xx,const void *yy)
{
node n1=*(node *)xx;
node n2=*(node *)yy;
return n1.d-n2.d;
}
int fa[110000];
int findfa(int x)
{
if(fa[x]==x) return x;
else
{
fa[x]=findfa(fa[x]);
return fa[x];
}
}
int main()
{
int n,x,y,d,ans,t;
while(scanf("%d\n",&n)!=EOF)
{
if(n==0) break;
len=0;
for(int i=1;i<n;i++)
{
char c[10];scanf("%s",c+1);
x=c[1]-'A'+1;
int k;scanf("%d",&k);
for(int j=1;j<=k;j++)
{
scanf("%s%d",c+1,&d);
y=c[1]-'A'+1;len++;
a[len].x=x; a[len].y=y; a[len].d=d;
}
}
qsort(a+1,len,sizeof(node),cmp);
for(int i=1;i<=n;i++) fa[i]=i;
ans=0; t=0;
for(int i=1;i<=len;i++)
{
int fx,fy;
fx=findfa(a[i].x); fy=findfa(a[i].y);
if(fx!=fy)
{
ans+=a[i].d;
fa[fy]=fx;
t++; if(t==n-1) break;
}
}
printf("%d\n",ans);
}
return 0;
}
2.pku2075
输入的人名有些许恶心(捂嘴偷笑),只需要将每个人互相建边就可以了,最后判断能否在给定的长度中完成。
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<ctime>
#include<set>
#include<map>
#include<iostream>
#include<algorithm>
using namespace std;
struct node
{
int x,y;double d;
char st[21];
}a[110000]; int len;
int cmp(const void *xx,const void *yy)
{
node n1=*(node *)xx;
node n2=*(node *)yy;
if(n1.d<n2.d) return -1;
if(n1.d>n2.d) return 1;
return 0;
}
int fa[110000];
int findfa(int x)
{
if(fa[x]==x) return x;
else
{
fa[x]=findfa(fa[x]);
return fa[x];
}
}
int main()
{
double sum=0;
scanf("%lf",&sum);
int n,m; scanf("%d",&n);
len=0;
for(int i=1;i<=n;i++)
{
scanf("%s",a[i].st+1);
}
scanf("%d",&m);
for(int i=1;i<=m;i++)
{
char s1[21],s2[21]; double aa;int x,y;
scanf("%s%s%lf",s1+1,s2+1,&aa);
for(int j=1;j<=n;j++)
{
if(strcmp(a[j].st+1,s1+1)==0) {x=j;break;}
}
for(int j=1;j<=n;j++)
{
if(strcmp(a[j].st+1,s2+1)==0) {y=j;break;}
}
len++; a[len].x=x; a[len].y=y; a[len].d=aa;
}
qsort(a+1,len,sizeof(node),cmp);
for(int i=1;i<=n;i++) fa[i]=i;
double ans=0.0;int t=0;
for(int i=1;i<=len;i++)
{
int fx,fy;
fx=findfa(a[i].x);
fy=findfa(a[i].y);
if(fa[fx]!=fa[fy])
{
ans+=a[i].d;
fa[fy]=fx;
t++;if(t==n-1) break;
}
}
if(sum>=ans) printf("Need %.1lf miles of cable\n",ans);
else printf("Not enough cable\n");
return 0;
}
3.pku2485
同样的做法,但自己和自己就不需要建边了。
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<ctime>
#include<set>
#include<map>
#include<iostream>
#include<algorithm>
using namespace std;
struct node
{
int x,y,d;
}a[110000]; int len;
int cmp(const void *xx,const void *yy)
{
node n1=*(node *)xx;
node n2=*(node *)yy;
if(n1.d<n2.d) return -1;
if(n1.d>n2.d) return 1;
return 0;
}
int fa[110000];
int findfa(int x)
{
if(fa[x]==x) return x;
else
{
fa[x]=findfa(fa[x]);
return fa[x];
}
}
int main()
{
int t;scanf("%d",&t);
while(t--)
{
int n;scanf("%d",&n);
len=0;
for(int i=1;i<=n;i++)
{
for(int j=1;j<=n;j++)
{
int x; scanf("%d",&x);
if(i!=j)
{
len++; a[len].x=i; a[len].y=j; a[len].d=x;
}
}
}
qsort(a+1,len,sizeof(node),cmp);
for(int i=1;i<=n;i++) fa[i]=i;
int ans=0,t=0,maxx=0;
for(int i=1;i<=len;i++)
{
int fx,fy;
fx=findfa(a[i].x);
fy=findfa(a[i].y);
if(fa[fx]!=fa[fy])
{
ans+=a[i].d;
maxx=max(maxx,a[i].d);
fa[fy]=fx;
t++; if(t==n-1) break;
}
}
printf("%d\n",maxx);
}
return 0;
}
4.pku2421
这题给出了一些已经连好的边,那么我们在建边的时候把这些边的边权赋值为0即可。
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<ctime>
#include<set>
#include<map>
#include<iostream>
#include<algorithm>
using namespace std;
struct node
{
int x,y,d;
}a[110000]; int len;
int cmp(const void *xx,const void *yy)
{
node n1=*(node *)xx;
node n2=*(node *)yy;
if(n1.d<n2.d) return -1;
if(n1.d>n2.d) return 1;
return 0;
}
int fa[110000];
int findfa(int x)
{
if(fa[x]==x) return x;
else
{
fa[x]=findfa(fa[x]);
return fa[x];
}
}
int main()
{
int n;scanf("%d",&n);
len=0;
for(int i=1;i<=n;i++)
{
for(int j=1;j<=n;j++)
{
int x; scanf("%d",&x);
if(i!=j)
{
len++; a[len].x=i; a[len].y=j; a[len].d=x;
}
}
}
int m;scanf("%d",&m);
for(int i=1;i<=m;i++)
{
int x,y; scanf("%d%d",&x,&y);
for(int j=1;j<=len;j++)
{
if(a[j].x==x && a[j].y==y) {a[j].d=0; break;}
}
}
qsort(a+1,len,sizeof(node),cmp);
for(int i=1;i<=n;i++) fa[i]=i;
int ans=0,t=0;
for(int i=1;i<=len;i++)
{
int fx,fy;
fx=findfa(a[i].x);
fy=findfa(a[i].y);
if(fa[fx]!=fa[fy])
{
ans+=a[i].d;
fa[fy]=fx;
t++; if(t==n-1) break;
}
}
printf("%d\n",ans);
return 0;
}
5.pku1789
这道题有点坑的QAQ看了好久的题目才理解了题意。。。
稍稍想一下就可以得出解法了,我们将两个字符串中不同的字符个数作为建边的权值就可以了。同样,自己和自己也无需建边。
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<ctime>
#include<set>
#include<map>
#include<iostream>
#include<algorithm>
using namespace std;
struct node
{
int x,y,d;char st[10];
}a[4100000];int len;
int cmp(const void *xx,const void *yy)
{
node n1=*(node *)xx;
node n2=*(node *)yy;
if(n1.d>n2.d) return 1;
if(n1.d<n2.d) return -1;
return 0;
}
int fa[4100000];
int findfa(int x)
{
if(fa[x]==x) return x;
else
{
fa[x]=findfa(fa[x]);
return fa[x];
}
}
int main()
{
int n;
while(scanf("%d",&n)!=EOF)
{
if(n==0) break;
for(int i=1;i<=n;i++) scanf("%s",a[i].st+1);
len=0;
for(int i=1;i<=n;i++)
{
for(int j=i+1;j<=n;j++)
{
int s=0;
for(int k=1;k<=7;k++)
{
if(a[i].st[k]!=a[j].st[k]) s++;
}
len++; a[len].x=i; a[len].y=j; a[len].d=s;
}
}
qsort(a+1,len,sizeof(node),cmp);
for(int i=1;i<=n;i++) fa[i]=i;
int ans=0,t=0;
for(int i=1;i<=len;i++)
{
int fx,fy;
fx=findfa(a[i].x);
fy=findfa(a[i].y);
if(fa[fx]!=fa[fy])
{
ans+=a[i].d;
fa[fy]=fx;
t++; if(t==n-1) break;
}
}
printf("The highest possible quality is 1/%d.\n",ans);
}
return 0;
}
如果认真做完上述的几道题,我们会发现最小生成树并不难,但却很有用。希望以后你在联赛中遇到这类型的题,可以快速的将它AC。
下面给出几道其他的练习题,I hope that you can improve your programming skills soon.
caioj1090
caioj1091
codevs3315
codevs2627
codevs1003
codevs1519