数位dp

emmmm一些技巧？？？

先上题目！

windy数 https://www.luogu.org/problem/P2657

数位dp最基础的入门了呢。

重要技巧 （划重点嘻嘻嘻）

1.如果要求求l-r的区间立的话，转化成[1,r]-[1,l-1]，即前缀和的形式，因为这样会很方便qwq

2. dp的话要从最高位开始，原因是你要判断是否超过了边界（就是如果前i位都和边界相同，那么i+1位置就要）但作为一个懒人，我们一般存数的时候，都是从最低为开始存贮。所以为了方便，我们记忆化搜索的时候也可以从高到低搞（最后到了0就弄出去！）

3.个人建议把所有的要用的东西（就是dfs的所有参数）都存下来，因为这样写真的超级超级方便呐！（开心）

4.别的没啥了嘛。。。哦哦哦一般要用的limit（判断前i位是否已经到达饱和）还有注意前导0啊啊啊啊我死了

代码来啦（我知道你只关注这个www）

#include <bits/stdc++.h>

using namespace std;
typedef long long ll;
int s[30];
ll f[30][15][2][2];
ll dfs(int pos,int pre,int limit,int rt)
{
//	cout<<pos<<endl;
	if(pos==0)
	{
//		cout<<"here"<<endl;
		return 1;	
	}
	if(f[pos][pre][limit][rt]!=-1) 
	{
		return f[pos][pre][limit][rt];	
	//	cout<<"here"<<endl;
	}
	int m=9;if(limit) m=s[pos];
	ll res=0;
	for(int i=0;i<=m;i++)
	{
		if(!rt&&abs(pre-i)<2) continue;
		res+=dfs(pos-1,i,(limit&&(i==m)),rt&&i==0);		
	}
//	cout<<res<<endl;
	f[pos][pre][limit][rt]=res;
	return res;
}
ll workt(ll x)
{
	memset(f,-1,sizeof(f));
	int cnt=0;while(x) {s[++cnt]=x%10;x/=10;}
//	cout<<cnt<<endl;
	return dfs(cnt,-2,1,1);
}
int main()
{
	ll A,B;scanf("%lld%lld",&A,&B);
	printf("%lld",workt(B)-workt(A-1));
//	cout<<workt(0)<<endl;
	return 0;
}

第二题（你觉得我会告诉你我只做了两题吗）

同类分布 https://www.luogu.org/problem/P4127

1.比刚才那题难一些（这是废话）

2.惊人发现：判断是否非0竟然比判断是否是-1快！（好蠢啊www）

3.似乎没啥了呢，思维难度高一些，代码没啥区别。。。（x是位置，y是limit，a是数字和，b是mod的余数）

代码（忽略我为了卡常的丑陋www不要欺负我）

#include <iostream>
#include <cstring>
#define ll long long
using namespace std;
ll f[27][3][300][300];
int vis[27][3][300][300];
int s[27];
int mod;
void chai(ll x)
{
	int cnt=0;
	while(x>0)
	{
		cnt++;
		s[cnt]=x%10;
		x/=10;
	}
}
ll find(int x,int y,int a,int b)
{
	if(x==0&&a==mod&&b==0) return 1;
	if(x==0) return 0;
	if(a>mod) return 0;
	if(vis[x][y][a][b]) 
		return f[x][y][a][b];
	vis[x][y][a][b]=1;
	int maxv=9;
	if(y==1) maxv=s[x];
	ll ans=0;
	for(int i=0;i<=maxv;i++)
	{
		ans+=find(x-1,y==1&&i==s[x],a+i,(b*10+i)%mod);
	}
	f[x][y][a][b]=ans;
	return ans;
}
ll work(ll x,int y)
{
	memset(s,0,sizeof(s));
	memset(vis,0,sizeof(vis));
	chai(x);
//	for(int i=1;i<=3;i++) cout<<s[i]<<" ";
	mod=y;
	return find(19,1,0,0);
}
int main()
{
	ll a,b;
	cin>>a>>b;
	ll ans=0;
	if(b==1000000000000000000ll) ans =29410615796612778 ;
	for(int i=1;i<=18*9;i++)
	{
		if (b == 1000000000000000000ll) ans = ans - work(a - 1, i);
		else ans+=work(b,i)-work(a-1,i);	
	}
	cout<<ans<<endl;
	return 0;
}

Scales poj3971 http://poj.org/problem?id=3971

很好的一题。

假设一个物品重量是w，那就是要满足w+x=y，其中x,y是能有n个砝码组成的重量。而且我们每一种砝码只有一个，在二进制中表示就是100000（一些0），而且这些1所在的位置各不相同。所以x,y的二进制表示的相同一位上最多只能有一个1。那我们很快发现还要考虑进位的问题，所以用dp[i][0/1]表示枚举到了w的第i位，这一位是否进位了。dp方程的话就分类讨论一下第i位是啥，是否进位四种情况，然后逐一去看上一位在什么情况下（是否进位）能办到。

#include <iostream>
#include <cstdio>
#include <string>
#include <cstring>

