yjqaa--BJ80中集训

给定一个非负整数x,你每次可以对x执行两种操作，要么令x=x−1，要么令x=x−2^i（前提条件是x&2^i≠0。

令f(x,y)表示从x变成y的最小操作次数。

求1<=x<=n,1<=y<=n,的所有f(x,y)的和对1000000007取模的结果。

输入格式

一行一个另一串表示询问的N的二进制表示（保证合法，没有前导0）

输出格式

一行一个整数表示答案对1000000007取模的结果。

样例输入

Copy

111

样例输出

Copy

43

数据范围

对于20%的数据，N的二进制表示不超过1010位。 对于50%的数据，N的二进制表示不超过20位。 对于100%的数据，N的二进制表示不超过500位。

这道题考试的时候没有抓住精髓，在错误的道路上越走越远

考试时一直在想对于每一个i如何快速求出所有j的贡献

但是实际上我们可以两个序列一起考虑

可以发现有一个贪心的变化策略，对于确定的i,j，是可以进行特定的转化的

将整个序列分为三段，分别表示需不需要使用-1，其中有一种是中间的过渡状态

整道题就是一道非常简单的数位dp了

#include<bits/stdc++.h>
using namespace std;
const int mod=1e9+7;
char s[505];
int dp[505][1005][2][2][3],n,temp1[505],temp2[505];
int moc(int x)
{
	if(x>=mod) return x-mod;
	return x;
}
int dfs(int depth,int top1,int top2,int kind,int val)
{
	if(depth==n+1)
	{
		if(kind==0||kind==2) return val;
		return 0;
	}
	if(dp[depth][val][top1][top2][kind]!=-1) return dp[depth][val][top1][top2][kind];
	int maxa1,maxa2;
	if(top1) maxa1=s[depth]-'0';
	else maxa1=1;
	int ans=0;
	for(int i=0;i<=maxa1;i++)
	{
		for(int j=0;j<=1;j++)
		{
			if((top2||!kind)&&j>i) continue;
			if(kind==1&&i>j) continue;
			temp1[depth]=i;
			temp2[depth]=j;
			if(kind==0&&i>j) val++;
			if((kind==2||kind==1)&&i==1) val++;
			if((kind==2||kind==1)&&j==0) val++;
			if(i>j&&!kind)
			ans=moc(ans+dfs(depth+1,top1&(i==maxa1),top2&(j==i),1,val));
			if(i<j&&kind==1)
			ans=moc(ans+dfs(depth+1,top1&(i==maxa1),top2&(j==i),2,val));
			if(!(i<j&&kind==1))
			ans=moc(ans+dfs(depth+1,top1&(i==maxa1),top2&(j==i),kind,val));
			if((kind==2||kind==1)&&i==1) val--;
			if((kind==2||kind==1)&&j==0) val--;
			if(i>j&&!kind) val--;
		}
	}
	return dp[depth][val][top1][top2][kind]=ans;
}
int main()
{
	scanf("%s",s+1);
	n=strlen(s+1);
	memset(dp,-1,sizeof(dp));
	cout<<dfs(1,1,1,0,0);
	return 0;
} 
/*
111
*/