結果

問題 No.493 とても長い数列と文字列(Long Long Sequence and a String）
ユーザー latte0119
提出日時 2017-03-10 23:17:32
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 2 ms / 800 ms
コード長 1,906 bytes
コンパイル時間 1,597 ms
コンパイル使用メモリ 164,260 KB
実行使用メモリ 6,948 KB
最終ジャッジ日時 2024-06-24 08:39:44
合計ジャッジ時間 4,056 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 115
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
using namespace std;

#define int long long
typedef pair<int,int>pint;
typedef vector<int>vint;
typedef vector<pint>vpint;
#define pb push_back
#define mp make_pair
#define fi first
#define se second
#define all(v) (v).begin(),(v).end()
#define rep(i,n) for(int i=0;i<(n);i++)
#define reps(i,f,n) for(int i=(f);i<(n);i++)
#define each(it,v) for(__typeof((v).begin()) it=(v).begin();it!=(v).end();it++)
template<class T,class U>void chmin(T &t,U f){if(t>f)t=f;}
template<class T,class U>void chmax(T &t,U f){if(t<f)t=f;}

const int mod=1000000007;
const int lim=60;

int mpow(int n,int m){
    int ret=1;
    while(m){
        if(m&1)ret=ret*n%mod;
        n=n*n%mod;
        m>>=1;
    }
    return ret;
}

int l[111111];
int s[111111];
int m[111111];


pint solve(int r,int k){
    if(r==l[k])return pint(s[k],m[k]);
    if(r<l[k-1]){
        return solve(r,k-1);
    }
    pint ret=pint(s[k-1],m[k-1]);
    r-=l[k-1];
    stringstream ss;ss<<k*k;
    string str;ss>>str;
    for(int i=0;r&&i<str.size();i++,r--){
        int d=str[i]-'0';
        if(d==0)d=10;
        ret.fi+=d;
        ret.se=ret.se*d%mod;
    }
    if(r==0)return ret;

    pint latte=solve(r,k-1);
    ret.fi+=latte.fi;
    ret.se=ret.se*latte.se%mod;
    return ret;
}

signed main(){
    l[1]=s[1]=m[1]=1;
    m[0]=1;
    for(int i=2;i<100;i++){
        stringstream ss;ss<<i*i;
        string str;ss>>str;
        l[i]=l[i-1]*2+str.size();
        s[i]=s[i-1]*2;
        m[i]=m[i-1]*m[i-1]%mod;
        for(int j=0;j<str.size();j++){
            int d=str[j]-'0';
            if(d==0)d=10;
            s[i]+=d;
            m[i]=m[i]*d%mod;
        }
    }

    int K,L,R;
    cin>>K>>L>>R;L--;
    if(K<lim&&l[K]<R){
        cout<<-1<<endl;
        return 0;
    }

    pint x,y;
    x=solve(L,lim);
    y=solve(R,lim);

    cout<<y.fi-x.fi<<" ";
    cout<<y.se*mpow(x.se,mod-2)%mod<<endl;
    return 0;
}
0