Sum of Log——数位dp

题面：
问题的关键点在于[i&j=0]⌊log2 (i+j)+1⌋
其中i&j=0表明了i+j不可能出现任何进位的情况贡献答案，log2 (i+j)+1表示i和j最先出现1的位数。

很显然这是一个数位dp题

#include<bits/stdc++.h>
#define int long long
using namespace std;
const int mod = 1e9+7;
int dp[50][2][2][2];//长度，limit_x,limit_y,是否有前导0
int x[50],y[50];
int dfs(int len,int limx,int limy,int lead){
	if(len<0)return 1;//边界处理
	if(dp[len][limx][limy][lead]!=-1)return dp[len][limx][limy][lead];
	int upx = limx?x[len]:1;
	int upy = limy?y[len]:1;
	int ans = 0,cnt=0;
	for(int i=0;i<=upx;i++)
	{
		for(int j=0;j<=upy;j++)
		{
			if(i&j)continue;
			if(lead && (i||j)) cnt = len+1;//有前导0并且x or y在当前位有一个为1，那么这个1就是最早出现的答案
			else cnt=1;
			ans = (ans+dfs(len-1,limx&&i==x[len],limy&&j==y[len],lead&&!i&&!j)*cnt%mod)%mod;
		}
	}
	return dp[len][limx][limy][lead] = ans;
}
int solve(int xx,int yy)
{
	memset(dp,-1,sizeof dp);
	memset(x,0,sizeof x);
	memset(y,0,sizeof y);
	int len = 0;
	while(xx||yy){
		x[len] = xx&1;
		y[len] = yy&1;
		xx>>=1;
		yy>>=1;
		len++;
	}
	return dfs(len-1,1,1,1);
}
signed main()
{
	int T;
	cin>>T;
	while(T--)
	{
		int a,b;
		cin>>a>>b;
		cout<<(solve(a,b)-1+mod)%mod<<"\n";
	}
}

ps1：最初的时候以为dp[x的位数][y的位数]，但后来发现x，y的每一位是同时考虑的，那么直接取max(x的位数,y的位数)就行。
ps2：关于卡常的问题，这题的状态数是50222=400，转移数是22=4，询问次数为1e5，也就是计算次数约为40041e5=1.6*10^8，但因为实际情况往往跑不满，所以600+ms可过。