bzoj1974&jzoj1523 [Sdoi2010]auction 代码拍卖会

168 篇文章 0 订阅

题意:当一个数的每一位都比前一位大时,我们把这种数叫做上升数。现在给你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;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值