using namespace std;
const int maxn=1000005;
int T,n,l,d,s[maxn],f[maxn][2];
//f[i][1]：第i位有进位 
string st;
void work()
{
	int len=st.size();
	st="&"+st;
	for(int i=len;i>=1;i--) s[i]=(st[len-i+1]-'0');
	for(int i=len+1;i<=n;i++) s[i]=0;
	memset(f,0,sizeof(f));f[0][0]=1;
	for(int i=1;i<=n;i++)
	{
		if(s[i]==1)
		{
			f[i][0]=f[i-1][0];
			f[i][1]=(f[i-1][0]+f[i-1][1])%d; 
		}
		else
		{
			f[i][0]=(f[i-1][1]+f[i-1][0])%d;
			f[i][1]=f[i-1][1];
		}
	}
	printf("%d\n",f[n][0]);
}
int main()
{
	
	scanf("%d",&T);
	for(int ti=1;ti<=T;ti++)
	{
		scanf("%d%d%d",&n,&l,&d);
		cin>>st;
		work();
	}
	return 0;
}
 

poj3252 Round numbers

这道题又是一个比较传统的数位dp惹，除了是二进制hhh 但还是因为一个小错调了一会儿，看注释啦。

#include <iostream>
#include <cstdio>
#include <cstring>
 

using namespace std;
const int maxn=105;
int s[maxn],cnt,f[maxn][maxn][maxn][2][2];
int dp(int pos,int num1,int num0,int limit,int rt)
{
//	cout<<pos<<" "<<num1<<" "<<num0<<endl;
	if(pos==0&&num1<=num0) return 1;
	if(pos==0) return 0; 
	if(f[pos][num1][num0][limit][rt]!=-1) return f[pos][num1][num0][limit][rt];
//	if(pos+num0<num1) return 0;
	int ans=0;
	if(!limit||s[pos]==1)
	{
		ans+=dp(pos-1,num1+1,num0,limit&(s[pos]==1),0);
	}
	if(rt) ans+=dp(pos-1,num1,num0,limit&(s[pos]==0),rt);
	//啊啊啊我怎么这么蠢如果有前导0,limit就不存在了啊 
	else ans+=dp(pos-1,num1,num0+1,limit&(s[pos]==0),0);
	f[pos][num1][num0][limit][rt]=ans;
	return ans;
}
int work(int x)
{
	memset(f,-1,sizeof(f));
	cnt=0;
	while(x)
	{
		s[++cnt]=x%2;
		x/=2;
	} 
	return dp(cnt,0,0,1,1);
}
int main()
{
	int l,r;scanf("%d%d",&l,&r);
	printf("%d",work(r)-work(l-1));
	return 0;
}

CF55B

sb卡常气死老娘了哼（我是不会告诉你我其实到现在还是没有过这道题哈哈哈）

有一些有趣，如果我们把每个数都存下来显然会爆空间。机智的某人（反正不是我）想到了1-9的最小公倍数只有2520，也就是说我们可以把一个很大的数mod2520对结果其实没啥影响。（你想如果一个数能整除某个数，那减掉这个数的一些倍数必定也没有问题）

#include <bits/stdc++.h>

using namespace std;
/*
真tm聪明：1-9的最小公倍数是2520 
*/
typedef long long ll;
int p[49]={0,1,2,3,4,5,6,7,8,9,10,12,14,15,18,20,21,24,28,
30,35,36,40,42,45,56,60,63,70,72,84,90,105,120,126,140,168,
180,210,252,280,315,360,420,504,630,840,1260,2520};
int s[25];ll f[25][2][50][2525];
bool vis[25][2][50][2525];
int ma[2520];

int gcd(int x,int y)
{
	if(x%y==0) return y;
	return gcd(y,x%y);
}
int cal(int x,int y)
{
	if(x<y) swap(x,y);
	if(x%y==0) return x;
	return x*y/gcd(x,y);
}
ll dp(int pos,int limit,int lcm,int sum)
{
//	cout<<pos<<" "<<limit<<" "<<lcm<<" "<<sum<<endl; 
	if(pos==0&&sum==0) return 1;
	if(pos==0&&sum%lcm==0) return 1;
	if(pos==0) return 0;
	if(vis[pos][limit][ma[lcm]][sum]) return f[pos][limit][ma[lcm]][sum];
	vis[pos][limit][ma[lcm]][sum]=1;
	int maxv=9;
	if(limit) maxv=s[pos];
	ll ans=0;
	for(int i=0;i<=maxv;i++)
	{
		if(i==0)
		{
			ans+=dp(pos-1,limit&(i==maxv),lcm,sum*10%2520);
		}
		else ans+=dp(pos-1,limit&(i==maxv),cal(lcm,i),(sum*10+i)%2520);
	} 
	f[pos][limit][ma[lcm]][sum]=ans;
	return ans;
}
ll work(ll x)
{
	memset(vis,0,sizeof(vis));
	memset(f,-1,sizeof(f));
	int cnt=0;
	while(x)
	{
		s[++cnt]=x%10;
		x/=10;	
	} 
	return dp(cnt,1,1,0);
}
int main()
{
	ios::sync_with_stdio(false);
	int cnt=0;
	int T;cin>>T;
	 for(int i=0;i<49;i++)   ma[p[i]]=i;
	for(int ti=1;ti<=T;ti++)
	{
		ll l,r;cin>>l>>r;
		cout<<work(r)-work(l-1)<<endl;	
	}
	return 0;
 }