[CSP-S模拟测试]:环(图论+期望)

若如初见. 提交于 2019-12-01 23:59:04

题目传送门(内部题79)


输入格式

  第一行读入两个整数$n,e$表示节点数及$cwystc$已确定的有向边边数。
  接下来$e$行,每行两个整数$x,y$描述$cwystc$确定的边。


输出格式

  输出一个整数表示期望陪妹子的天数。


样例

见下发文件


数据范围与提示

  对于$30\%$的数据:$n\leqslant 300$
  对于另外i$20\%$的数据:$e=0$
  对于另外$20\%$的数据:$e=\frac{n\times (n-1)}{2}$
  对于$100\%$的数据:$n\leqslant 100,000,e\leqslant 100,000$


题解

先转化一下题意,即让我们求竞赛图的期望最小环的个数。

而对于一张竞赛图,其要么是一张拓扑图,要么存在一个三元环。

因为该图不存在环,于是问题转化为三元环计数。

考虑容斥,不妨设$p[i]$表示点$i$的出度,那么对于三元组$(i,x,y)\{x\in p[i],y\in p[i],x\neq y\}$不组成三元环,而且这样的三元组只会在第$i$点枚举一次。

不妨再设$q[i]$表示连接点$i$未确定的边数,那么期望要减去:

$$\frac{p[i]\times (p[i]-1)}{2}+p[i]\times \frac{q[i]}{2}+\frac{q[i]\times (1-q[i])}{2}\times \frac{1}{4}$$

时间复杂度:$\Theta(n+e)$。

期望得分:$100$分。

实际得分:$100$分。


代码时刻

#include<bits/stdc++.h>
using namespace std;
const int mod=1000000007;
const int two=500000004;
const int six=166666668;
const int eig=125000001;
int n,e;
int du[100001],p[100001],q[100001];
long long ans;
long long qpow(long long x,long long y)
{
	long long res=1;
	while(y)
	{
		if(y&1)res=res*x%mod;
		x=x*x%mod;
		y>>=1;
	}
	return res;
}
int main()
{
	scanf("%d%d",&n,&e);
	while(e--)
	{
		int x,y;
		scanf("%d%d",&x,&y);
		du[x]++;du[y]++;p[x]++;
	}
	ans=1LL*n*(n-1)%mod*(n-2)%mod*six%mod;
	for(int i=1;i<=n;i++)
	{
		q[i]=n-du[i]-1;
		ans=(ans-1LL*p[i]*(p[i]-1)%mod*two%mod-1LL*p[i]*q[i]%mod*two%mod-1LL*q[i]*(q[i]-1)%mod*eig%mod+mod)%mod;
	}
	printf("%lld\n",ans);
	return 0;
}

rp++

标签
易学教程内所有资源均来自网络或用户发布的内容,如有违反法律规定的内容欢迎反馈
该文章没有解决你所遇到的问题?点击提问,说说你的问题,让更多的人一起探讨吧!