[DP] BZOJ 1974 [Sdoi2010]auction 代码拍卖会 & 51Nod 1261 上升数

tls的题解
关键在于 对于一个上升数 比如 1123445
我们习惯性的把他纵向切割 1 1 2 3 4 4 5
但是只要我们转变一下 横向切割 变成 1111111+11111+1111+111+1
就会发现一些很好的性质 首先 最多被切割成9段 然后这些形似$1111\cdots$的数取模肯定是有循环节的 把余数相同的合并成一组 做一个分组背包
记得因为$10^{k−1} \over 9$是必选一个的 最后特殊处理下
sdoi的坑点 K可以等于1

#include<cstdio>
#include<cstdlib>
#include<algorithm>
using namespace std;
typedef long long ll;

const int P=999911659;
const int N=505;

ll n; int K;
ll cnt[N]; int vst[N];
int lst[N],pnt;
int f[N][15][N];
int C[N][15],inv[15];

int main(){
  freopen("t.in","r",stdin);
  freopen("t.out","w",stdout);
  scanf("%lld%d",&n,&K);
  ll t=1%K,r,x,y,xn,i;
  for (i=1;!vst[t] && i<=n;i++,t=(t*10+1)%K) vst[t]=i,lst[++pnt]=t;
  if (i==n+1){
    for (int i=1;i<=pnt;i++) cnt[lst[i]]=1;
    xn=lst[pnt];
  }else{
    t=vst[t]-1,r=pnt-t;
    for (int i=1;i<=t;i++) cnt[lst[i]]++;
    x=(n-t)/r,y=(n-t)%r;
    for (int i=1;i<=r;i++){
      cnt[lst[t+i]]+=x;
      if (i<=y) cnt[lst[t+i]]++;
    }
    y=y==0?r:y;
    xn=lst[t+y];
  }
  inv[1]=1; for (int i=2;i<=9;i++) inv[i]=(ll)(P-P/i)*inv[P%i]%P;
  for (int i=0;i<K;i++){
    C[i][0]=1;
    if (cnt[i])
      for (int t=1;t<=9;t++)
    C[i][t]=(ll)C[i][t-1]*((cnt[i]+t-1)%P)%P*inv[t]%P;
  }
  for (int j=0;j<=9;j++)
    f[0][j][0]=C[0][j];
  for (int i=0;i<K-1;i++)
    for (int j=0;j<=9;j++)
      for (int k=0;k<K;k++)
    if (f[i][j][k])
      for (int t=0;t<=9 && j+t<=9;t++)
        (f[i+1][j+t][(k+(i+1)*t)%K]+=(ll)f[i][j][k]*C[i+1][t]%P)%=P;
  int Ans=0;
  for (int j=0;j<=8;j++)
    (Ans+=f[K-1][j][(K-xn)%K])%=P;
  printf("%d\n",Ans);
  return 0;
}