创作背景
(神秘的声音):学最小生成树了
最小生成树原理
一个无向图有n个结点和m条边
其中连接结点u和结点v的边的权值为wu,v(不含自环)
从这m条边中选择n-1条边,使得这n个点连通,这样得到的树称为该无向图的(一颗)生成树。
问如何选边使得生成树的权值和最小。
Prim算法(无限接近贪心)
Prim算法简介
一个无向图有n个结点和m条边
其中连接结点u和结点v的边的权值为wu,v(不含自环)
从这m条边中选择n-1条边,使得这n个点连通,这样得到的树称为该无向图的(一颗)生成树。
问如何选边使得生成树的权值和最小。
模拟讲解
看不懂?来,模拟一下 你就知道
首先,我们确定一个点1,准备从1开始
此刻,将与1联通的3个节点更新权值,准备下一步的贪心。
好的,现在将1存进来,根据贪心,3节点的权值最小,选择,其他节点的值也要更新
好的,现在将3存进来,根据贪心,6节点的权值最小,选择,其他节点的值也要更新
好的,现在将6存进来,根据贪心,4节点的权值最小,选择,其他节点的值也要更新(貌似没有其他节点了)
好了,现在把4存进来,没有节点了,怎么办,是时候叫出超级飞侠了 回溯
根据贪心,2节点的权值最小,选择。
好了,现在把2存进来了,只剩下5了,选择,存入。
Code
#include<bits/stdc++.h>
using namespace std;
typedef pair<int, int> pr;
const int maxn=10005;
int m,n,minn[maxn],u,v,w;
bool vis[maxn];
vector <pr> g[10001];
struct g
{
int from, to, dis;
}G[2*maxn];
int prim()
{
memset(vis, false, sizeof(vis));
for(int i=1;i<=n;i++)
minn[i]=1e9;
minn[1]=0;
for(int i = 1;i < n; i++)
{
int u=-1;
for (int j = 1; j <= n; j++)
if(!vis[j] && (u == -1 || minn[j] < minn[u]))
u=j;
vis[u]= true;
for(int j=0;j<g[u].size();j++)
{
int v = g[u][j].first, w = g[u][j].second;
if (!vis[v] && w < minn[v])
minn[v] = w;
}
}
int sum = 0;
for (int i = 1; i <= n; i++)
sum += minn[i];
return sum;
}
int main()
{
cin>>n>>m;
for(int j=1;j<=m;j++)
{
cin>>u>>v>>w;
g[u].push_back(pr(v,w));
g[v].push_back(pr(u,w));
}
int ans=prim();
cout<<ans;
}
Kruskal算法
Kruskal算法实现
看看简介:
首先将边全部按边权升序排序,并设包含原图全部顶点却不包含任何边的图为G,然后按顺序将边加入图G,当轮到边p时,如果p的两端已经在同一个连通块,就忽略边p,否则就加入边p使两端的顶点连通。加入的边即为最小生成树的边。
是不是很生涩很难懂?其实就3个词:
边权升序排序:
s
o
r
t
sort
sort
两个点是否在一个连通块:
宽
搜
,
并
查
集
宽搜,并查集
宽搜,并查集
并查集忘记的戳这里复习一下
所以:
#include<cstdio>
#include<iostream>
#include<algorithm>
using namespace std;
int ans,sum=0,m,x,i,j;
struct point
{
int x;
int y;
int v;
};
point a[99001];
int f[1001];
int find(int x)
{
if(f[x]!=x) f[x]=find(f[x]);
return f[x];
}
void join(int x,int y)
{
int fa=find(x);
int fb=find(y);
if(fa != fb) f[fa]=fb;
}
int cmp(const point &x,const point &y)
{
if(x.v<y.v) return 1;
else return 0;
}
inline int read()
{
int x=0,f=1;char ch=getchar();
while (!isdigit(ch)){if (ch=='-') f=-1;ch=getchar();}
while (isdigit(ch)){x=x*10+ch-48;ch=getchar();}
return x*f;
}
int main()
{
while(1)
{
int n=read();
for(i=1;i<=n;+i++)
for(j=1;j<=n;j++)
{
int x=read();
if(x != 0)
{
m++;
a[m].x=i;a[m].y=j;a[m].v=x;
}
}
for(i=1;i<=n;i++)
f[i]=i;
sort(a+1,a+m+1,cmp);
for(i=1;i<=m;i++)
{
if(find(a[i].x) !=find(a[i].y))
{
join(a[i].x,a[i].y);
sum+=a[i].v;
ans++;
}
if(ans==n-1) break;
}
cout<<sum;
sum=0;
ans=0;
}
return 0;
}
The End
额……这种算法难度应该也还好吧……溜了溜了~