[AGC032F]One Third

Description

一个长度为1的环,随机角度切n刀,对于每一块长度为x的,取最接近1/3的,即|x-1/3|最小的那一块
问这个最小值的期望
n<=10^6

Solution

感觉还是有一些不明觉厉
首先,对于一次在x的划分,在x画红线,在x+π/3画蓝线,在x-π/3画绿线,问题变成任意两条异色线之间的夹角的最小值
由于所有线是对称?的所以可以只看某一段
也就是说,令第一次划分的x为0,x+π/3为1/3,问题变成在(0,1/3)中撒n-1个随机颜色得随机点,其中相邻异色点距离的最小值的期望
先不考虑异色,问题是在(0,1)中撒n-1个点,距离最小值的期望
考虑概率函数P(mn>=x),每一段都>=x的话我们可以从每一段中挖去x,剩余的点在(1-nx)中随机,所以P(mn>=x)=(1-nx)^(n-1)
积分一下可以得到最小值的期望为1/n^2
考虑第k小的期望E[k],E[k]-E[k-1]可以把所有<k的线段删去并把所有>=k的线段减去E[k-1]
E[k]E[k1]=1(Nk+1)2(1i=1k1E[i]E[k1](Nk+1))E[k]-E[k-1]={1\over (N-k+1)^2}(1-\sum_{i=1}^{k-1}E[i]-E[k-1]*(N-k+1))
归纳可得E[k]E[k1]=1N(Nk+1)E[k]-E[k-1]={1\over N(N-k+1)}
设P1(k)表示至少前k段颜色相同的概率,显然P1(k)=(1/3)^k
那么答案为i=1np1(i1)(E[i]E[i1])=i=1n13i1N(Ni+1)\sum_{i=1}^{n}p1(i-1)(E[i]-E[i-1])=\sum_{i=1}^{n}{1\over 3^{i-1}N(N-i+1)}
由于我们是当总长为1算所以最后还需要再乘以1/3

Code

#include <cstdio>
#include <cstring>
#include <algorithm>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
using namespace std;

typedef long long ll;

const int N=1e6+5,Mo=1e9+7,inv3=(Mo+1)/3;

int n,inv[N];

int main() {
	scanf("%d",&n);
	inv[0]=inv[1]=1;fo(i,2,n) inv[i]=(ll)(Mo-Mo/i)*inv[Mo%i]%Mo;
	int ans=0;
	for(int pw=inv3,i=1;i<=n;i++,pw=(ll)pw*inv3%Mo) (ans+=(ll)inv[n-i+1]*pw%Mo)%=Mo;
	printf("%d\n",(ll)ans*inv[n]%Mo);
	return 0;
}
©️2020 CSDN 皮肤主题: 编程工作室 设计师: CSDN官方博客 返回首页
实付0元
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值