【2020ccpc长春站】D题 Meaningless Sequence（dp，数学） 详解

题目链接：Meaningless Sequence

分析：
首先打表找规律：

#include <cstdio>
#include <iostream>
#include <string.h>
#include <algorithm>
using namespace std;

const int mod=1e9+7;

int main()
{
    int c=2;
    int a[105];
    a[0]=1;
    cout<<"a[0]="<<1<<endl;
    for(int i=1;i<=100;i++){
        a[i]=0;
        for(int j=0;j<i;j++)
              a[i]=max(a[i],a[i&j]);
        a[i]*=c;
        printf("a[%d]=%d\n",i,a[i]);
    }
    return 0;
}

可见$a_n$的值由其二进制表示中1的个数决定，设$a_n$的二进制表示中1的个数为x个，则$a_n$的值为$c^x$

然后考虑如何求${\sum }_{i=1}^n{a}_i$的和

先考虑下标0000~1111的所有数的和：

运用组合数学的思想：
（相当于从四个盒子中选x个放1的方案数）

①二进制表示下含0个1的数，值为$c^0$，有$C_4^0=1$个
②二进制表示下含1个1的数，值为$c^1$，有$C_4^1=4$个
③二进制表示下含2个1的数，值为$c^2$，有$C_4^2=6$个
④二进制表示下含3个1的数，值为$c^3$，有$C_4^3=4$个
⑤二进制表示下含4个1的数，值为$c^4$，有$C_4^4=1$个

则总的值为$C_4^0\cdot c^0+C_4^1\cdot c^1+C_4^2\cdot c^2+C_4^3\cdot c^3+C_4^4\cdot c^4$

很容易就能发现这个值就是$(1+c{)}^4$的二项展开式

综上，可以知道如果题目给的n二进制表示正好是1111111……是很容易求解的，答案就是$(1+c{)}^n$

但对于更一般的11101这种数应该怎么求解呢？

这么考虑：

假设我们已经知道了n为1101的答案，则在其前面加上一个1，也就是11101的答案是可以很快知道的:

因为所谓n为1101的答案，即下标i从0000 ~ 1101的$a_i$的和，而11101的答案，就是下标i从00000~11101的$a_i$的和

我们将00000~11101这些数划分为两部分，一部分为0000 ~ 1111，根据之前的推导值就是$(1+c{)}^n$，很容易求得

而另外一部分就是10000~11101，这部分数相当于0000 ~ 1101的每个数前面加个1推得的，每个数的二进制表示都多了一个1，那也就是每个数对应的$a_i$都乘上了1个c。因为之前假设我们已经知道了n为1101的答案，称其为x，则下标10000 ~11101的这些数就是$x\cdot c$

由此，我们得到了一个可以用于递推的表达式：
对于$n=(a_0{a}_1{a}_2\cdot \cdot \cdot a_k{)}_2$，如果我们知道了$m=(a_1{a}_2\cdot \cdot \cdot a_k{)}_2$的值x

则：

$n=x\cdot c+(1+c{)}^{k-1}$

代码：

从n的低位不断往高位推即可

#include <iostream>
#include <cstdio>
#include <algorithm>
using namespace std;

const long long mod=1e9+7;

string s;
long long c;


long long quickpow(long long base,long long power){
    long long ret=1;
    while(power){
        if(power%2)   ret=ret*base%mod;
        base=base*base%mod;
        power/=2;
    }
    return ret;
}


int main()
{
    ios::sync_with_stdio(false);
    cin>>s>>c;
    long long ans=1;
    for(int i=s.size()-1;i>=0;i--)
        if(s[i]=='1')
              ans=(ans*c%mod+quickpow(c+1,s.size()-i-1))%mod;
    cout<<ans<<endl;
    return 0;
}