[hdu3943]K-th Nya Number

　　挺正常的一道模板题。

　　f[i][j][k]表示i位的数，有j个4，k个7的方案数。

　　具体实现的话...我写了发二分答案。。需要注意的是二分时应该是mid=L+(R-L)/2。。不然分分钟爆longlong（unsigned long long党自行退散

　　其实也可以从左端点开始慢慢爬。。。但总觉得比较蛋疼所以没敢写

　　由网上题解可得，其实还可以确定答案的位数后，从高位往低位一个一个试= =...复杂度会比二分答案的少个log

　　需要注意一下对0的特判。。

　　对于从左端点爬到根再爬到右端点的奇怪姿势一直不敢碰。。药丸的节奏啊= =。。GDOI应该不会出这么sxbk的东西吧（手动立flag？

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
struct poi{
    ll k;int id;
}q[233];ll ans[233];
ll f[20][21][21],ten[20];bool u[20][21][21];
int i,j,x,y,X,Y,num;
ll l,r,tmp;
int s[23],len;

bool cmp(poi a,poi b){return a.k<b.k;}
inline ll get(ll x){
    for(tmp=x,len=0;tmp;tmp/=10)s[++len]=tmp%10;
    if(!x)s[len=1]=0;
    ll ans=0;register int i,j;int sx=X,sy=Y;
    if(X+Y>len)return 0;
    for(i=sx+sy;i<len;i++)
        for(j=1;j<=9;j++)
            if((sx||j!=4)&&(sy||j!=7))ans+=f[i-1][sx-(j==4)][sy-(j==7)];
//    printf("     %lld
",ans);
    for(i=1;i<s[len];i++)
        if((sx||i!=4)&&(sy||i!=7))ans+=f[len-1][sx-(i==4)][sy-(i==7)];
    sx-=s[len]==4,sy-=s[len]==7;
    for(i=len-1;i&&sx>=0&&sy>=0;i--){
        for(j=0;j<s[i];j++)
            if((sx||j!=4)&&(sy||j!=7))ans+=f[i-1][sx-(j==4)][sy-(j==7)];
        sx-=s[i]==4,sy-=s[i]==7;
    }
    if(!sx&&!sy&&x)ans++;
//    printf("     %lld
",ans);
    return ans+(!X&&!Y);
}
int main(){
    f[1][1][0]=f[1][0][1]=1;f[1][0][0]=8;f[0][0][0]=1;
    for(i=ten[0]=1;i<=18;i++)ten[i]=ten[i-1]*10;
    for(i=2;i<=19;i++){
        for(x=1;x<=i;x++)for(y=i-x;y;y--)
            f[i][x][y]=x<=y?(f[i-1][x][y]*8+f[i-1][x-1][y]+f[i-1][x][y-1]):f[i][y][x]
        ;//    ,printf("  %d %d %d   %lld
",i,x,y,f[i][x][y]);
        for(x=1;x<=i;x++)f[i][x][0]=f[i][0][x]=f[i-1][x-1][0]+f[i-1][x][0]*8;
        f[i][0][0]=f[i-1][0][0]<<3;
    //    if(i<=3)
    //    for(x=0;x<=i;x++)for(y=0;y<=i-x;y++)printf("  %d %d %d  %lld
",i,x,y,f[i][x][y]);
    }
    int T,TT;scanf("%d",&TT);
    for(T=1;T<=TT;T++){
        scanf("%lld%lld%d%d",&l,&r,&X,&Y);scanf("%d",&num);
        for(i=1;i<=num;i++)scanf("%lld",&q[i].k),q[i].id=i;
        
        sort(q+1,q+1+num,cmp);
        ans[0]=l;
        ll lnum=get(l),mx=get(r)-lnum;
        
    //    printf("    lnum:%lld   mx:%lld
",lnum,mx);
        printf("Case #%d:
",T);
        for(i=1;i<=num;i++){
            ll L=ans[q[i-1].id],R=r,mid;
            if(q[i].k==q[i-1].k){ans[q[i].id]=L;continue;}
            if(q[i].k>mx){ans[q[i].id]=-1;continue;}
            while(L<R){
                mid=L+((R-L)>>1);//printf("   %lld %lld  
",L,R);
                if(get(mid)>=q[i].k+lnum)R=mid;else L=mid+1;
            }//puts("");
            ans[q[i].id]=L;
        }
        for(i=1;i<=num;i++)
            if(ans[i]==-1)puts("Nya!");else printf("%lld
",ans[i]);
    }
    return 0;
}

　　话说网上有些题解遇到极限数据会炸>_<。。