BZOJ2698: 染色

求出m次染色后被染成白色的格子个数的期望
相当于求每个格子被染白的概率，每个格子带权乘上他的概率，和就是期望
求一个格子m次被染白的概率，可以求1次被染白的概率p，$1-(1-p)^m$（1-m次都没被染的概率）
求一次被染白的概率，就是求有多少个不同的区间覆盖了这个格子，因为区间的总的种类数是易得的

那么有多少个区间覆盖了一个格子，是可以用DP求的（好像有O(1)的推法，一个公式什么的？我不会，有会的评论教一下我qaq？）

我做这题的时候zz了，写了一个树状数组求，具体说，就是走到一个点的时候，枚举以他为左端点的 k <script type="math/tex" id="MathJax-Element-54">k</script>个区间，对于每个区间在他的左端点加一个+1标记，右端点+1的地方加-1标记，标记的前缀和就是这个点被多少个区间覆盖，然后发现，+1标记打在同一个地方可以直接加，-1标记是区间加，可以差分后用树状数组维护（树状数组小常数过去了，线段树不知道行不行）

code：

#include<set>
#include<map>
#include<deque>
#include<queue>
#include<stack>
#include<cmath>
#include<ctime>
#include<bitset>
#include<string>
#include<vector>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<climits>
#include<complex>
#include<iostream>
#include<algorithm>
#define ll long long
#define lowbit(x) x&(-x)
using namespace std;

const int maxn = 1100000;

int n,m,s,t,c[maxn];

long double pw(long double t,int k)
{
    long double re=1.0,tmp=t;
    for(;k;k>>=1,tmp=tmp*tmp)
        if(k&1) re*=tmp;
    return re;
}

void upd(int x,const int v){for(;x<=n;x+=lowbit(x)) c[x]+=v;}
int query(int x){int re=0;for(;x;x-=lowbit(x))re+=c[x];return re;}

int main()
{
    scanf("%d%d%d%d",&n,&m,&s,&t); int u=n-s+1;
    ll num=0,tmp=0;
    for(int i=1;i<=u;i++)
    {
        int r=min(n,i+t-1);
        num+=(ll)(r-i-s+2);
    }

    if(num==0) { printf("0.000\n"); return 0; }
    long double ans=0.0,al=num;
    for(int i=1;i<=n;i++)
    {
        int l=i+s-1,r=min(n,i+t-1);
        if(l<=n) tmp+=(ll)(r-l+1); 
        upd(i+s,-1); upd(i+t+1,1);
        tmp+=(ll)query(i);
        long double k=tmp/al;
        ans+=1.0-pw(1-k,m);
    }
    printf("%.3Lf\n",ans);

    return 0;
}