CodeForces - 55D Beautiful numbers(数位dp)

题目链接:点击这里

题目大意:
[ l , r ] [l,r] [l,r] 有多少种漂亮数字,其定义为能被他所有非零数字整除的数字

题目分析:
数位dp

如果一个数能整除的他的所有非零数位,相当于他能整除所有非零数位的最小公倍数,所以在记忆化是我们可以维护一下数位前缀的 l c m lcm lcm
我们发现 l c m ( 1 , 2 , . . . , 9 ) = 2520 lcm(1,2,...,9)=2520 lcm(1,2,...,9)=2520 ,所以可以利用取余的性质:若 a m o d    b = 0 a\mod b=0 amodb=0 c m o d    b = 0 c\mod b=0 cmodb=0 那么 ( a m o d    c ) m o d    b = 0 (a\mod c) \mod b=0 (amodc)modb=0 简化状态,所以我们可以设计状态: d p [ i ] [ j ] [ k ] dp[i][j][k] dp[i][j][k] 表示处理了前 i i i 位,前面数的 l c m lcm lcm j j j ,和是 k k k ,由于 1 1 1~ 9 9 9 组成的 l c m lcm lcm 只有 48 48 48 个,所以可以先离散化一下,这样数组的第二维只需要开 48 48 48 左右即可;数组第三维利用前面所述的性质开 2520 2520 2520 左右即可

具体细节见代码:

//#pragma GCC optimize(2)
//#pragma GCC optimize("Ofast","inline","-ffast-math")
//#pragma GCC target("avx,sse2,sse3,sse4,mmx")
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#include<set>
#include<map>
#include<stack>
#include<queue>
#define ll long long
#define inf 0x3f3f3f3f
#define int  ll
#define endl '\n'
#define IOS ios::sync_with_stdio(0); cin.tie(0); cout.tie(0)
using namespace std;
int read()
{
	int res = 0,flag = 1;
	char ch = getchar();
	while(ch<'0' || ch>'9')
	{
		if(ch == '-') flag = -1;
		ch = getchar();
	}
	while(ch>='0' && ch<='9')
	{
		res = (res<<3)+(res<<1)+(ch^48);//res*10+ch-'0';
		ch = getchar();
	}
	return res*flag;
}
const int maxn = 1e6+5;
const int mod = 1e9+7;
const double pi = acos(-1);
const double eps = 1e-8;
int n,k,a[maxn],dp[25][50][2555],mp[maxn],cnt;
int _lcm(int a,int b)
{
	return a/__gcd(a,b)*b;
}
const int maxlcm = 2520;
int dfs(int pos,int lcm,int sum,bool f1)
{
	if(pos == -1) return sum%lcm==0;
	if(!f1 && dp[pos][mp[lcm]][sum] != -1) return dp[pos][mp[lcm]][sum];
	int up = f1 ? a[pos] : 9,res = 0;
	for(int i = 0;i <= up;i++)
	{
		int cursum = (sum*10+i)%maxlcm;
		int curlcm = i ? _lcm(lcm,i) : lcm;
		res += dfs(pos-1,curlcm,cursum,f1&&i==up);
	}
	if(!f1) dp[pos][mp[lcm]][sum] = res;
	return res;
}
int calc(int x)
{
	int cnt = 0;
	while(x)
	{
		a[cnt++] = x%10;
		x /= 10;
	}
	return dfs(cnt-1,1,0,1);
}
signed main() {
	for(int i = 1;i <= maxlcm;i++)
		if(maxlcm%i == 0) mp[i] = ++cnt;
	memset(dp,-1,sizeof(dp));
	int t = read();
	while(t--)
	{
		int l = read(),r = read();
		printf("%lld\n",calc(r)-calc(l-1));
	}
	return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值