zoukankan      html  css  js  c++  java
  • codeforce 1073E. Segment Sum

    看到这个就是数位DP了,然而细节极多,对于i=1状态直接判了,还有最后一位直接算了

    设f[i][zt][0/1]表示枚举到第i位,用了那些数字,是否有前导0(前导0不计入数字,否则就不知道后面有没有0了)的数的和,g是数的个数

    转移看代码

    #include<cstdio>
    #include<iostream>
    #include<cstring>
    #include<cstdlib>
    #include<algorithm>
    #include<cmath>
    using namespace std;
    typedef long long LL;
    const LL mod=998244353;
    
    int K,n,a[30];
    LL mi[30],f[30][2100][2],g[30][2100][2];int o[2100];
    LL calc(LL G)
    {
        memset(a,0,sizeof(a));n=0;
        LL d=G;
        while(d>0)a[++n]=d%10,d/=10;
        
        LL ret=0;int z=0;d=0;
        for(int i=n;i>=2;i--)
        {
            d=(d*10)%mod;
            for(int k=0;k<a[i];k++)
            {
                LL t=ret;
                if(k==0&&i==n)
                {
                    for(int zt=(1<<10)-1;zt>=0;zt--)
                        if(o[zt]<=K)
                        {
                            ret=(ret+f[i-1][zt][0]+f[i-1][zt][1])%mod;
                        }
                }
                else
                {
                    for(int zt=(1<<10)-1;zt>=0;zt--)
                    {
                        if(o[zt|z|(1<<k)]<=K)
                            ret=(ret+(mi[i-1]*(d+k)%mod*g[i-1][zt][0]%mod)+f[i-1][zt][0])%mod;
                        if(o[zt|z|(1<<k)|1]<=K)
                            ret=(ret+(mi[i-1]*(d+k)%mod*g[i-1][zt][1]%mod)+f[i-1][zt][1])%mod;
                    }
                }
                int p;p++;
            }
            z|=(1<<a[i]);d=(d+a[i])%mod;
        }
        for(int k=0;k<=a[1];k++)
            if(o[z|(1<<k)]<=K)ret+=(d*10+k)%mod;
        return ret;
    }
    
    int main()
    {
        mi[0]=1;for(int i=1;i<=22;i++)mi[i]=mi[i-1]*10%mod;
        memset(f,0,sizeof(f));
        memset(g,0,sizeof(g));
        g[1][0][1]=1;
        for(int k=1;k<=9;k++)f[1][1<<k][0]=k,g[1][1<<k][0]=1;
        for(int i=2;i<=20;i++)
            for(int zt=(1<<10)-1;zt>=0;zt--)
            {
                f[i][zt][1]=(f[i][zt][1]+f[i-1][zt][0]+f[i-1][zt][1])%mod;
                g[i][zt][1]=(g[i][zt][1]+g[i-1][zt][0]+g[i-1][zt][1])%mod;
                //k==0;
                
                for(int k=1;k<=9;k++)
                    if(zt&(1<<k))
                    {
                        f[i][zt][0]=(f[i][zt][0]+(mi[i-1]*k%mod*g[i-1][zt^(1<<k)][0]%mod)+f[i-1][zt^(1<<k)][0])%mod;
                        g[i][zt][0]=(g[i][zt][0]+g[i-1][zt^(1<<k)][0])%mod;
                        
                        f[i][zt][0]=(f[i][zt][0]+(mi[i-1]*k%mod*(g[i-1][zt][0])%mod)+f[i-1][zt][0])%mod;
                        g[i][zt][0]=(g[i][zt][0]+g[i-1][zt][0])%mod;
                        if(zt&1)
                        {
                            f[i][zt][0]=(f[i][zt][0]+(mi[i-1]*k%mod*(g[i-1][zt^(1<<k)][1]+g[i-1][zt^(1<<k)^1][1])%mod)+f[i-1][zt^(1<<k)][1]+f[i-1][zt^(1<<k)^1][1])%mod;
                            g[i][zt][0]=(g[i][zt][0]+g[i-1][zt^(1<<k)][1]+g[i-1][zt^(1<<k)^1][1])%mod;
                            
                            f[i][zt][0]=(f[i][zt][0]+(mi[i-1]*k%mod*(g[i-1][zt][1]+g[i-1][zt^1][1])%mod)+f[i-1][zt][1]+f[i-1][zt^1][1])%mod;
                            g[i][zt][0]=(g[i][zt][0]+g[i-1][zt][1]+g[i-1][zt^1][1])%mod;
                        }
                    }
            }
        memset(o,0,sizeof(o));
        for(int zt=(1<<10)-1;zt>=0;zt--)
            for(int k=0;k<=9;k++)
                if(zt&(1<<k))o[zt]++;
                
        LL L,R;
        scanf("%lld%lld%d",&L,&R,&K);
        printf("%lld
    ",((calc(R)-calc(L-1))%mod+mod)%mod);
        return 0;
    }
  • 相关阅读:
    30天敏捷结果(17):找出高效时间,并利用它来处理重要事情
    需求:需求获取技术之 观察
    30天敏捷结果(13):如何对事情付诸行动
    2010年11月blog汇总:敏捷练习、需求和建模
    MDSF:Eclipse MDD Day学习
    需求:结合TOGAF做好需求获取工作
    MDSF:架构工具简要功能说明
    信息系统开发平台OpenExpressApp:【OpenTest】 之 框架实现说明
    吃:第二次去吃香草香草
    30天敏捷生活(911):调整30天的节奏
  • 原文地址:https://www.cnblogs.com/AKCqhzdy/p/9860813.html
Copyright © 2011-2022 走看看