給定X和Y的值，一個數字用十進位制表示，對於數位4出現恰好X次，7恰好Y次的數字稱為nya數。

然後給P和Q， 後面是一系列查詢， 每個查詢K就是在區間(P，Q]上找到第K個nya數，不存在則輸出Nya!（= = 這是尖叫的意思麼）

因為數字的範圍比較大，又是跟數位有關，所以採取數位DP的演算法。

dp[i][j][k]，代表長度為i的數字中有j個4和k個7的個數。

先用遞推的方法預處理出dp陣列。

dp[i][j][k] += dp[i-1][j][k]*8，第i位不是4和7的情況；

dp[i][j+1][k] += dp[i-1][j][k]，第i位為4；

dp[i][j][k+1] += dp[i-1][j][k]，第i位為7；

然後讀入p和q的時候先計算出對應的區間裡nya數的個數tmp，然後讀入k，如果k超過tmp就輸出Nya!，否則二分出答案。

#include<cstdio>
#include<cstring>
#define LL long long
#define ull unsigned long long
ull dp[20][21][21];
void init(){
    memset(dp,0,sizeof(dp));
    dp[0][0][0]=1;
    for(int i=1; i<20; i++){
        for(int j=0; j<=20; j++){
            for(int k=0; k<=20; k++)    dp[i][j][k]+=dp[i-1][j][k]*8;
        }
        for(int j=0; j<20; j++){
            for(int k=0; k<=20; k++){
                dp[i][j+1][k]+=dp[i-1][j][k];
                dp[i][k][j+1]+=dp[i-1][k][j];
            }
        }
    }
}
int x, y;
ull query(ull a){
    int buf[20];
    int len=0;
    while(a>0){
        buf[len++]=(int)(a%10);
        a/=10;
    }
    LL sum=0;
    int j=x, k=y;
    for(int i=len-1; i>=0; i--){
        if(buf[i]<=4){
            sum += (LL)buf[i]*dp[i][j][k];
            if(buf[i]==4)   j--;
        }
        else if(buf[i]<=7){
            sum += (LL)(buf[i]-1)*dp[i][j][k];
            if(j)   sum += dp[i][j-1][k];
            if(buf[i]==7)   k--;
        }
        else{
            sum += (LL)(buf[i]-2)*dp[i][j][k];
            if(j)   sum += dp[i][j-1][k];
            if(k)   sum += dp[i][j][k-1];
        }
        if(j<0 || k<0)  break;
    }
    return sum;
}
int t, n;
ull p, q, k, tmp, cur;
int main(){
    init();
    scanf("%d", &t);
    for(int ct=1; ct<=t; ct++){
        printf("Case #%d:\n", ct);
        scanf("%I64u %I64u %d %d", &p, &q, &x, &y);
        tmp = query(q+1)-query(p+1);
        scanf("%d", &n);
        ull low, top, mid;
        while(n--){
            scanf("%I64u", &k);
            if(k>tmp){
                puts("Nya!");
                continue;
            }
            low=p+1; top=q;
            while(low<top){
                mid = (low+top)/2;
                cur = query(mid+1)-query(p+1);
                if(cur<k)   low=mid+1;
                else    top=mid;
            }
            printf("%I64u\n", low);
        }
    }
    return 0;
}