[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[k-1]=\frac{1}{(N-k+1{)}^2}(1-{\sum }_{i=1}^{k-1}E[i]-E[k-1]\ast (N-k+1))$
归纳可得$E[k]-E[k-1]=\frac{1}{N(N-k+1)}$
设P1(k)表示至少前k段颜色相同的概率，显然P1(k)=(1/3)^k
那么答案为${\sum }_{i=1}^{n}p1(i-1)(E[i]-E[i-1])={\sum }_{i=1}^n\frac{1}{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;
}