计蒜客 And and Pair 题解

题目链接

分析：

（将一个二进制数的每一位$1$的下标作元素构成一个集合）
$i$&$n$=$i$，那么就是说i中的$1$是$n$中的子集。$i$的个数即为$n$的子集的个数。$i$&$j$=$0$，那么说$i$与$j$没有共同位上的$1$，即交集为空集。故对于一个确定的$i>1$，j的个数有$2^m$个($m$：$i$中最高位$1$后的$0$的个数)

所以可以想到从后往前枚举，即从小到大枚举i的最高位的位数，例如$n$=$11010$
先枚举$i$=$1$ $0$，再枚举$i$=$1$ $0$ $\_$ $0$的所有情况（ $\_$ 表示可以填入$0$和$1$的空位），然后是$i$= $1$ $\_$ $0$ $\_$ $0$。
分析到了这里就很容易想到动态规划的方法了。定义$dp[u]$为第$u$个枚举状态（即从后往前i的最高位取第u个1）时的$ij$pair的数目，那么可得：
$$dp\left[u\right]=\sum _{v=1}^{u-1}dp\left[v\right]\cdot 2^{k_{u,v}-1}+2^l$$
$k_{u,v}$：从第$v$个$1$到第$u$个$1$间的下标差
$l$ ：第$u$个$1$后面的二进制数的个数

公式的具体解释：
以上面的 $n=11010$ 为例,$i$= $1$ $\_$ $0$ $\_$ $0$（$u=3$）可以分解为
$i$= $1$ [$1$ $0$ $\_$ $0$] ($u=2$)
和
$i$= $1$ $0$ $0$ [$1$ $0$] ($u=1$)
和
$i$= $1$ $0$ $0$ $0$ $0$
其中[]中的部分分别是$u=2,1$时的情况,这样我们就把$u=3$的情况分成了$u=2,1$的情况和最高位全是0的情况。而$u=3$和$u=1$中间相差了两个0,这两个0在$u=1$的基础上又给j带来了两个不确定位，所以加上的是$dp[1]\cdot 2^2$,这也就得到了公式的第一部分${\sum }_{v=1}^{u-1}dp\left[v\right]\cdot 2^{k-1}$。而第二部分很简单，就是最高位全是0的情况。
虽然得到了公式，但是这个公式的带来的时间复杂度（最坏为$O(|S{|}^2)$）仍然不理想。不过容易发现是可以转化成为递推式的。
$$dp\left[u+1\right]=\sum _{v=1}^{u}dp\left[v\right]\cdot 2^{k_{u+1,v}-1}+2^{l+k_{u+1,u}}$$

$$=(\sum _{v=1}^{u}dp\left[v\right]\cdot 2^{k_{u,v}-1}+2^l)\cdot 2^{k_{u+1,u}}$$

$$=(\sum _{v=1}^{u-1}dp\left[v\right]\cdot 2^{k_{u,v}-1}+2^l+dp\left[u\right]\cdot 2^{-1})\cdot 2^{k_{u+1,u}}$$

$$=(dp\left[u\right]+dp\left[u\right]\cdot 2^{-1})\cdot 2^{k_{u+1,u}}$$

$$=3\cdot 2^{k_{u+1,u}-1}\cdot dp[u]$$
第一次用LaTex，等号不会对齐[捂脸]
这下我们就得到了简洁的递推式 :

$$dp\left[u+1\right]=3\cdot 2^{k_{u+1,u}-1}\cdot dp[u]$$

时间复杂度降低到了$O(|S|)$。不计字符串的存储空间，空间复杂度为$O(1)$。
注意：上面我们求的情况中没有包含$i=0$的情况，所以记得最后将计算结果加 $1$

AC代码：

#include<iostream>
#include<string>
using namespace std;
const int MAX_N=1e5+10;
const int MOD=1e9+7;
typedef long long ll;
int main()
{
	ll sum=0,now=0,last=0;
	string s;
	int t;
	cin>>t;
	while(t--)
	{
		cin>>s;
		int i;
		for(i=s.size()-1,now=1;i>=0&&s[i]!='1';i--,now=(now*2)%MOD);//从后往前找到第一个1
		if(i<0){cout<<1<<'\n';continue;}//n=0
		sum=(now+1)%MOD;
        last=now;
        now=1;
		for(i--;i>=0;i--)
			if(s[i]=='1')
			{
				last=(now*3*last)%MOD;
				sum=(last+sum)%MOD;
				now=1;
			}
			else now=(now*2)%MOD;
		cout<<sum%MOD<<'\n';
	}
	return 0;
}