思路:
最小树形图就是在有向图中,以一个点为根的有向生成树,并且边权值和最小。
思路大体上是贪心。
1、贪心地选取每个点最小的入边。
2、存在环则缩点,将其余边权值减去这个边终点的最小入边值。因为这相当于删除这个边的原来入边,再加上另外一条边。
3、重复1,2直到无环。
复杂度为
O
(
N
M
)
O(NM)
O(NM)。Tarjan更快的算法还不会。
代码:
#include<bits/stdc++.h>
#define pii pair<int,int>
#define ll long long
#define cl(x,y) memset(x,y,sizeof(x))
#define loop(x,y,z) for(x=y;x<=z;x++)
#define reve(x,y,z) for(x=y;x>=z;x--)
#define ct cerr<<"Time elapsed:"<<1.0*clock()/CLOCKS_PER_SEC<<"s.\n";
#define mp make_pair
#define pb push_back
#define fi first
#define se second
#define all(x) x.begin(),x.end()
#define lson x<<1,l,mid
#define rson x<<1|1,mid+1,r
#define INF 1e18
const int N=1e6+10;
const int mod=1e9+7;
const int inf=0x3f3f3f3f;
const double eps=1e-8;
const double pi=acos(-1);
using namespace std;
struct edge
{
int u,v,w;
}e[N];
int in[N],pre[N],id[N];//in最小入边权,pre对应起点,id结点所属环的编号
int n,m,vis[N];
int ZL(int rt)
{
int ans=0,i;
while(1)
{
for(i=1;i<=n;i++)
{
in[i]=inf;
vis[i]=id[i]=0;
}
for(i=1;i<=m;i++)
{
int u=e[i].u,v=e[i].v;
if(e[i].w<in[v] && u!=v)//遍历所有边,找每个边的最小入边
{
pre[v]=u;
in[v]=e[i].w;
}
}
for(i=1;i<=n;i++)
if(i!=rt && in[i]==inf)
return -1;//除根之外仍存在孤立点
//找环
int cnt=0;//记录环数
in[rt]=0;
for(i=1;i<=n;i++)
{
ans+=in[i];
int v=i,u;
while(v!=rt && vis[v]!=i && !id[v])//找环
//遍历到根,找到相同标记的点,结点属于其余环
{
vis[v]=i;
v=pre[v];
}
if(v!=rt && !id[v])//环上的点标记为同一点
{
id[v]=++cnt;
for(u=pre[v];u!=v;u=pre[u])
id[u]=cnt;
}
}
if(cnt==0)//无环
break;
for(i=1;i<=n;i++)
if(!id[i])
id[i]=++cnt;
for(i=1;i<=m;i++)//缩点,修改边权
{
int u=e[i].u,v=e[i].v;
e[i].u=id[u];
e[i].v=id[v];
if(id[u]!=id[v])
e[i].w-=in[v];
}
rt=id[rt];
n=cnt;
}
return ans;
}
int main()
{
ios::sync_with_stdio(false);
cin.tie(0);cout.tie(0);
int i,rt;
cin>>n>>m>>rt;
for(i=1;i<=m;i++)
{
cin>>e[i].u>>e[i].v>>e[i].w;
if(e[i].u==e[i].v)
e[i].w=inf;//去自环
}
cout<<ZL(rt)<<endl;
return 0;
}