分治算法

最新推荐文章于 2024-04-11 09:50:43 发布

原创最新推荐文章于 2024-04-11 09:50:43 发布 · 112 阅读

0 ·

CC 4.0 BY-SA版权

c++ 同时被 2 个专栏收录

186 篇文章

订阅专栏

算法

166 篇文章

订阅专栏

本文深入探讨了一种基于分治思想的大规模模拟算法实现，通过具体代码示例，详细讲解了如何利用递归和结构化数据处理复杂网络问题，特别关注于节点间的权重计算和路径选择，为理解高级算法提供了一个实践视角。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

可能我太菜了，众多题解的t0*t0 +t1*t2*2始终没看懂

分治大模拟好！

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int INF = 1000000000;
const int M=20010;
struct Node{ int next,to,w;}e[M*2];
int n,head[M*2],siz[M],cnt,sum,mp,rt,v[M],dis[M],judge[5];
ll res=0,ans=0,tem;
ll gcd(ll a,ll b){
    if(a<b)swap(a,b);
	return b?gcd(b,a%b):a;
}
void add(int x,int to,int w){
    e[cnt].to=to;
    e[cnt].w=w;
    e[cnt].next=head[x];
    head[x]=cnt++;
}
void getrt(int x,int fa){
    siz[x]=1;
    int mm=0;
    for(int to,i=head[x];~i;i=e[i].next){
        to=e[i].to;
        if(v[to]||fa==to)continue;
        getrt(to,x);
        siz[x]+=siz[to];
        mm=max(mm,siz[to]);
    }
    mm=max(sum-siz[x],mm);
    if(mp>mm)mp=mm,rt=x;
}
void dfs(int x,int fa){
    siz[x]=1;
    int mm=0;
    for(int to,i=head[x];~i;i=e[i].next){
        to=e[i].to;
        if(v[to]||fa==to)continue;
        dfs(to,x);
        siz[x]+=siz[to];
    }
}
void getdis(int x,int fa,int w){
     for(int to,i=head[x];~i;i=e[i].next){
        to=e[i].to;
        if(v[to]||to==fa)continue;
        dis[++dis[0]]=(e[i].w+w)%3;
        getdis(to,x,dis[dis[0]]);
    }
}

void work(int x){
    v[x]=1;
    judge[0]=1;
    for(int to,i=head[x];~i;i=e[i].next){
        to=e[i].to;
        if(v[to])continue;

        dis[dis[0]=1]=e[i].w;
        getdis(to,x,dis[dis[0]]);

        for(int j=dis[0];j;j--)res+=judge[(3-dis[j])%3];
        for(int j=dis[0];j;j--)judge[dis[j]]++;
    }
    memset(judge,0,sizeof(judge));
}

void slove(int x){
    work(x);
    for(int to,i=head[x];~i;i=e[i].next){
        to=e[i].to;
        if(v[to])continue;
        sum=siz[to];mp=INF;
        getrt(to,0);
        dfs(rt,0);
        slove(rt);
    }
}
int main(){
    cnt=1;memset(head,-1,sizeof(head));
    scanf("%d",&n);
    for(int x,y,z,i=1;i<n;i++){
        scanf("%d%d%d",&x,&y,&z);
        z%=3;
        add(x,y,z);
        add(y,x,z);
    }
    mp=INF;sum=n;

    getrt(1,0);
    dfs(rt,0);

    slove(rt);
    res=res+res+n;
    ans=n*n;
    tem=gcd(ans,res);
    cout<<res/tem<<"/"<<ans/tem;
    return 0;
}