GZM1 模拟赛 Day25

W

一、题目及数据范围

点此看题

二、解法

一看就是$dp$，考试时候想A，结果写炸了。
我们定义$dp[i]$为处理到第i位的最小花费（处理指把原序列变成全0序列，我们从后往前处理），则有：
$dp[i]=min(dp[j]+sum1(j,i],dp[j]+sum0(j,i]+2)$
我们维护后缀和$cnt$，由于我们只需要找$cnt0[j]+dp[j]$和$cnt1[j]+dp[j]$，的最小值，在$dp$中$O(1)$维护即可。

#include <cstdio>
#include <cstring>
const int MAXN = 1000005;
int n,dp[MAXN],cnt0[MAXN],cnt1[MAXN];
char str[MAXN];
int min(int x,int y)
{
    return x<y?x:y;
}
int main()
{
    scanf("%s",str+1);
    n=strlen(str+1);
    for(int i=n; i>=1; i--)
    {
        cnt0[i]=cnt0[i+1]+(str[i]=='0');
        cnt1[i]=cnt1[i+1]+(str[i]=='1');
    }
    int min1=0,min2=0;
    for(int i=n; i>=1; i--)
    {
        dp[i]=min(min1+cnt1[i],min2+cnt0[i]+2);
        min1=min(min1,dp[i]-cnt1[i]);
        min2=min(min2,dp[i]-cnt0[i]);
    }
    printf("%d\n",dp[1]);
}

inint

一、题目及数据范围

题目描述 从起点1开始，每次选择当前数的任意一位上加上去，问得到n的最小步数以及方案数。多组数据。

二、解法

对于20%的数据，直接暴力$dp$跑。
后面就不行了，我们考虑正解。
我们需要一种高效的在线算法，可以考虑预处理。
由于我们的加法都是针对个位的，所以相对其他位，个位更难控制。我们可以在此基础上，我们可以考虑数位$dp$，把个位的变化特殊考虑。
定义$f[S][d][i][j]$，把数$x$分成两部分，第一部分是$\frac{x}{10^{d+1}}$，另一部分是$x\mathrm{m}\mathrm{o}\mathrm{d}10^{d+1}$，第一部分数的组成状压成$S$，我们要把个位$i$变成$j$，并把$10^{d+1}$上进一位。
定义$fs[S][d][i][j]$为在最小步数小的方案数。
仿照倍增，我们先预处理处个位的情况，也就是$f[S][0][i][j]$，再从$f[d]$转移到$f[d+1]$，每次在$d$位上跳满10下，就可以进位了。
处理答案时，我们从高到低，依次凑给出的数位，最后对个位进行修正即可。

#include <cstdio>
#include <cstring>
#define LL long long
const int MAXN = 512;
const int MOD = 1e9+7;
const LL INF = 0x3f3f3f3f3f3f3f3fll;
LL read()
{
    LL x=0,flag=1;
    char c;
    while((c=getchar())<'0' || c>'9') if(c=='-') flag=-1;
    while(c>='0' && c<='9') x=(x<<3)+(x<<1)+(c^48),c=getchar();
    return x*flag;
}
LL f[MAXN+10][15][12][12],fs[MAXN][15][12][12];
LL tmp[15][15],tmps[15][15],T,n;
bool vis[20];
void get(int S)
{
    for(int i=1; i<=9; i++)
        if(S&(1<<i-1))
            vis[i]=1;
}
void prepare()
{
    memset(f,0x3f,sizeof f);
    for(int S=0; S<MAXN; S++)
        for(int i=9; i>=0; i--)
        {
            memset(vis,0,sizeof vis);
            get(S);
            vis[i]=1;
            for(int j=0; j<=8; j++)
                if(vis[10-i+j]) f[S][0][i][j]=1,fs[S][0][i][j]=1;
                else
                {
                    for(int k=i+1; k<=9; k++)
                        if(vis[k-i])
                        {
                            if(f[S][0][k][j]+1<f[S][0][i][j])
                                f[S][0][i][j]=f[S][0][k][j]+1,fs[S][0][i][j]=0;
                            if(f[S][0][k][j]+1==f[S][0][i][j])
                                fs[S][0][i][j]+=fs[S][0][k][j];
                        }
                }
        }
    for(int i=1; i<=11; i++)
        for(int S=0; S<MAXN; S++)
            for(int s=0; s<=8; s++)
            {
                memset(tmp,0x3f,sizeof tmp);
                tmp[0][s]=0;
                tmps[0][s]=1;
                for(int j=0; j<=9; j++)
                    for(int a=0; a<=8; a++)
                        for(int b=0; b<=8; b++)
                        {
                            int nS=j?S|(1<<j-1):S;
                            if(tmp[j+1][b]>tmp[j][a]+f[nS][i-1][a][b])
                                tmp[j+1][b]=tmp[j][a]+f[nS][i-1][a][b],tmps[j+1][b]=0;
                            if(tmp[j+1][b]==tmp[j][a]+f[nS][i-1][a][b])
                                tmps[j+1][b]=(tmps[j+1][b]+tmps[j][a]*fs[nS][i-1][a][b])%MOD;
                        }
                for(int t=0; t<=8; t++)
                    f[S][i][s][t]=tmp[10][t],fs[S][i][s][t]=tmps[10][t];
            }
}
int main()
{
    prepare();
    T=read();
    while(T--)
    {
        n=read();
        LL tmp[2][15],tmps[2][15];
        memset(tmp,0x3f,sizeof tmp);
        int t=0,S=0;
        tmp[0][1]=0,tmps[0][1]=1;
        for(LL i=1000000000000ll,ti=11; i>=10; i/=10,ti--)
        {
            for(int j=0; j<(n/i)%10; j++)
            {
                t=!t;
                memset(tmp[t],0x3f,sizeof tmp[t]);
                int nS=j?S|(1<<j-1):S;
                for(int a=0; a<=8; a++)
                    for(int b=0; b<=8; b++)
                    {
                        if(tmp[t][b]>tmp[!t][a]+f[nS][ti][a][b])
                            tmp[t][b]=tmp[!t][a]+f[nS][ti][a][b],tmps[t][b]=0;
                        if(tmp[t][b]==tmp[!t][a]+f[nS][ti][a][b])
                            tmps[t][b]=(tmps[t][b]+tmps[!t][a]*fs[nS][ti][a][b])%MOD;
                    }
            }
            if((n/i)%10) S|=1<<((n/i)%10-1);
        }
        for(int i=0; i<=9; i++)
        {
            memset(vis,0,sizeof vis);
            get(S);
            vis[i]=1;
            for(int a=1; a<=9; a++)
                if(vis[a] && i+a<=9)
                {
                    int j=i+a;
                    if(tmp[t][j]>tmp[t][i]+1)
                        tmp[t][j]=tmp[t][i]+1,tmps[t][j]=0;
                    if(tmp[t][j]==tmp[t][i]+1)
                        tmps[t][j]=(tmps[t][j]+tmps[t][i])%MOD;
                }
        }
        if(tmp[t][n%10]==INF) printf("IMPOSSIBLE\n");
        else printf("%lld %lld\n",tmp[t][n%10],tmps[t][n%10]);
    }
}