熊猫（i）

题目描述

熊猫喜欢吃数，熊猫对与每个数都有他独特的评价。具体来说，熊猫对数 xx 的评价是个四元组 (a, b, c, d)(a,b,c,d)，计算方式如下：

首先将 xx 写成二进制形式（不含前导零），然后将这作为一个字符串，那么 a =a= 子串 “00”“00” 的数量，b =b= 子串 “01”“01” 的数量，c =c= 子串 “10”“10” 的数量，d =d= 子串 “11”“11” 的数量。

现在熊猫想吃掉区间 [A,B] 之间所有评价是 (a, b, c, d)(a,b,c,d) 的数，由于你被熊猫抓了起来，所以他要你回答他能吃到多少个数。如果你回答不出，他就只好吃你。

输入格式

第一行一个 0101 串为 AA 的二进制表示。

第一行一个 0101 串为 BB 的二进制表示。

第三行四个数 a, b, c, da,b,c,d。

输出格式

一个数表示答案，对 10000000071000000007 取模。

样例

样例输入 1

10
1001
0 0 1 1

样例输出 1

1

数据范围与提示

子任务一（13pts）：B ≤ 2^{20}B≤220。

子任务二（15pts）：B ≤ 2^{100}B≤2100。

子任务三（19pts）：B ≤ 2^{1000}B≤21000。

子任务四（15pts）：A = 1, B = 2^k − 1A=1,B=2k−1。

子任务五（38pts）：无特殊限制。

对于所有的数据，1 ≤ A ≤ B ≤ 2^{100000}, a + b + c + d ≤ 10^51≤A≤B≤2100000,a+b+c+d≤105。

来源

2019年1月雅礼集训

 

---

solution

我们可以用abcd算出区间长，需要多少0（s0），这些0分成几段（n0），需要多少1（s1），这些1分成几段（n1），

我们先不考虑限制。

考虑插板法，0的方案为C（s0-1,n0-1）

方案数C（s0-1,n0-1）*C（s1-1,n1-1）

有限制的就从第1为开始像数位dp一样记上下界有没有顶满。

如果都没顶满就直接算。

这样效率就2n了

#include<cstdio>
#include<iostream>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<cmath>
#define maxn 100005
#define mod 1000000007
#define ll long long
using namespace std;

int n,la,lb,a,b,c,d,n0,n1,s0,s1;
char A[maxn],B[maxn];
ll h[maxn],ny[maxn],ans;
ll work(ll a,int num){
    ll ans=;
    while(num){
        if(num&)ans=ans*a;
        a=a*a;a%=mod;ans%=mod;num>>=;
    }
    return ans;
}
void pre(){
    n=;

    h[]=;for(int i=;i<=n;i++)h[i]=h[i-]*i%mod;
    ny[n]=work(h[n],mod-);
    for(int i=n-;i>=;i--){
        ny[i]=ny[i+]*(i+)%mod;
    }
}
void Pd(){
    if(n>la){
        A[]='';for(int i=;i<=n;i++)A[i]='';
    }
    if(n<la){puts("");exit();}
    if(n>lb){puts("");exit();}
    if(n<lb){
        for(int i=;i<=n;i++)B[i]='';
    }
}
ll C(int N,int M){
    if(N<||M<||N<M)return ;
    return h[N]*ny[M]%mod*ny[N-M]%mod;
}
void calc(int n0,int s0,int n1,int s1){
    ll w0=C(s0-,n0-),w1=C(s1-,n1-);
    if(s0==&&n0==)ans+=w1;
    if(s1==&&n1==)ans+=w0;
    ans=ans+w0*w1%mod;ans%=mod;
}
void solve(int i,int la,bool fa,bool fb,int n0,int s0,int n1,int s1){
    if((n0==&&s0>)||(n1==&&s1>))return;
    if(i==n+){
        if(n0==&&la==&&!n1&&!s1&&!s0)ans++;
        if(n1==&&la==&&!n0&&!s0&&!s1)ans++;
        ans%=mod;
        return;
    }
    if(!fa&&!fb){
        calc(n0,s0,n1,s1);
        if(la==)calc(n0-,s0,n1,s1);
        if(la==)calc(n0,s0,n1-,s1);
        return;
    }

    if(s0<||n0<||s1<||n1<)return;
    if(i==)solve(i+,,fa,fb,n0,s0,n1,s1-);
    else {

        if(!fa||A[i]==''){
            if(la==)solve(i+,,fa,fb&(B[i]==''),n0,s0-,n1-,s1);
            else solve(i+,,fa,fb&(B[i]==''),n0,s0-,n1,s1);
        }
        if(!fb||B[i]==''){
            if(la==)solve(i+,,fa&(A[i]==''),fb,n0-,s0,n1,s1-);
            else solve(i+,,fa&(A[i]==''),fb,n0,s0,n1,s1-);
        }
    }
}

int main()
{
    pre();
    scanf(" %s",A+);la=strlen(A+);
    scanf(" %s",B+);lb=strlen(B+);

    scanf("%d%d%d%d",&a,&b,&c,&d);
    n=a+b+c+d+;
    Pd();ans=;
    if(b==c){
        n0=b;n1=n0+;
        s0=a+n0,s1=d+n1;
        solve(,,,,n0,s0,n1,s1);
        printf("%lld\n",ans);
    }
    else if(b+==c){
        n0=n1=c;
        s0=a+n0,s1=d+n1;
        solve(,,,,n0,s0,n1,s1);
        printf("%lld\n",ans);
    }
    else puts("");
    return ;
}