ZOJ 3962

就是统计1~n中出现的各个数字的次数，当然是在16进制下。

不过有个区间问题的小技巧，统计从 [x,y] 可以转换成 从 [1,y] 减去 [1,x-1]。

不过要分类讨论一下，因为有可能会出现溢出，从ffffffff +1 得到 00000000 就是溢出了。

因为 n < 1e9 所以只会溢出一次。

// 290ms

#include<cstdio>
#include<cmath>
#include<cctype>

const int cost[] = {,,,,,,,,
                    ,,,,,,,};

long long solve(long long n, int k) {
    long long i = ;   // 位数
    long long cnt = ;
    long long low = ;
    long long tmp = n;
    while(n) {
        long long t = n % ;
        n /= ;
        low = tmp % i;
        if(t < k) {
            cnt += n*i;
        } else if(t > k) {
            cnt += (+n)*i;
        } else {
            cnt += n*i;
            cnt += low+;
        }
        i *= ;
    }
    return cnt;
}
int Hex(char c){
    if(isdigit(c))
        return c-'';
    return c - 'A' + ;
}
int main() {
    int t;
    scanf("%d",&t);
    for(int i=;i<t;i++){
        long long n;
        long long x=;
        char str[];
        scanf("%lld%s",&n,str);
        for(int i=;i<;i++)
            x = (x<<) + Hex(str[i]);
        long long cnt = ;
        long long ans = ;

        if((n-)+x <= 4294967295ll){ // 没有溢出
            long long tmp;
            for(int i=;i<;i++){
                tmp = solve((n-1ll)+x,i);
                cnt += tmp;
                ans += tmp*cost[i];
            }
            for(int i=;i<;i++){
                tmp = solve(x-,i);
                cnt -= tmp;
                ans -= tmp*cost[i];
            }
            ans += (n*-cnt)*cost[];
        }else {
            long long tmp;
            for(int i=;i<;i++){
                tmp = solve(4294967295ll,i);
                cnt += tmp;
                ans += tmp*cost[i];
            }
            for(int i=;i<;i++){
                tmp = solve(x-,i);
                cnt -= tmp;
                ans -= tmp*cost[i];
            }
            for(int i=;i<;i++){
                tmp = solve((n-2ll)+x-4294967295ll,i);
                cnt += tmp;
                ans += tmp*cost[i];
            }
            ans += (n*-cnt)*cost[];
        }
        printf("%lld\n",ans);
    }
    return ;
}

搜了一下题解发现可以用数位dp做。。

当时倒是也想到了数位dp但是因为做过的都是不定长的，所以就没有想到可以定长为8。

// 250ms

#include<cstdio>
#include<cmath>
#include<cstring>
#include<cctype>

int digit[];
long long dp[][];
const int cost[] = {,,,,,,,,
                    ,,,,,,,};

long long dfs(int d,int sum,bool shangxian){
    if(d == )return sum;
    if(!shangxian && dp[d][sum] != -)
        return dp[d][sum];
    int maxn = shangxian?digit[d]:;
    long long tmp = ;
    for(int i=;i<=maxn;i++){
        tmp += dfs(d-,sum+cost[i],shangxian && i==maxn);
    }
    return shangxian?tmp:dp[d][sum]=tmp;
}
long long solve(long long x){
    int k = ;
    if(x<)return ;
    memset(digit,,sizeof(digit));
    while(x){
        digit[++k] = x % ;
        x /= ;
    }
    return dfs(,,);
}
int Hex(char c){
    if(isdigit(c))
        return c-'';
    return c - 'A' + ;
}
int main() {
    int t;
    memset(dp,-,sizeof(dp));
    scanf("%d",&t);
    while(t--){
        long long n;
        long long x=;
        char str[];
        scanf("%lld%s",&n,str);
        for(int i=;i<;i++)
            x = (x<<) + Hex(str[i]);
        // 4294967295ll
        if(x+n- <= 4294967295ll)
            printf("%lld\n",solve(x+n-)-solve(x-));
        else
            printf("%lld\n",
                   solve(4294967295ll)-solve(x-)+solve(x+n--4294967295ll));
    }
    return ;
}