結果

問題 No.493 とても長い数列と文字列(Long Long Sequence and a String）
ユーザー e869120
提出日時 2017-03-08 20:40:27
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 3 ms / 800 ms
コード長 2,276 bytes
コンパイル時間 842 ms
コンパイル使用メモリ 72,056 KB
実行使用メモリ 6,948 KB
最終ジャッジ日時 2024-06-23 19:32:50
合計ジャッジ時間 3,430 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 115
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

#include<iostream>
#include<algorithm>
#include<string>
using namespace std;
long long K, L, R, U, V; string Y;
long long sums(long long P, long long Q) {
    if (P >= (1LL << Q))return (1LL << 62);
    if (P < 0)return -(1LL << 62);
    long long ret = 0, U = P;
    for (int i = 1; i <= Q; i++) {
        ret += ((U + 1) / 2) * to_string(i*i).size();
        U -= (U + 1) / 2;
    }
    return ret;
}
long long modpow(long long a, long long b, long long m) {
    long long p = 1, q = a;
    for (int i = 0; i < 63; i++) {
        if ((b / (1LL << i)) % 2 == 1)p *= q;
        q *= q;
        p %= m; q %= m;
    }
    return p;
}
void search() {
    //左の境界を探索
    long long L1 = -1, R1 = (1LL << K) + 1, M1;
    while (true) {
        M1 = (L1 + R1) / 2;
        long long Y1 = sums(M1 - 1, K), Y2 = sums(M1, K);
        if (Y1 < L && Y2 >= L) {
            U = M1; U++;
            if (Y2 - L >= 1) {
                long long I = M1, J = 1;
                while (I % 2 == 0) { I /= 2; J++; }
                Y += to_string(J*J).substr(to_string(J*J).size() - (Y2 - L), (Y2 - L));
            }
            break;
        }
        if (Y2 < L)L1 = M1;
        if (Y1 >= L)R1 = M1;
    }
    //右の境界を探索
    L1 = -1; R1 = (1LL << K) + 1;
    while (true) {
        M1 = (L1 + R1) / 2;
        long long Y1 = sums(M1, K), Y2 = sums(M1 + 1, K);
        if (Y1 <= R && Y2 > R) {
            V = M1;
            if (R - Y1 >= 1) {
                long long I = M1, J = 1; I++;
                while (I % 2 == 0) { I /= 2; J++; }
                Y += to_string(J*J).substr(0, (R - Y1));
            }
            break;
        }
        if (Y2 <= R)L1 = M1;
        if (Y1 > R)R1 = M1;
    }
}
int main() {
    cin >> K >> L >> R; L--; if (K >= 62)K = 62;
    if (R > sums((1LL << K) - 1, K)) { cout << "-1" << endl; return 0; }
    search();
    //端数部分の和と積を求める
    long long wa = 0, seki = 1;
    for (int i = 0; i < Y.size(); i++) {
        long long G = (Y[i] - '0'); if (G == 0)G = 10;
        wa += G; seki *= G;
    }
    //区間の部分の積と和を求める
    U--;
    for (int i = 1; i <= K; i++) {
        long long U1 = (U + 1) / 2;
        long long U2 = (V + 1) / 2;
        string T = to_string(i*i);
        long long w1 = 0, w2 = 1;
        for (int j = 0; j < T.size(); j++) {
            long long G = (T[j] - '0'); if (G == 0)G = 10;
            w1 += G; w2 *= G;
        }
        //cout << U2 - U1 << ' ' << U << ' ' << V << endl;
        wa += w1*(U2 - U1);
        seki *= modpow(w2, U2 - U1, 1000000007);
        seki %= 1000000007;
        U -= U1; V -= U2;
    }
    //出力
    cout << wa << ' ' << seki << endl;
    return 0;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0