bzoj 1974: [Sdoi2010]代码拍卖会

Description

随着iPig在P++语言上的造诣日益提升，他形成了自己一套完整的代
码库。猪王国想参加POI的童鞋们都争先恐后问iPig索要代码库。iPi
g不想把代码库给所有想要的小猪，只想给其中的一部分既关系好又
肯出钱的小猪，于是他决定举行了一个超大型拍卖会。 在拍卖会上
，所有的N头小猪将会按照和iPig的好感度从低到高，从左到右地在i
Pig面前站成一排。每个小猪身上都有9猪币（与人民币汇率不明），
从最左边开始，每个小猪依次举起一块牌子，上面写上想付出的买代
码库的猪币数量（1到9之间的一个整数）。大家都知道，如果自己付
的钱比左边的猪少，肯定得不到梦寐以求的代码库，因此从第二只起
，每只猪出的钱都大于等于左边猪出的价钱。最终出的钱最多的小猪
（们）会得到iPig的代码库真传，向着保送PKU（Pig Kingdom Unive
rsity）的梦想前进。 iPig对自己想到的这个点子感到十分满意，在
去现场的路上，iPig就在想象拍卖会上会出现的场景，例如一共会出
现多少种出价情况之类的问题，但这些问题都太简单了，iPig早已不
敢兴趣了，他想要去研究更加困难的问题。iPig发现如果他从台上往
下看，所有小猪举的牌子从左到右将会正好构成一个N位的整数，他
现在想要挑战的问题是所有可能构成的整数中能正好被P整除的有多
少个。由于答案过大，他只想要知道答案mod 999911659就行了。
Input
一行：两个数N（1≤N≤10^18）、P(1≤P≤500)，用一个空格分开。
Output
一行：一个数，表示答案除以999911659的余数。
Sample Input
2 3
Sample Output
15
样例解释
方案可以是：12 15 18 24 27 33 36 39 45 48 57 66 69 78 99，共15种。

解题报告:

这题真的可以玩死人啊,见过的最恶心,没有之一,首先要明白,这题中因为存在数列递增的性质,所以可以拆成1,11,111,1111,11111...中任意选9个以下的数字的和构成的数字,这也是此题给我的最大启发.
然后就是DP,首先看数据范围我们只能从p下手,我们先预处理出1,11,111...中cnt[i]表示mod P后==i的数的个数,这个可以发现是循环的,处理完这里你就胜利了,DP十分简单,定义f[i][j][k] 表示已经处理到了前cnt[i],和mod P为j,一共用了k个数字的方案,容易看出f[i+1][(j+l×i)][k+l]=f[i][j][k]×C[cnt[i]+l-1][l],C[cnt[i]+l-1][l] 是从cnt[i]个中选出l个的可重组合方案

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
long long n;
int k,f[505][10][505],dep[505],inv[10],c[505][10];
long long cnt[505];
long long p=999911659;
int main()
{    scanf("%lld%d",&n,&k);
    int s=0;
    cnt[0]++;
    int beg,len;
    if(n<=k)for(int i=1;i<=n;i++)cnt[s=(s*10+1)%k]++,cnt[s]%=p;
    else
    {    for(int i=1;i<=k+1;i++)
        {    s=(s*10+1)%k;
            if(cnt[s])
            {    beg=dep[s];
                len=i-dep[s];
                break;
            }
            cnt[s]++;
            dep[s]=i;
        }
        for(int i=0;i<k;i++)if(cnt[i]>0&&dep[i]>=beg)
        {    cnt[i]=(n-beg+1)/(long long)len%p;
            if((dep[i]-beg+1)<=(n-beg+1)%(long long)len)cnt[i]=(cnt[i]+1)%p;
            if((dep[i]-beg+1)%(long long)len==(n-beg+1)%(long long)len)s=i;
        }
    }
    inv[1]=1;
    for(int i=2;i<=8;i++)inv[i]=p-(long long)(p/i)*inv[p%i]%p;
    memset(f,0,sizeof 0);
    memset(c,0,sizeof c);
    for(int i=0;i<k;i++)
    {    c[i][0]=1;
        for(int j=1;j<=8;j++)
        {    c[i][j]=cnt[i]*c[i][j-1]%p*inv[j]%p;
            cnt[i]++;
            cnt[i]%=p;
        }
    }
    f[0][0][s]=1;
    for(int i=0;i<k;i++)
    {    for(int j=0;j<9;j++)
        {    for(int x=0;x<k;x++)
            {    for(int y=0;y<=j;y++)
                {    f[i+1][j][x]=(f[i+1][j][x]+(long long)f[i][j-y][(x-y*i%k+k)%k]*c[i][y]%p)%p;
                }
            }
        }
    }
    printf("%d\n",f[k][8][0]);
    return 0;
}