Cod sursa(job #413465)

Utilizator freak93Adrian Budau freak93 Data 8 martie 2010 17:13:33
Problema Cifre Scor 100
Compilator cpp Status done
Runda Arhiva de probleme Marime 2.64 kb
#include<fstream>

using namespace std;

const char iname[]="cifre.in";
const char oname[]="cifre.out";
const int maxn=15;

ifstream f(iname);
ofstream g(oname);

long long a,b,c,k,i,j,cnt[maxn][maxn],cnt0[maxn][maxn],p,pows[maxn],rez,pows2[maxn];

void precalc()
{
    cnt[1][1]=1;
    cnt[1][0]=9;
    for(i=2;i<=12;++i)
        for(j=0;j<=i;++j)
            cnt[i][j]=cnt[i-1][j]*9+(j>0?cnt[i-1][j-1]:0);
    cnt0[1][1]=1;
    cnt0[1][0]=9;
    for(i=2;i<=12;++i)
        for(j=0;j<=i;++j)
            if(c==0)
                cnt0[i][j]=9*cnt[i-1][j]+(j>=0?cnt0[i-1][j]:0);
            else
                cnt0[i][j]=8*cnt[i-1][j]+(j>0?cnt[i-1][j-1]:0)+cnt0[i-1][j];
    pows[0]=1;
    pows2[0]=1;
    for(p=1,i=0;b/p>=10;p*=10,++i,pows[i]=p,pows2[i]=pows2[i-1]*9);
}

long long comp(long long x,long long type,long long pow,long long many)
{
    if(x<0)
        return 0;
    if(many<0)
        return 0;
    if(pow==0)
        if(x>=c)
            if(many==1)
                return 1;
            else
                if(many==0)
                    return x;
                else
                    return 0;
        else
            if(many==0)
                return x+1;
            else
                return 0;
    if(type==0)
        if(c==0)
            if(x<pows[pow])
                return comp(x,type,pow-1,many);
            else
                return (x/pows[pow]-1)*cnt[pow][many]+cnt0[pow][many]+comp(x%pows[pow],1,pow-1,many);
        else
        {
            if(x<pows[pow])
                return comp(x,type,pow-1,many);
            long long y=(x/pows[pow]-1);
            if(y>=c)
                return (y-1)*cnt[pow][many]+cnt[pow][many-1]+cnt0[pow][many]+comp(x%pows[pow],1,pow-1,many);
            if(y+1==c)
                return y*cnt[pow][many]+cnt0[pow][many]+comp(x%pows[pow],1,pow-1,many-1);
            return y*cnt[pow][many]+cnt0[pow][many]+comp(x%pows[pow],1,pow-1,many);
        }
    long long y=x/pows[pow];
    if(y==0)
            return comp(x,type,pow-1,many-(c==0));
    --y;
    if(y>=c)
        return y*cnt[pow][many]+cnt[pow][many-1]+comp(x%pows[pow],1,pow-1,many);
    if(y+1==c)
        return (y+1)*cnt[pow][many]+comp(x%pows[pow],1,pow-1,many-1);
    return (y+1)*cnt[pow][many]+comp(x%pows[pow],1,pow-1,many);
}

long double aux;

int main()
{
    f>>a>>b>>c>>k;
    precalc();
    for(j=k;j<=i+1;++j)
    {
        rez+=comp(b,0,i,j);
        rez-=comp(a-1,0,i,j);
    }
    g.setf(ios::fixed,ios::floatfield);
    g.precision(4);
    aux=rez*1.0;
    aux/=1.0*(b-a+1);
    g<<aux<<"\n";

    f.close();
    g.close();

    return 0;
}