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

来源：互联网发布：ios数据恢复工具编辑：程序博客网时间：2024/05/17 03:07

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

0 0