Codeforces Round #268 (Div. 1) C. Hack it! 数位dp + 构造数位

30 篇文章 0 订阅

传送门

文章目录

题意:

定义 f ( x ) f(x) f(x)表示 x x x的十进制下数位和,现在给你 a a a,让你选一个区间 [ l , r ] [l,r] [l,r],满足 ∑ i = l r f ( i )   m o d   a = 0 \sum_{i=l}^rf(i)\bmod a=0 i=lrf(i)moda=0

1 ≤ a ≤ 1 e 18 1\le a\le1e18 1a1e18

要求输出的 l , r l,r l,r满足 1 ≤ l ≤ r ≤ 1 0 200 1\le l\le r\le 10^{200} 1lr10200

思路:

考虑这样一个事情,对于 x x x,我们找一个 1 0 k > x 10^k>x 10k>x,那么此时 f ( 1 0 k + x ) − f ( x ) = 1 f(10^k+x)-f(x)=1 f(10k+x)f(x)=1

有了这个之后,我们自然的想构造出来 a a a个这样的对,这样模 a a a就是 0 0 0了。

取一个比 a a a大的 10 10 10的幂次 1 0 19 10^{19} 1019,所以就会想到构造这样的区间 [ a + 1 , 1 0 19 + a ] [a+1,10^{19}+a] [a+1,1019+a],转换成前缀和的形式 [ 1 , 1 0 19 + a ] − [ 1 , a ] [1,10^{19}+a]-[1,a] [1,1019+a][1,a],显然这是不行的,虽然这样保证了 [ 1 , a ] , [ 1 0 19 + 1 , 1 0 19 + a ] [1,a],[10^{19}+1,10^{19}+a] [1,a],[1019+1,1019+a]两个区间对应数 f f f做差为 1 1 1,并且有 a a a个这样的数,但是忽略了 [ a + 1 , 1 0 19 ] [a+1,10^{19}] [a+1,1019]这段区间对答案的贡献,那么既然去不掉这段区间,那么不妨将这段区间加入答案中,也就是我们上面说的 a a a对这样的数,假设多出来 k k k,那么我们只需要添加 a − k a-k ak对即可,考虑这些多出来的怎么计算。

我们考虑求出来 ∑ i = 1 1 0 19 f ( i )   m o d   a \sum_{i=1}^{10^{19}}f(i)\bmod a i=11019f(i)moda,设其为 k k k,换句话说,我们现在需要加上 a − k a-k ak 1 1 1就可以达到要求了。考虑将原本区间 [ 1 , 1 0 19 ] [1,10^{19}] [1,1019]不断右移,比如右移一次 [ 2 , 1 0 19 + 1 ] [2,10^{19}+1] [2,1019+1],不难发现,这个时候原本的 1 1 1的贡献变成了 1 + 1 0 19 1+10^{19} 1+1019的贡献,增加了 1 1 1,所以我们右移 a − k a-k ak次,得到区间 [ a − k + 1 , 1 0 19 + a − k ] [a-k+1,10^{19}+a-k] [ak+1,1019+ak],答案就是这个区间啦。

取一个比较大的 10 10 10的幂次,防止移动过程中超过这个幂次即可。

// Problem: C. Hack it!
// Contest: Codeforces - Codeforces Round #268 (Div. 1)
// URL: https://codeforces.com/problemset/problem/468/C
// Memory Limit: 256 MB
// Time Limit: 1000 ms
// 
// Powered by CP Editor (https://cpeditor.org)

//#pragma GCC optimize("Ofast,no-stack-protector,unroll-loops,fast-math")
//#pragma GCC target("sse,sse2,sse3,ssse3,sse4.1,sse4.2,avx,avx2,popcnt,tune=native")
//#pragma GCC optimize(2)
#include<cstdio>
#include<iostream>
#include<string>
#include<cstring>
#include<map>
#include<cmath>
#include<cctype>
#include<vector>
#include<set>
#include<queue>
#include<algorithm>
#include<sstream>
#include<ctime>
#include<cstdlib>
#include<random>
#include<cassert>
#define X first
#define Y second
#define L (u<<1)
#define R (u<<1|1)
#define pb push_back
#define mk make_pair
#define Mid ((tr[u].l+tr[u].r)>>1)
#define Len(u) (tr[u].r-tr[u].l+1)
#define random(a,b) ((a)+rand()%((b)-(a)+1))
#define db puts("---")
using namespace std;

//void rd_cre() { freopen("d://dp//data.txt","w",stdout); srand(time(NULL)); }
//void rd_ac() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//AC.txt","w",stdout); }
//void rd_wa() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//WA.txt","w",stdout); }

typedef long long LL;
typedef unsigned long long ULL;
typedef pair<int,int> PII;

const int N=1000010,INF=0x3f3f3f3f;
const double eps=1e-6;

LL f[30][10*20],mod;
int a[30];

LL dp(int pos,LL sum,int flag) {
	if(pos==0) return sum%mod;
	if(f[pos][sum]!=-1&&flag) return f[pos][sum];
	int x=flag? 9:a[pos];
	LL ans=0;
	for(int i=0;i<=x;i++) (ans+=dp(pos-1,sum+i,flag||i<x))%=mod;
	if(flag) f[pos][sum]=ans;
	return ans; 
}

LL solve() {
	// int tot=0;
	// while(x) a[++tot]=x%10,x/=10;
	// return dp(tot,0,0);
	for(int i=1;i<=19;i++) a[i]=0;
	a[20]=1;
	return dp(20,0,0);
}

inline __int128 read()
{
    __int128 x=0,f=1;
    char ch=getchar();
    while(ch<'0'||ch>'9')
    {
        if(ch=='-')
            f=-1;
        ch=getchar();
    }
    while(ch>='0'&&ch<='9')
    {
        x=x*10+ch-'0';
        ch=getchar();
    }
    return x*f;
}

inline void write(__int128 x)
{
    if(x<0)
    {
        putchar('-');
        x=-x;
    }
    if(x>9)
        write(x/10);
    putchar(x%10+'0');
}


int main()
{
//	ios::sync_with_stdio(false);
//	cin.tie(0);

	memset(f,-1,sizeof(f));
	scanf("%lld",&mod);
	LL k=solve()%mod;
	write((__int128)mod-k+1); putchar(' ');
	write((__int128)100000000000000000*100+mod-k);
	puts("");
	


	return 0;
}
/*

*/









  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值