【gym102394B】Binary Numbers（DP）

题意：From https://blog.csdn.net/m0_37809890/article/details/102886956

思路：

可以发现转移就是右上角的一个区间前缀和

std只要开1倍空间，个人习惯开2倍滚动数组， 1倍有时候没想清楚的前提下初始化或者求前缀和顺序手抖就会gg，反正时空只有常数倍差距

 #include<bits/stdc++.h>
 using namespace std;
 typedef long long ll;
 typedef unsigned int uint;
 typedef unsigned long long ull;
 typedef long double ld;
 typedef pair<int,int> PII;
 typedef pair<ll,ll> Pll;
 typedef vector<int> VI;
 typedef vector<PII> VII;
 typedef pair<ll,ll>P;
 #define N  1<<20
 #define M  1000000
 #define INF 1e9
 #define fi first
 #define se second
 #define MP make_pair
 #define pb push_back
 #define pi acos(-1)
 #define mem(a,b) memset(a,b,sizeof(a))
 #define rep(i,a,b) for(int i=(int)a;i<=(int)b;i++)
 #define per(i,a,b) for(int i=(int)a;i>=(int)b;i--)
 #define lowbit(x) x&(-x)
 #define Rand (rand()*(1<<16)+rand())
 #define id(x) ((x)<=B?(x):m-n/(x)+1)
 #define ls p<<1
 #define rs p<<1|1
 #define fors(i) for(auto i:e[x]) if(i!=p)

 const int MOD=1e8+,inv2=(MOD+)/;
       double eps=1e-;
       int dx[]={-,,,};
       int dy[]={,,-,};

 ll s[][][],dp[][][];
 int l[N],r[N],m,n;

 int read()
 {
    int v=,f=;
    char c=getchar();
    while(c<||<c) {if(c=='-') f=-; c=getchar();}
    while(<=c&&c<=) v=(v<<)+v+v+c-,c=getchar();
    return v*f;
 }

 ll readll()
 {
    ll v=,f=;
    char c=getchar();
    while(c<||<c) {if(c=='-') f=-; c=getchar();}
    while(<=c&&c<=) v=(v<<)+v+v+c-,c=getchar();
    return v*f;
 }

 void add(ll &a,ll b)
 {
     a+=b;
     if(a>=MOD) a-=MOD;
     if(a<) a+=MOD;
 }

 int lcp(int a,int b)
 {
     int res=;
     per(i,m-,)
       if((a>>i&)==(b>>i&)) res++;
        else break;
     return res;
 }

 void solve()
 {
     m=read(),n=read();
     rep(i,,n) l[i]=read(),r[i]=read();
     rep(v,,)
      rep(j,,m+)
       rep(k,,m+) s[v][j][k]=dp[v][j][k]=;
     int v=;
     dp[v][][m]=s[v][][m]=;
     rep(i,,m)
      per(j,m,)
      {
          add(s[v][i][j],s[v][i][j+]);
          if(i)
          {
             add(s[v][i][j],s[v][i-][j]);
             add(s[v][i][j],-s[v][i-][j+]);
          }
      }
      int A,B,C,D;
      rep(i,,n)
      {
          v^=;
          rep(j,,m)
           rep(k,,m) s[v][j][k]=dp[v][j][k]=;
          rep(j,l[i],r[i])
          {
             A=lcp(j,l[i]),B=lcp(j,r[i]);
             if(i<n) C=lcp(j,l[i+]);
              else C=;
             if(i>) D=lcp(j,r[i-]);
              else D=;
             dp[v][C][B]=(dp[v][C][B]+s[v^][A][D]*j)%MOD;
          }
          rep(i,,m)
           per(j,m,)
           {
                 s[v][i][j]=dp[v][i][j];
                 add(s[v][i][j],s[v][i][j+]);
                 if(i)
                 {
                     add(s[v][i][j],s[v][i-][j]);
                     add(s[v][i][j],-s[v][i-][j+]);
                 }
           }
      }
      printf("%I64d\n",s[v][m][]);
 }

 int main()
 {
     //freopen("1.in","r",stdin);
     //freopen("1.out","w",stdout);
     int cas=read();
     while(cas--) solve();
     return ;
 }