https://www.luogu.org/recordnew/lists?uid=56840
题解
首先可以发现我们对于每种颜色的鱼，长一点的能够覆盖的方案已定完全包含短一点的方案。
所以我们可以只对每种颜色最长的鱼计算贡献。
然后有一个(naive)的想法，我们从按照最长的鱼的长度小到大枚举每种颜色，然后算出那条最长的鱼能够包含的方案。
这样会算重。
那么我们还有一个(naive)的想法，我们可以在枚举的时候，只维护出比在a这种颜色前面的颜色的所有方案。
这样会算少。
考虑在什么情况吗，没有被算到。
对于两种颜色(a,b)，我们本该在枚举a的时候枚举(a,b)这种集合，结果由于排序所以没有枚举到，但在枚举b的时候因为长度原因没有够到a。
所以我们在枚举a的时候，当我们没有拿走能够拿走的所有a的时候可以按照上面的方法做，但是如果a的全部拿走了，我们需要把所有颜色b算出来，这些b可以拿到的颜色个数是和a一样的。
这样的话，a的数量是now+1的，枚举到b是最多枚举到now，所以不会算重。
我也不知道我在写啥
代码
#include<bits/stdc++.h>
#define N 500009
using namespace std;
typedef long long ll;
int n,k,mod,cnt[N];
int id[N],pos[N],ct[N];
int now[N],nxt[N];
int mx[N],_mx[N];
ll ans;
inline ll rd(){
ll x=0;char c=getchar();bool f=0;
while(!isdigit(c)){if(c=='-')f=1;c=getchar();}
while(isdigit(c)){x=(x<<1)+(x<<3)+(c^48);c=getchar();}
return f?-x:x;
}
inline void MOD(ll &x){x=x>=mod?x-mod:x;}
inline bool cmp(int a,int b){return mx[a]<mx[b];}
struct node{
int len,co;
inline bool operator <(const node &b)const{
return len<b.len;
}
}a[N];
struct segment{
int tr[N<<2];
void upd(int cnt,int l,int r,int x){
if(l==r){tr[cnt]++;return;}
int mid=(l+r)>>1;
if(mid>=x)upd(cnt<<1,l,mid,x);
else upd(cnt<<1|1,mid+1,r,x);
tr[cnt]=1ll*tr[cnt<<1]*tr[cnt<<1|1]%mod;
}
void build(int cnt,int l,int r){
tr[cnt]=1;
if(l==r)return;
int mid=(l+r)>>1;
build(cnt<<1,l,mid);build(cnt<<1|1,mid+1,r);
}
ll query(int cnt,int l,int r,int L,int R){
if(L>R)return 1;
if(l>=L&&r<=R)return tr[cnt];
int mid=(l+r)>>1;
if(mid>=L&&mid<R)return query(cnt<<1,l,mid,L,R)*query(cnt<<1|1,mid+1,r,L,R)%mod;
if(mid>=L)return query(cnt<<1,l,mid,L,R);
if(mid<R)return query(cnt<<1|1,mid+1,r,L,R);
}
inline void init(int n){build(1,1,n);}
}T;
inline int efs(int num,int l,int r){
int ans=l;
while(l<=r){
int mid=(l+r)>>1;
if(_mx[mid]<num)ans=mid,l=mid+1;
else r=mid-1;
}
return ans;
}
int main(){
n=rd();
k=rd();
mod=rd();
for(int i=1;i<=n;++i){
a[i].len=rd();a[i].co=rd();
mx[a[i].co]=max(mx[a[i].co],a[i].len);
}
sort(a+1,a+n+1);
for(int i=1;i<=k;++i)id[i]=i;
sort(id+1,id+k+1,cmp);
for(int i=1;i<=k;++i)pos[id[i]]=i,_mx[i]=mx[id[i]];
for(int i=n;i>=1;--i){
a[i].co=pos[a[i].co];
cnt[a[i].co]++;
nxt[i]=now[a[i].co];
now[a[i].co]=i;
}
T.init(k);
int p=1;
for(int i=1;i<=n;++i){
cnt[a[i].co]--;
while(a[p].len*2<=a[i].len)now[a[p].co]=nxt[now[a[p].co]],T.upd(1,1,k,a[p].co),ct[a[p].co]++,p++;
if(!cnt[a[i].co]){
int ps=efs(a[now[a[i].co]].len*2,a[i].co,k);
ll x=T.query(1,1,k,1,a[i].co-1);
MOD(ans+=x*ct[a[i].co]%mod);
MOD(ans+=x*T.query(1,1,k,a[i].co+1,ps)%mod);
}
}
cout<<ans;
return 0;
}