题意:当一个数的每一位都比前一位大时,我们把这种数叫做上升数。现在给你n和k,让你求n以内有多少个上升数能整除,n<=10^18,k<=500。
分析:大概能想到背包,但是细节很复杂所以就炸了。。。不得已膜了一波题解。。
每一个上升数都可以分解为1,11,111,1111之类的数。那么我们可以找出所有的这种树,把他们分成k组,做分组背包。
令f[i][j][k]表示前i组里选j个数的和在模K意义下是k的方案数,由于同一组里选m个数组成可重集的方案数是C(X−1+m,m),其中X表示组里不同数字的个数,所以转移是枚举当前组里选多少个数进行转移。
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
using namespace std;
typedef long long ll;
ll n;
int m;
const int maxn=1e5,mo=999911659;
int a[maxn],pos[maxn],s[maxn],val,inv[maxn],f[11][maxn],ans;
int main()
{
scanf("%lld%d",&n,&m);
a[1]=1%m;
pos[a[1]]=1;
fo(i,2,m+1)
{
a[i]=(a[i-1]*10+1)%m;
if (pos[a[i]])
{
int l=pos[a[i]],r=i,len=r-l;
for(int j=1;j<l&&j<=n;++j)
++s[a[j]];
for(int j=0;j<len&&l+j<=n;++j)
s[a[l+j]]=(s[a[l+j]]+(n-l-j)/len+1)%mo;
if (n<l)val=a[n];
else val=a[l+(n-l)%len];
--s[val];
if (s[val]<0)s[val]+=mo;
break;
}
pos[a[i]]=i;
}
inv[1]=1;
fo(i,2,9)
inv[i]=mo-mo/i*(ll)inv[mo%i]%mo;
f[0][0]=1;
fo(i,0,m-1)
{
if (!s[i])continue;
fd(j,9,1)
{
int tmp=1;
for(int k=1;k<=j&&tmp;++k)
{
tmp=(ll)tmp*(s[i]-1+k)%mo*inv[k]%mo;
fo(l,0,m-1)
f[j][(l+k*i)%m]=(f[j][(l+k*i)%m]+(ll)tmp*f[j-k][l])%mo;
}
}
}
fo(i,1,9)
{
int res=val*i%m;
if (res)res=m-res;
for(int j=0;i+j<10;++j)
{
ans+=f[j][res];
if (ans>=mo)ans-=mo;
}
}
printf("%d\n",ans);
return 0;
}