• 【bzoj 4675】 点对游戏


    题目

    发现一个人如果最终拿走了(k)个点,那么这个人的答案就是

    [frac{inom{n-2}{k-2}sum_{i=1}^{n}sum_{j=1}^{n}[dis(i,j)in M]}{inom{n}{k}} ]

    考虑每一个点对对答案的贡献,我们枚举一个点对,之后对于剩下的(n-2)个点我们随便选出(k-2)对点组成我们所选择的(k)个点,只需要看看这个点对的距离是否属于(M)就好了

    组合数划开发现是(frac{k imes(k-1)}{n imes (n-1)}),于是我们可以直接长链剖分求后面那个柿子就好了

    代码

    #include<algorithm>
    #include<iostream>
    #include<cstring>
    #include<cstdio>
    #define re register
    #define LL long long
    #define max(a,b) ((a)>(b)?(a):(b))
    #define min(a,b) ((a)<(b)?(a):(b))
    const int maxn=50005;
    inline int read() {
    	char c=getchar();int x=0;while(c<'0'||x>'9') c=getchar();
    	while(c>='0'&&c<='9') x=(x<<3)+(x<<1)+c-48,c=getchar();return x;
    }
    struct E{int v,nxt;}e[maxn<<1];
    int len[maxn],n,m,son[maxn];
    int deep[maxn],head[maxn];
    int tax[maxn],*id=tax,*f[maxn];
    LL ans;
    int num,a[15];
    inline void add(int x,int y) {
    	e[++num].v=y;e[num].nxt=head[x];head[x]=num;
    }
    void dfs1(int x) {
    	for(re int i=head[x];i;i=e[i].nxt) {
    		if(deep[e[i].v]) continue;
    		deep[e[i].v]=deep[x]+1;
    		dfs1(e[i].v);
    		if(len[e[i].v]>len[son[x]]) son[x]=e[i].v;
    	}
    	len[x]=len[son[x]]+1;
    }
    void dfs(int x) {
    	f[x][0]=1;
    	if(son[x]) {
    		f[son[x]]=f[x]+1,dfs(son[x]);
    		for(re int k=1;k<=m;k++)
    			if(a[k]<len[x]) ans+=f[x][a[k]];	
    	}
    	for(re int i=head[x];i;i=e[i].nxt) {
    		if(deep[e[i].v]<deep[x]||e[i].v==son[x]) continue;
    		f[e[i].v]=id;id+=len[e[i].v];dfs(e[i].v);
    		for(re int j=0;j<len[e[i].v];j++) 
    			for(re int k=1;k<=m;k++)
    				if(a[k]-j-1>=0&&a[k]-j-1<len[x]) ans+=f[e[i].v][j]*f[x][a[k]-j-1];
    		for(re int j=0;j<len[e[i].v];j++)
    			f[x][j+1]+=f[e[i].v][j];
    	}
    }
    inline double C(double n,double m) {
    	double t=m*(m-1),b=n*(n-1);
    	return t/b;
    }
    int main() {
    	n=read(),m=read();
    	for(re int i=1;i<=m;i++) a[i]=read();
    	for(re int x,y,i=1;i<n;i++) 
    		x=read(),y=read(),add(x,y),add(y,x);
    	deep[1]=1,dfs1(1);
    	f[1]=id;id+=len[1];dfs(1);
    	int k[3];int h=n,now=0;k[0]=k[1]=k[2]=0;
    	while(h) {k[now]++;now++;h--;now%=3;}
    	for(re int i=0;i<3;i++)
    	if(k[i]>=2) printf("%.2lf
    ",C(n,k[i])*ans);
    		else puts("0.00");
    	return 0;
    }
    
  • 相关阅读:
    启动 Eclipse 弹出“Failed to load the JNI shared library jvm.dll”错误的解决方法!
    处理textarea的空格和换行
    微信内置页实现自动拨号到页面的功能
    兼容ie的渐变css
    H5压缩图片上传(FileReader +canvas)
    css样式写一个小三角
    js禁止选中图片和文字
    js解决浏览器打印自动分页的问题
    IE8透明度问题
    在css加载完毕后自动判断页面是否加入css或js文件
  • 原文地址:https://www.cnblogs.com/asuldb/p/10683128.html
Copyright © 2020-2023  润新知