LightOJ 1068 Investigation (数位DP)

题意求出区间[A,B]内能被K整除且各位数字之和也能被K整除的数的个数。

分析:A,B最多有10位    各位数字之和不会超过90,那么当 k >= 90时,为0.

dp[len][y1][y2]    第len位   y1这个数除k的余数  y2位各个位和除k的余数.



#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<string>
#include<vector>
#include<cctype>
#include<set>
#include<map>
#include<queue>
#include<stack>
#include<iomanip>
#include<sstream>
#include<limits>
#define ll long long
#define inf 0x3f3f3f3f
using namespace std;
const int N= 2000;
const int maxn = 1e3+10;
const int maxm = 2000000001;
const int MOD = 1000;
int dp[20][100][100],k,dig[maxn];
int dfs(int len ,int y1, int y2, int ok)
{
    if (len == 0) return y1 == 0 && y2 == 0;
    if(!ok && dp[len][y1][y2] != -1) return dp[len][y1][y2];
    int ans = 0, res = ok? dig[len]:9;    //如果上一位到最大的数则这一位不能取9.
    for(int i = 0; i <= res; i++)
      ans += dfs(len-1,(y1*10+i) % k,(y2+i)%k,ok && i==res);
    if (!ok) dp[len][y1][y2] = ans;    //这一位到了最大数,则不能赋值,因为不能为后面所用。
    return ans;
}
int work(int n)
{
    int cnt = 1;
    while(n)
    {
        dig[cnt++] = n %10;
        n /= 10;
    }
    return dfs(cnt-1,0,0,1);
}
int main(){
#ifdef LOCAL
	freopen("C:\\Users\\lanjiaming\\Desktop\\acm\\in.txt","r",stdin);
	//freopen("output.txt","w",stdout);
#endif
ios_base::sync_with_stdio(0);
    int T,kase = 0;
    cin>>T;
    while(T--)
    {
        cout<<"Case "<<++kase<<": ";
        int a,b;
        cin>>a>>b>>k;
        if (k > 90)
        {
            cout<<0<<endl;
            continue;
        }
        memset(dp,-1,sizeof(dp));
        cout<<work(b) - work(a-1)<<endl;
    }
	return 0;
}





  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
Sigma函数是指一个数字的所有因子之和。给定一个数字n,需要求出有多少个数字的Sigma函数是偶数。\[2\] 为了解决这个问题,可以先筛选出n范围内的素数(范围在10^6即可),然后对n进行素因子分解。对于每个因子,如果它的Sigma函数中连乘的每一项都是偶数,那么整个Sigma函数就是偶数。具体实现中,可以判断每个因子的平方根是否为偶数,如果是偶数,则减去(平方根+1)/2。\[1\] 另外,还可以使用O(1)的做法来解决这个问题。根据观察,所有的完全平方数及其两倍的值都会导致Sigma函数为偶数。因此,可以直接计算n的平方根,然后减去(平方根+1)/2即可得到结果。\[3\] #### 引用[.reference_title] - *1* [Sigma Function](https://blog.csdn.net/PNAN222/article/details/50938232)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v91^control_2,239^v3^insert_chatgpt"}} ] [.reference_item] - *2* *3* [【LightOJ1336】Sigma Function(数论)](https://blog.csdn.net/qq_30974369/article/details/79009498)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v91^control_2,239^v3^insert_chatgpt"}} ] [.reference_item] [ .reference_list ]
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值