【概率DP】小象涂色

【题意】
n个箱子编号为1至n，m种颜色编号为0至m-1，涂色k次，每个箱子初始颜色为1。
每次涂色给出一个区间（l，r）,其中每个箱子都有50%几率被涂上随机一种颜色。涂色后箱子颜色变为（原颜色*新颜色%c）。
求箱子颜色和期望。
多组数据。
T<=10，n/k<=50，m<=100。
【思路】
鲁莽的看见期望就认为是期望，然后想着怎么倒着推，然后就没有然后了。
膜题解才知道这是一道概率题，期望=该情况概率*该情况权值。
f[i][j]表示一个箱子在i次涂色的区间内经过若干次涂色后颜色后变为j的概率。
我们又将这个箱子放进一次涂色的区间内，有以下情况：
它在涂色的区间内，却没有被涂色。几率是50%，转移到f[i+1][j]，f[i+1][j]+=50%*f[i][j]。
它被涂上x色，几率是0.5/c，转移到f[i+1][j*x%c]，f[i+1][j*x%c]+=0.5/c*f[i][j]。
【代码】

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#define MAXN 60
#define MAXM 110
using namespace std;
inline int read()
{
    int ans=0,f=1;char x=getchar();
    while(x<'0'||x>'9')
    {
        if(x=='-')f=-1;
        x=getchar();
    }
    while(x>='0'&&x<='9')
    {
        ans=ans*10+x-'0';
        x=getchar();
    }
    return ans*f;
}
int n,m,k;
int a[MAXN];
double f[MAXN][MAXM];
void ex(int x,int y)
{
    f[x+1][y]+=0.5*f[x][y];
    for(int i=0;i<m;i++)
        f[x+1][y*i%m]+=f[x][y]*0.5/m;
    return ;
}
void dp()
{
    memset(f,0,sizeof(f));
    f[0][1] = 1.0;
    ex(0,1);
    for(int i=1;i<=k;i++)
        for(int j=0;j<m;j++)
            ex(i,j);
    return ;
}
int main()
{
    int T=read();
    while(T--)
    {
        n=read();m=read();k=read();
        dp();
        memset(a,0,sizeof(a));
        for(int i=1;i<=k;i++)
        {
            int l=read(),r=read();
            a[l]++;a[r+1]--;
        }
        for(int i=1;i<=n;i++)a[i]+=a[i-1];
        double ans=0.0;
        for(int i=1;i<=n;i++)
            for(int j=0;j<m;j++)
                ans+=f[a[i]][j]*(double)j;
        printf("%.9lf\n",ans);
    }
    return 0;
}