# 江西CCPC省赛-Rng(概率+逆元）

江西CCPC省赛-Rng(概率+逆元）

题意：

• 给出一个n，在[1,n]之间选一个R1，在[1,R1]之间选一个L1，得到区间[L1,R1],同理获取区间[L2,R2]，问两个区间相交的概率对1e9+7取模。

思路1:

• 分情况讨论，条件概率的使用
  
  AcCode

#include <bits/stdc++.h>
using namespace std;
#define fre freopen("C:\\Users\\22765\\Desktop\\in.txt","r",stdin);
#define ms(a) memset((a),0,sizeof(a))
#define re(i,a,b) for(int (i)=(a);(i)<(b);(i)++)
#define sf(x) scanf("%d",&(x))
#define rg register
#define il inline
typedef long long LL;
const int inf=(0x7f7f7f7f);
const int maxn=1e6+5;
const int mod=1e9+7;

LL inv[maxn];
LL sum_inv[maxn];

//线性递推求逆元
void init(){
    inv[1]=1;
    re(i,2,maxn){
        inv[i]=((mod-mod/i)*inv[mod%i])%mod;
    }   
}

//逆元前缀和处理
void getsum() {
    re(i,1,maxn){
        sum_inv[i]=sum_inv[i-1]+inv[i]%mod;
    }   
}

int main(){
    LL n;
    init();getsum();
//  for(int i=1;i<=20;i++)cout<<inv[i]<<" "<<sum_inv[i]<<endl;
    while(cin>>n){
        LL ans=(n+3)*inv[n]%mod*inv[4]%mod;
        re(i,1,n+1){
            ans=ans+i*(sum_inv[n]-sum_inv[i]+mod)%mod*inv[n]%mod*inv[n]%mod;
            ans%=mod;
        }
        cout<<ans<<endl;
    }
    return 0;
}

思路2：

• 直接求相交比较麻烦，可以通过先求不相交，再用1减去即可，正难则反定理，比上一种方法简单很多（逃
  
  其中分子为不相交时可能出现的位置，分母为两个端点在没有限制情况下所有可能出现的位置。

AcCode:

#include <bits/stdc++.h>
using namespace std;
#define fre freopen("C:\\Users\\22765\\Desktop\\in.txt","r",stdin);
#define ms(a) memset((a),0,sizeof(a))
#define re(i,a,b) for(int (i)=(a);(i)<(b);(i)++)
#define sf(x) scanf("%d",&(x))
#define rg register
#define il inline
typedef long long LL;
const int inf=(0x7f7f7f7f);
const int maxn=3000;
LL q=1000000007;
LL qpow(LL a,LL b){
    LL ans=1;
    a%=q;
    while(b){
        if(b&1)ans=ans*a%q;
        b>>=1;
        a=a*a%q;
    } 
    return ans;
}
LL inv(LL a){
    if(a==1)return 1;
    return qpow(a,q-2);
}

int main(){
    LL n;
    while(cin>>n){
        cout<<(n+1)*inv(2*n)%q<<endl; 
    }
    return 0;
}