Codeforces 908G New Year and Original Order 数位dp

最新推荐文章于 2021-01-29 10:53:55 发布

_beginend

最新推荐文章于 2021-01-29 10:53:55 发布

阅读量156

点赞数 1

分类专栏：数位dp

本文链接：https://blog.csdn.net/qq_33229466/article/details/95318198

版权

数位dp 专栏收录该内容

11 篇文章 0 订阅

订阅专栏

该博客介绍了Codeforces 908G问题的解决方案，涉及数位动态规划（dp）的应用。博主详细解析了题意，提出通过dp计算满足特定条件的数的数量，其中dp状态表示为第i位至少为j的数的个数，并讨论了如何进行状态转移。

摘要由CSDN通过智能技术生成

题意

定义 $S (n)$ 为 $n$ 在十进制表示下每一位数字从小到大排序后拼接得到的数字。
求 $\sum_{i=1}^nS(i)$
$n\le10^{700}$

分析

用dp求出 $f_{i,j}$ 表示有多少个 $x$ 满足 $S (x)$ 的第 $i$ 位至少为 $j$ 。
dp的时候设 $f_{i,j,0/1}$ 表示考虑了前 $i$ 位数字，有 $j$ 位不大于当前枚举的数，前 $i$ 位是否严格小于 $n$ 。转移显然。

代码

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>

typedef long long LL;

const int N=705;
const int MOD=1000000007;

int n,f[N][N][2],g[15][N],pre[N],a[N];
char str[N];

void solve(int now)
{
	memset(f,0,sizeof(f));
	f[0][0][1]=1;
	for (int i=0;i<n;i++)
		for (int j=0;j<=i;j++)
		{
			if (f[i][j][0])
			{
				for (int k=0;k<=now;k++) (f[i+1][j+1][0]+=f[i][j][0])%=MOD;
				for (int k=now+1;k<=9;k++) (f[i+1][j][0]+=f[i][j][0])%=MOD;
			}
			if (f[i][j][1])
			{
				for (int k=0;k<a[i+1];k++)
					if (k<=now) (f[i+1][j+1][0]+=f[i][j][1])%=MOD;
					else (f[i+1][j][0]+=f[i][j][1])%=MOD;
				if (a[i+1]<=now) (f[i+1][j+1][1]+=f[i][j][1])%=MOD;
				else (f[i+1][j][1]+=f[i][j][1])%=MOD;
			}
		}
	for (int i=0;i<n;i++) g[now+1][i+1]=(f[n][i][0]+f[n][i][1])%MOD,(g[now+1][i+1]+=g[now+1][i])%=MOD;
}

void calc()
{
	int ans=0;
	for (int i=1;i<=9;i++)
		for (int j=1;j<=n;j++)
			(ans+=(LL)i*pre[j]%MOD*(g[i][j]-g[i+1][j])%MOD)%=MOD;
	ans+=ans<0?MOD:0;
	printf("%d\n",ans);
}

int main()
{
	scanf("%s",str+1);
	n=strlen(str+1);
	for (int i=1;i<=n;i++) a[i]=str[i]-'0';
	pre[1]=1;
	for (int i=2;i<=n;i++) pre[i]=(LL)pre[i-1]*10%MOD;
	std::reverse(pre+1,pre+n+1);
	for (int i=0;i<9;i++) solve(i);
	calc();
	return 0;
}