做完之后感觉更像是一道找规律的题。
思路:可以把一个数分成若干个2的i次方,比如一个数的二进制表示为1011,则这个数可以分成1000 + 10 + 1,然后求1011的种类数就可以从这三个数的种类数出发,
打表发现1 、 10 、100这些数的种类数都是3的i次方+2(i是二进制的位数)。
a + b = a ^ b就意味着a和b的二进制表示中不会有一位同时为1,所以可以从a + b的和入手,如果a + b = 0的话a + b只有一种情况,a + b = 1的话,a + b就有两种情况,a + b = 101的话,就有四种情况,所以我们可以发现情况数是2的j次方(j是二进制表示中1的个数)。
1000的计算过程中其实包含了10的计算过程,所以再往后计算10时其实就不是计算10、而是计算(1000的种类数) * (10的个数)了,所以每往后计算一个1就要多乘一个2
还有0建议最后再计算,不然过程中0会出现难以计算的重叠区域)。
#include<bits/stdc++.h>
#define IOS ios::sync_with_stdio(0);cin.tie(0);cout.tie(0);
#define endl '\n'
using namespace std;
typedef pair<int, int> PII;
typedef long long ll;
const int mod = 1e9 + 7, N = 100010;
ll w[N];
ll qmi(ll a, ll k)
{
ll res = 1;
while(k)
{
if(k & 1)res = res * a % mod;
k >>= 1;
a = a * a % mod;
}
return res;
}
int main()
{
IOS
for(int i = 1; i < N; i ++)
{
ll tmp = qmi(3, i - 1) + 1;
w[i] = tmp % mod;
}
string s;
cin >> s;
ll ans = 0, one = 0;
for(int i = 0; i < s.size(); i ++)
{
if(s[i] == '1')
{
int x = s.size() - i;
ans = (ans + w[x] * qmi(2, one) % mod ) % mod;
one ++;
}
}
cout << (ans + 1) % mod;
return 0;
}