期望 【bzoj2698】 染色

题目大意：
有n个格子，初始都是黑色的。
有M次操作，每次随机挑选一个长度范围在s到t之间的区间进行染色，把选定区间内的格子都染白。
询问染白格子数量的期望。
1 ≤ S ≤ T ≤ N ≤ 1000000，0 ≤ M ≤ 1000000

题目分析：
因为期望具有线性性，所以我们把每个格子为白色的期望（此题中就相当誉概率）加在一起就是答案。
但是一个格子被染白的概率不好求，所以我们转化成不被染白的概率。
一个格子不被染白的概率等于所有染不到它的方案数除以总的方案数。

我们选择用一个a数组来记录，a[i]表示我们有i个格子，我们能选出的不同的区间个数是多少，这个O（n）即可算出，其实也可以找规律套公式O（1）算出来，但是这样思维量比较小，而且时间复杂度可以接受。

一个格子不被染白的概率就是在它左边的部分选区间的总数加上在它右边部分选区间的总数除以全部格子内选区间的总数。（就是在它左边随便选，在它右边随便选，但是不能跨过它）

这样我们得到了一个格子一次染色中不被染白的概率，m次操作就是它的m次方，快速幂就好了。
然后用1减去我们得到的答案，就是这个格子被染白的概率。
每个格子都这样算，然后把答案加在一起输出即可。

代码如下：

#include<cstdio>
#define N 1200000
using namespace std;
int n,m,s,t;
double a[N];
double ans,cs,k;
int main()
{
    scanf("%d%d%d%d",&n,&m,&s,&t);
    for(int i=s;i<=t;i++)   a[i]=a[i-1]+i-s+1;
    for(int i=t+1;i<=n;i++) a[i]=a[i-1]+t-s+1;
    for(int i=1;i<=n;i++)
    {
        k=(a[i-1]+a[n-i])/a[n];
        cs=1;
        int p=m;
        while(p>0)
        {
            if(p&1) cs*=k;
            k*=k;
            p>>=1;
        }
        ans+=1-cs;
    }
    printf("%.3lf",ans);
    return 0;
}