CodeForces 288 D.Polo the Penguin and Trees(组合数学)

Description

给出一棵 n 个节点的树,问四元组(a,b,c,d)个数使其满足:

1. 1a<bn,1c<dn

2.树上从 a b的简单路径与 c d的简单路径不经过相同节点

Input

第一行一整数 n 表示点数,之后n1行每行两个整数 u,v 表示一条树边 (1n80000)

Output

输出满足条件的四元组个数

Sample Input

4
1 2
2 3
3 4

Sample Output

2

Solution

满足第一个条件的四元组数量为 C2nC2n ,考虑求不满足第二个条件的四元组数量,即两条路径有公共点,考虑每个点 u 作为两条路径交点中深度最小的公共点对不合法答案的贡献,记size[u]为以 u 为根的子树个数,num1为经过点 u 且一端在u的子树中,另一端不在 u 子树中的路径数量,num2表示经过点 u 且两端都在u的子树中的路径数量,这两个值都比较好求

num2=C2size[u]fa[v]=uC2size[v]

num1=C2nfa[v]=uC2sizevC2size[fa[u]]num2

由于不可能存在两条路径使得这两条路径都是一端在 u 子树中一端不在,因为这样的话这两条路径公共点中深度最小的不会是u,至少也是 fa[u] ,故这两条路径必然是两种情况:一条路径一段在 u 子树一段不在另一条路径两端都在,或者是两条路径都是两端都在u子树中,故方案数为 2num1num2+num2num2

Code

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#include<queue>
#include<map>
#include<set>
#include<ctime>
using namespace std;
typedef long long ll;
typedef pair<int,int>P;
const int INF=0x3f3f3f3f,maxn=80005;
int n,size[maxn],fa[maxn];
vector<int>g[maxn]; 
void dfs(int u,int f)
{
    size[u]=1;
    fa[u]=f;
    for(int i=0;i<g[u].size();i++)
    {
        int v=g[u][i];
        if(v==f)continue;
        dfs(v,u);
        size[u]+=size[v];
    }
}
ll C(int n)
{
    return (ll)n*(n-1)/2;
}
int main()
{
    scanf("%d",&n);
    for(int i=1;i<n;i++)
    {
        int u,v;
        scanf("%d%d",&u,&v);
        g[u].push_back(v),g[v].push_back(u);
    }
    dfs(1,0);
    ll ans=C(n)*C(n);

    for(int u=1;u<=n;u++)
    {
        ll num1=C(n),num2=C(size[u]);
        for(int i=0;i<g[u].size();i++)
        {
            int v=g[u][i];
            if(v!=fa[u])num1-=C(size[v]),num2-=C(size[v]);
            else num1-=C(n-size[u]);
        }
        ans-=(2ll*num1*num2-num2*num2);
    } 
    printf("%I64d\n",ans);
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值