Codeforces 55D 数位DP

本文介绍了一种算法,用于计算指定范围内漂亮数的数量。漂亮数是指可以被其所有非零数字整除的数字。文章提供了详细的算法实现过程,并通过递归函数进行状态转移计算。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

/*
找出l到r内有多少个漂亮数
漂亮数的定义是该数能被每个数位上的非零数字整除
如250是漂亮数 250%2=250%5=0

写出的bug
1.注意返回dp值得调节是flag要为false
2.下一个 状态的计算也错了。不应该写入是拆出的数位,应该是枚举的i
*/
#include <cstring>
#include <cstdio>
#include <algorithm>
#include <cmath>
#include <iostream>
#include <set>
using namespace std;
#define ll long long
int gcd(int a,int b){return b==0?a:gcd(b,a%b);}
int id[6666];
int last=0;
int get(int a,int b){
    return a/gcd(a,b)*b;
}
ll dp[20][2520][66];
int q[666];
int tail=0;
ll dfs(int now,int mod,int lcm,bool flag){
    if(now<0) return mod%lcm ==0;
    if(!flag && dp[now][mod][id[lcm]] != -1) return dp[now][mod][id[lcm]];
    int limit = flag?q[now]:9;
    ll num=0;
    for(int i=0;i<=limit;i++){
        if(i!=0)
            num += dfs(now-1, (mod*10+i)%2520, get(lcm,i), flag&&(i==limit));
        else
            num += dfs(now-1, mod*10%2520, lcm, flag&&(i==limit));
    }
    if(!flag) dp[now][mod][id[lcm]] = num;
    return num;
}

ll W(ll x){
    for(tail=0;x;x/=10){
        q[tail++]=x%10;
    }
    return dfs(tail-1,0,1,true);
}
int main()
{
    //freopen("in.txt","r",stdin);
    memset(dp,-1,sizeof(dp));
    int t;
    scanf("%d",&t);
    set<int > ss;
    for(int st=0;st<(1<<9);st++){
        int LCM=1;
        for(int i=0;i<9;i++){
            if(st&(1<<i)) LCM=get(LCM,i+1);
        }
        ss.insert(LCM);
    }
    for(int x : ss){
        id[x] = last++;
    }
    while(t--){
        ll l,r;
        scanf("%I64d%I64d",&l,&r);
        //cout<<W(r)<<endl;
        //cout<<W(l-1)<<endl;
        printf("%I64d\n",W(r) - W(l-1));
    }
    return 0;
}

区间DP是一种动态规划的方法,用于解决区间范围内的问题。在Codeforces竞赛中,区间DP经常被用于解决一些复杂的字符串或序列相关的问题。 在区间DP中,dp[i][j]表示第一个序列前i个元素和第二个序列前j个元素的最优解。具体的转移方程会根据具体的问题而变化,但是通常会涉及到比较两个序列的元素是否相等,然后根据不同的情况进行状态转移。 对于区间长度为1的情况,可以先进行初始化,然后再通过枚举区间长度和区间左端点,计算出dp[i][j]的值。 以下是一个示例代码,展示了如何使用区间DP来解决一个字符串匹配的问题: #include <cstdio> #include <cstring> #include <string> #include <iostream> #include <algorithm> using namespace std; const int maxn=510; const int inf=0x3f3f3f3f; int n,dp[maxn][maxn]; char s[maxn]; int main() { scanf("%d", &n); scanf("%s", s + 1); for(int i = 1; i <= n; i++) dp[i][i] = 1; for(int i = 1; i <= n; i++) { if(s[i] == s[i - 1]) dp[i][i - 1] = 1; else dp[i][i - 1] = 2; } for(int len = 3; len <= n; len++) { int r; for(int l = 1; l + len - 1 <= n; l++) { r = l + len - 1; dp[l][r] = inf; if(s[l] == s[r]) dp[l][r] = min(dp[l + 1][r], dp[l][r - 1]); else { for(int k = l; k <= r; k++) { dp[l][r] = min(dp[l][r], dp[l][k] + dp[k + 1][r]); } } } } printf("%d\n", dp[n]); return 0; } 希望这个例子能帮助你理解区间DP的基本思想和应用方法。如果你还有其他问题,请随时提问。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值