LOJ#2331. 「清華集訓 2017」某位歌姬的故事

將序列離散化後，可以給每個點確定一個取值的上界wi$w_i$

對於限制(lj,rj,cj)$(l_j,r_j,c_j)$ ，只有[lj,rj]$[l_j,r_j]$ 中wi=cj$w_i=c_j$ 的點能貢獻
對於一個cj$c_j$ ，將所有wi=cj$w_i=c_j$ 的點拿出來，令f[i][j]$f[i][j]$ 表示滿足了前i個區間，最後一個權值取到了wi$w_i$ 的點i$i$ 是j$j$ 的方案數，轉移用前綴和優化到O(1)

不同的cj$c_j$ 互相獨立，分別dp後把貢獻乘起來

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
using namespace std;

const int maxn = 2005;
const int mod  = 998244353;
inline void add(int &a,const int &b){a+=b;if(a>=mod)a-=mod;}

int pw(int x,int k)
{
    int re=1;
    for(;k;k>>=1,x=(ll)x*x%mod) if(k&1)
        re=(ll)re*x%mod;
    return re;
}
int inv(int x){ return pw(x,mod-2); }

int n,Q,A;

struct Op
{
    int l,r,c;
    friend inline bool operator <(const Op x,const Op y)
    {
        if(x.c==y.c) return x.l==y.l?x.r<y.r:x.l<y.l;
        else return x.c<y.c;
    }
}op[maxn];

int t[maxn],tp;
struct Point
{
    int l,r;
}p[maxn],np[maxn]; int cnt,N;
map<int,int>mp;
void Trans()
{
    mp.clear(); tp=0;
    for(int i=1;i<=Q;i++) t[++tp]=op[i].l,t[++tp]=op[i].r;
    sort(t+1,t+tp+1); t[0]=0; cnt=0;
    for(int i=1;i<=tp;i++)
    {
        if(t[i]-1>t[i-1]) p[++cnt]=(Point){t[i-1]+1,t[i]-1};
        if(i==1||t[i]!=t[i-1]) p[++cnt]=(Point){t[i],t[i]},mp[t[i]]=cnt;
    }
    if(p[cnt].r!=n) p[cnt+1]=(Point){p[cnt].r+1,n},++cnt;
    for(int i=1;i<=Q;i++) op[i].l=mp[op[i].l],op[i].r=mp[op[i].r];
}

struct node
{
    int pos,c;
    friend inline bool operator <(const node x,const node y){return x.pos<y.pos;}
}a[maxn]; int an;
multiset<int>S;
multiset<int>::iterator it;
int U[maxn];
void Checku()
{
    an=0;
    for(int i=1;i<=Q;i++) a[++an]=(node){op[i].l,op[i].c},a[++an]=(node){op[i].r+1,-op[i].c};
    sort(a+1,a+an+1);
    int nowa=1; S.clear();
    for(int i=1;i<=cnt;i++)
    {
        while(nowa<=an&&a[nowa].pos==i)
        {
            if(a[nowa].c<0) it=S.find(-a[nowa].c),S.erase(it);
            else S.insert(a[nowa].c);
            nowa++;
        }
        if(S.empty()) U[i]=A;
        else it=S.begin(),U[i]=(*it);
    }
}

int ok0[maxn],ok1[maxn],ok2[maxn];
int f[maxn],v[maxn];
int dp()
{
    memset(v,0,sizeof v);
    int ans=1;

    sort(op+1,op+Q+1);
    for(int i=1;i<=Q;)
    {
        int j; for(j=i+1;j<=Q&&op[j].c==op[i].c;j++);j--;
        int c=op[i].c;

        N=0; for(int k=1;k<=cnt;k++) if(U[k]==c) np[++N]=p[k],v[k]=1,f[N]=0;
        for(int k=1;k<=N;k++)
        {
            int siz=np[k].r-np[k].l+1;
            ok0[k]=pw(c-1,siz);
            ok2[k]=pw(c,siz);
            ok1[k]=(ok2[k]-ok0[k]+mod)%mod;
        }

        f[0]=1;
        for(;i<=j;i++)
        {
            int s=0,k;
            for(k=0;np[k].r<p[op[i].l].l&&k<=N;k++)
                s=(ll)s*ok2[k]%mod,add(s,f[k]),f[k]=0;
            for(;np[k].r<=p[op[i].r].r&&k<=N;k++)
            {
                add(f[k],(ll)s*ok1[k]%mod);
                s=(ll)s*ok0[k]%mod;
            }
            for(;k<=N;k++) f[k]=0;
        }
        int s=0;
        for(int k=N,temp=1;k>=0;k--)
        {
            add(s,(ll)f[k]*temp%mod);
            temp=(ll)temp*ok2[k]%mod;
        }
        ans=(ll)ans*s%mod;
    }

    for(int i=1;i<=cnt;i++) if(!v[i]) ans=(ll)ans*pw(A,p[i].r-p[i].l+1)%mod;
    return ans;
}

int main()
{
    int tcase; scanf("%d",&tcase);
    while(tcase--)
    {
        scanf("%d%d%d",&n,&Q,&A);
        for(int i=1;i<=Q;i++) scanf("%d%d%d",&op[i].l,&op[i].r,&op[i].c);
        Trans();
        Checku();

        printf("%d\n",dp());
    }

    return 0;
}