hdu3652

数位dp

分析题目需求,需要求数字中含有13且能被13整除的数字,是否含有13很容易判断,数位dp时记录一下是否含有13即可,若含有,则不论首位,如不含,则首位为1且次位为3的才满足题意,但是确定能否被13整除就优点麻烦了,按照一般数位dp的方法,由于从高位到低位进行统计,并不知道某个数能否被13整除,无法确定该数字是否满足题意。

再分析,对任意一个数字对13取余均只有13中情况,那么在统计时,记录高位对13取余的结果f,再加上低位对13取余等于13-f的数的个数,不就是需要的对13取余等于0的数字个数吗?照这个思路,增加一维记录对13取余之后为0,1,2.。。。。12的数字个数,便能轻易统计出能被13整除的数的个数。

/*
数位dp,保留余数是为了能与高位数一起确定是否能被13整除,记录是否含13是为了确定在前面不含13的情况下
能否直接累加数据,前两项含义和一般数位dp相同,i位数首位为t的记录。 
*/
#include 
#include 
#include 

using namespace std;

int dp[12][12][2][15];

int main()
{
    memset(dp,0,sizeof(dp));
    int q=10;
    for (int i=0;i<=9;i++)
        dp[1][i][0][i]=1;
    for (int i=2;i<=10;i++)
    {
        for (int t=0;t<=9;t++)
        {
            for (int k=0;k<=9;k++)
            {
                int s=(t*q)%13;
                for (int j=0;j<=12;j++)
                {
                    dp[i][t][1][j]+=dp[i-1][k][1][(13-s+j)%13];                  
                    if (t==1&&k==3)
                        dp[i][t][1][j]+=dp[i-1][k][0][(13-s+j)%13];
                    else
                        dp[i][t][0][j]+=dp[i-1][k][0][(13-s+j)%13];
                }
            }
        }
        q=q*10;
    } 
    int a;
    while (cin>>a)
    {
        int sum=0;
        a=a+1;
        string s;
        int r=1;
        while (a>0)
        {
            int t=a%10;
            char a1=t+48;
            s=a1+s;
            a=a/10;
            if (a>0)
                r=r*10;
        }
        int l=s.length()-1;
        int m=0;
        int p=0;
        for (int i=0;i<=l;i++)
        {
            int c=s[i]-48;
            int f=m%13;
            for (int t=0;t


你可能感兴趣的:(动态规划)