luogu P1340 兽径管理

分析

题意非常明确，就是在不断加边的过程中跑MST（求最小生成树）
Kruskal或prim
由于题是要不断加边，若用prim就必须在线操作，而kruskal则可以离线一遍sort就行
于是这道题的最优解应该就是kruskal了

code

#include<bits/stdc++.h>
using namespace std;
#define loop(i,start,end) for(int i=start;i<=end;++i)
#define anti_loop(i,start,end) for(int i=start;i>=end;--i)
#define clean(arry,num); memset(arry,num,sizeof(arry));
#define asc(a) (a-'A')
const int maxw=6000+10;
const int maxn=200+10;
struct ed{int s;int e;int v;int ord;}edge[maxw];
int fa[maxn];
int n,w;
inline int read()
{
   int ans=0;
   char r=getchar();bool neg=false;
   while(r<'0'||r>'9'){if(r=='-')neg=true;r=getchar();}
   while(r>='0'&&r<='9'){ans=ans*10+r-'0';r=getchar();}
   return (neg)?-ans:ans;
}
int findfa(int u)
{
    return fa[u]=(fa[u]==u)?u:findfa(fa[u]);
}
inline void sum(int a,int b)
{
    fa[findfa(a)]=findfa(b);
}
inline bool cmp(ed a,ed b)
{
    return a.v<b.v;
}
void write(int x){
    if(x<0){putchar('-');x=-x;}
    if(x>9)write(x/10);
    putchar(x%10+'0');
}
inline void kruskal(int cnt)
{
    int finish_e=0;
    int val=0;
    loop(i,1,w)
    {
        int u=edge[i].s,v=edge[i].e,vval=edge[i].v;
        if(findfa(u)==findfa(v)||edge[i].ord>cnt)continue;
        sum(u,v);
        val+=vval;
        ++finish_e;
    }
    if(finish_e+1<n){write(-1);printf("\n");}
    else {write(val);printf("\n");}
}
int main()
{
    //freopen("datain.txt","r",stdin);
    n=read();w=read();
    loop(cnt,1,w)
    {
        int s,e,v;s=read();e=read();v=read();
        edge[cnt].s=s;
        edge[cnt].e=e;
        edge[cnt].v=v;
        edge[cnt].ord=cnt;
    }
    sort(edge+1,edge+1+w,cmp);
    loop(cnt,1,w)
    {
        if(cnt+1<n){write(-1);printf("\n");continue;}
        loop(i,1,n)fa[i]=i;
        kruskal(cnt);
    }
    return 0;
}

学到的东西

• 回顾MST相关算法

1.kruskal
贪边，对边排序，然后用并查集维护入选边所在的集合，选到（n-1）条边时，求得当下MST

Kruskal算法的时间复杂度
Kruskal算法每次要从都要从剩余的边中选取一个最小的边。通常我们要先对边按权值从小到大排序，这一步的时间复杂度为为O(|Elog|E|)。Kruskal算法的实现通常使用并查集，来快速判断两个顶点是否属于同一个集合。最坏的情况可能要枚举完所有的边，此时要循环|E|次，所以这一步的时间复杂度为O(|E|α(V))，其中α为反Ackermann函数，其增长非常慢，我们可以视为常数。所以Kruskal算法的时间复杂度为O(|Elog|E|)。

2.prim
prim算法以某个点为起点，向四周的未知领域拓展，利用距离信息，以三角形不等式为判断依据更新拓展到的未知点距离，选入n个点时，完成算法

prim算法的时间复杂度

由于prim算法是以图的顶点起始点，每次选取所有顶点上对应的权值最小的边进行构建，所以prim的时间开销和边无关，对于定点数为n时prim的时间复杂度为 O(n^2)，所以prim算法更适合求解边数很多的稠密图的MST。