HDU 3652 B-number

HDU 3652 B-number
题 意:让你求1~n含连续的13并且能被13整出的数
输入范围:
1<=n<=1e9
输入样例:

13
100
200
1000

输入样例:

1
1
2
2

思 路:这题在加了一纬储存 mod 这个地方用秦九韶算法。 这个是一个关键点不只可以算多项式乘法，也可以算取模。

#include<bits/stdc++.h>
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
using namespace std;
const int maxn = 25;
typedef long long ll;
ll dp[maxn][maxn][3];
//dp[pos][mod][sta]
//nstate 0代表  的包含13  1代表以1结尾 2代表包含13
int a[maxn];
int r;
ll dfs(int pos,int mod,int sta,bool limit){
    ll tmp = 0;
    if(pos <= -1){
        return sta == 2 && mod == 0;
    }
    if(!limit && dp[pos][mod][sta] !=-1)return dp[pos][mod][sta];
    int up = limit?a[pos]:9;
    for(int i=0;i<=up;i++){
        int nmod = (mod*10+i)%13;
        int state = sta;
        if(sta == 0 && i == 1){
            state = 1;
        }else if(sta == 1 && i!=1 && i!=3){
            state = 0;
        }else if(sta == 1 && i == 3){
            state = 2;
        }
        tmp += dfs(pos-1,nmod,state,limit&&i==a[pos]);
    }
    if(!limit){
        dp[pos][mod][sta] = tmp;
    }
    return tmp;
}
ll solve(int x){
    int pos = 0;
    while(x){
        a[pos++] = x%10;
        x/=10;
    }
    ll tmp = dfs(pos-1,0,0,true);
    return tmp;
}
void init(){
    memset(dp,-1,sizeof(dp));
}
int main(){
    while(~scanf("%d",&r)){
        init();
        printf("%lld\n",solve(r));
    }
    return 0;
}
//