洛谷3600，大力期望DP

传送门
这题从下午做到晚上，真是欲仙欲死。
题解区里有题解，就不说了，一开始有三个地方写挂。
1.把包含于被包含的关系想反了;
2.没有考虑一个元素不在任何询问里的情况
3.处理2问题我用的是排序，而且cmp是return l<x.l;
感觉我好菜啊。

#include<cstdio>
#include<cstring> 
#include<algorithm>
typedef long long ll;
const int N=2010;
const ll mo=666623333;
int n,m,i,j,f[N],x,q,k,p,ans,g[N],s,s1,s2,pp[N],inp[N];
struct seg{
    int l,r;
    bool operator<(const seg&x)const{
        return l<x.l || (l==x.l && r<x.r);
    }
}a[N],b[N];
bool bb[N];
ll pow(ll x,int y){
    if(!y)return 1;
    ll u=pow(x,y>>1);
    return y&1?u*u%mo*x%mo:u*u%mo; 
}
int main(){
    //freopen("a.txt","r",stdin);
    scanf("%d%d%d",&n,&x,&q);
    for(i=1;i<=q;++i){
        scanf("%d%d",&a[i].l,&a[i].r);
        for(j=1;j<i;++j)
            if(a[j].l<=a[i].l && a[i].r<=a[j].r)bb[j]=1;
                else if(a[i].l<=a[j].l && a[j].r<=a[i].r)bb[i]=1;
    }
    for(i=1;i<=q;++i)if(!bb[i])a[++m]=a[i];
    std::sort(a+1,a+m+1);
    for(i=1;i<=n;++i)
        for(j=1;j<=m;++j)
            if(a[j].l<=i && i<=a[j].r)b[i].l?b[i].r=j:b[i].l=b[i].r=j;
    std::sort(b+1,b+n+1);
    for(i=1;i<=n && !b[i].l;++i);
    for(j=i;j<=n;++j)b[j-i+1]=b[j];
    n-=i-1;
    for(i=*inp=*pp=1;i<x;++i){
        p=1ll*i*pow(x,mo-2)%mo;
        for(pp[j=1]=1+mo-p;j<=n;++j)pp[j]=1ll*pp[1]*pp[j-1]%mo;
        for(inp[j=1]=pow(pp[1],mo-2);j<=n;++j)inp[j]=1ll*inp[j-1]*inp[1]%mo;
        memset(g,0,sizeof g);
        for(k=s=0,j=1;j<=n;++j){
            while(b[k].r<b[j].l-1)s=(s+mo-g[k++])%mo;
            f[j]=1ll*p*pp[j-1]%mo*(s+(b[j].l==1))%mo;
            s=(s+(g[j]=1ll*f[j]*inp[j]%mo))%mo;
        }
        for(s1=0,j=1;j<=n;++j)
            if(b[j].r==m)s1=(s1+1ll*f[j]*pp[n-j])%mo;
        ans=(ans+1ll*(s1+mo-s2)*i)%mo;
        s2=s1;
    }
    printf("%lld\n",(ans+1ll*(1+mo-s2)*x%mo)%mo);
    return 0;
}