Phone Call HDU - 6074

http://acm.hdu.edu.cn/showproblem.php?pid=6074

看完题解感觉很奇妙的一道题

和求最小生成树很相像! 先对m个操作按权值排序 可以想到每一次将abcd这四个点表示的两条链暴力合并在一起 最后看1所在连通块的总权值是多少 但是需要优化

可以开两个数组 f[i]就是普通并查集的代表元素 up[i]就是用来找i之上的一个与i不联通的点的 也是用并查集来维护 和f数组一样 up数组一开始也是指向自己 然后对每一次操作 求ab的lca 将a到lca这条链上的点都通过f合并在一起 同时up也通过路径压缩不断更新 这样每一次向上爬链时都是并查集路径压缩的复杂度 b到lca以及cd路径的操作也是一样的 最后再把两条链并查集合并一下就好

 

#include <cstdio>
#include <cmath>
#include <cstring>
#include <algorithm>
using namespace std;
#define ll long long

struct node1
{
    int v;
    int next;
};

struct node2
{
    int a;
    int b;
    int c;
    int d;
    ll w;
};

node1 edge[200010];
node2 pre[100010];
ll val[100010];
int dp[100010][20];
int first[100010],f[100010],up[100010],cnt[100010],deep[100010];
int n,m,num;

void addedge(int u,int v)
{
    edge[num].v=v;
    edge[num].next=first[u];
    first[u]=num++;
}

bool cmp(node2 n1,node2 n2)
{
    return n1.w<n2.w;
}

void dfs(int cur,int fa)
{
    int i,v;
    for(i=first[cur];i!=-1;i=edge[i].next)
    {
        v=edge[i].v;
        if(v!=fa)
        {
            dp[v][0]=cur;
            deep[v]=deep[cur]+1;
            dfs(v,cur);
        }
    }
    return;
}

void init()
{
    int i,j;
    dp[1][0]=0;
    deep[1]=1;
    dfs(1,-1);
    for(j=1;(1<<j)<=n;j++)
    {
        for(i=1;i<=n;i++)
        {
            dp[i][j]=dp[dp[i][j-1]][j-1];
        }
    }
    return;
}

int getlca(int u,int v)
{
    int i;
    if(deep[u]<deep[v]) swap(u,v);
    for(i=log2(n);i>=0;i--)
    {
        if(deep[dp[u][i]]>=deep[v])
        {
            u=dp[u][i];
        }
    }
    if(u==v) return u;
    for(i=log2(n);i>=0;i--)
    {
        if(dp[u][i]!=dp[v][i])
        {
            u=dp[u][i];
            v=dp[v][i];
        }
    }
    return dp[u][0];
}

int getf(int p)
{
    if(f[p]==p) return p;
    else
    {
        f[p]=getf(f[p]);
        return f[p];
    }
}

int getup(int p)
{
    if(up[p]==p) return p;
    else
    {
        up[p]=getup(up[p]);
        return up[p];
    }
}

void uniteII(int u,int v,ll w)
{
    int fu,fv;
    fu=getf(u),fv=getf(v);
    if(fu!=fv)
    {
        f[fv]=fu;
        cnt[fu]+=cnt[fv];
        val[fu]+=(val[fv]+w);
    }
}

void uniteI(int u,int v,ll w)
{
    while(1)
    {
        u=getup(u);
        if(deep[u]<=deep[v]) return;
        uniteII(u,dp[u][0],w);
        up[u]=dp[u][0];
    }
}

void solve(int a,int b,int c,int d,ll w)
{
    int lca;
    lca=getlca(a,b);
    uniteI(a,lca,w);
    uniteI(b,lca,w);
    lca=getlca(c,d);
    uniteI(c,lca,w);
    uniteI(d,lca,w);
    uniteII(a,c,w);
}

int main()
{
    int t,i,u,v;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d%d",&n,&m);
        memset(first,-1,sizeof(first));
        num=0;
        for(i=1;i<=n-1;i++)
        {
            scanf("%d%d",&u,&v);
            addedge(u,v);
            addedge(v,u);
        }
        init();
        for(i=1;i<=m;i++) scanf("%d%d%d%d%lld",&pre[i].a,&pre[i].b,&pre[i].c,&pre[i].d,&pre[i].w);
        sort(pre+1,pre+m+1,cmp);
        for(i=1;i<=n;i++)
        {
            f[i]=i;
            up[i]=i;
            cnt[i]=1;
            val[i]=0;
        }
        for(i=1;i<=m;i++) solve(pre[i].a,pre[i].b,pre[i].c,pre[i].d,pre[i].w);
        printf("%d %lld\n",cnt[getf(1)],val[getf(1)]);
    }
    return 0;
}

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值