UVa 10069 - Distinct Subsequences

这一题的状态方程不是特别难想,还可以用滚动数组优化,就是高精度略嫌麻烦，我把答案转换成字符串来输出了

#include <iostream>
    #include <cstdio>
    #include <cstdlib>
    #include <cstring>
    #include <ctime>
    using namespace std;
    const int xmax=10000+10;
    const int zmax=100+10;
    const int maxn=30;
    const int base=10000;
    char x[xmax],z[zmax];
    int f[xmax][maxn],last[xmax][maxn];//last记录上一次的d[x][z]
    int n;
    void add_bignum(int a[],int b[],int c[])
    {
        for(int i=maxn-1;i>=1;i--)
            {
                int sum=b[i]+c[i];
                a[i]+=sum%base;//本位和
                a[i-1]+=sum/base;//进位
            }
        a[0]+=b[0]+c[0];
    }
    void print_ans(int d[])
    {
        int i,len=0,rlen=0;
        char s[120],str[5];//把有效位转换成字符串;
        for(i=0;i<maxn;i++)
            if(d[i])//找到第一个不为零的数组下标i，并把它转换成字符串插入s
            {
                sprintf(str,"%d",d[i]);
                rlen=strlen(str);
                for(int j=0;j<rlen;j++)
                    s[len++]=str[j];
                    s[len]='\0';
                break;
            }
        for(int j=i+1;j<maxn;j++)
            {
                sprintf(str,"%d",d[j]);
                rlen=strlen(str);
                for(int k=0;k<4-rlen;k++)//处理前导零
                    s[len++]='0';
                for(int m=0;m<rlen;m++)
                    s[len++]=str[m];
            }
        s[len]='\0';
        printf("%s\n",s);
    }
    int main()
    {
        //freopen("in.txt","r",stdin);
        //freopen("myout.txt","w",stdout);
        scanf("%d",&n);
        while(n--)
        {
            scanf("%s %s",x,z);
            int xlen=strlen(x),zlen=strlen(z);
            int count=0;
            memset(f,0,sizeof(f));
            memset(last,0,sizeof(last));
            for(int j=xlen-1;j>=0;j--)
               {
                   if(z[zlen-1]==x[j])
                {
                    count++;
                }
                    f[j][maxn-1]=count;
               }
    //        for(int j=0;j<xlen;j++)
    //            printf("%d ",f[j]);
    //            printf("\n");
            for(int i=zlen-2;i>=0;i--)
                for(int j=xlen-1;j>=0;j--)
                    {
                        if(z[i]==x[j])//j=xlen-1,i<zlen-1?
                        {
                            //last[j]=f[j];滚动数组，如果没有高精度的话，这样就可以了，实际上f[j]=d[z+1][x]  
                            memcpy(last[j],f[j],maxn*sizeof(int));
                            //f[j]=last[j+1]+f[j+1];//f[j]=d[z+1][x+1]+d[z][x+1]
                            //f[j]=[0];
                            memset(f[j],0,maxn*sizeof(int));
    //                        for(int k=0;k<maxn;k++)
    //                            printf("%d ",f[j][k]);
    //                            printf("\n");
                            add_bignum(f[j],last[j+1],f[j+1]);
                        }
                        else
                        {
                            //last[j]=f[j];//这个last真不知道为什么会有、、、这几天看一下
                            memcpy(last[j],f[j],maxn*sizeof(int));
                            //f[j]=f[j+1];//d[z][x+1]
                            memcpy(f[j],f[j+1],maxn*sizeof(int));
                        }
                    }
    //        for(int i=0;i<maxn;i++)
    //        printf("%d ",f[0][i]);
    //        printf("\n");
              print_ans(f[0]);
        }
        return 0;
    }