#196892 (C++14) No.493 とても長い数列と文字列(Long Long Sequence and a String）

提出ソース

結果

問題	No.493 とても長い数列と文字列(Long Long Sequence and a String）
ユーザー	Pachicobue
提出日時	2017-08-13 21:32:07
言語	C++14 (gcc 12.3.0 + boost 1.83.0)
結果	AC
実行時間	2 ms / 800 ms
コード長	3,599 bytes
コンパイル時間	1,595 ms
コンパイル使用メモリ	174,748 KB
実行使用メモリ	5,248 KB
最終ジャッジ日時	2024-10-13 09:15:44
合計ジャッジ時間	4,121 ms
ジャッジサーバーID （参考情報）	judge3 / judge4

このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	2 ms 5,248 KB
testcase_01	AC	2 ms 5,248 KB
testcase_02	AC	2 ms 5,248 KB
testcase_03	AC	2 ms 5,248 KB
testcase_04	AC	2 ms 5,248 KB
testcase_05	AC	1 ms 5,248 KB
testcase_06	AC	2 ms 5,248 KB
testcase_07	AC	1 ms 5,248 KB
testcase_08	AC	2 ms 5,248 KB
testcase_09	AC	1 ms 5,248 KB
testcase_10	AC	2 ms 5,248 KB
testcase_11	AC	2 ms 5,248 KB
testcase_12	AC	2 ms 5,248 KB
testcase_13	AC	1 ms 5,248 KB
testcase_14	AC	1 ms 5,248 KB
testcase_15	AC	2 ms 5,248 KB
testcase_16	AC	2 ms 5,248 KB
testcase_17	AC	2 ms 5,248 KB
testcase_18	AC	2 ms 5,248 KB
testcase_19	AC	2 ms 5,248 KB
testcase_20	AC	1 ms 5,248 KB
testcase_21	AC	2 ms 5,248 KB
testcase_22	AC	2 ms 5,248 KB
testcase_23	AC	2 ms 5,248 KB
testcase_24	AC	2 ms 5,248 KB
testcase_25	AC	1 ms 5,248 KB
testcase_26	AC	2 ms 5,248 KB
testcase_27	AC	1 ms 5,248 KB
testcase_28	AC	1 ms 5,248 KB
testcase_29	AC	1 ms 5,248 KB
testcase_30	AC	2 ms 5,248 KB
testcase_31	AC	2 ms 5,248 KB
testcase_32	AC	2 ms 5,248 KB
testcase_33	AC	2 ms 5,248 KB
testcase_34	AC	2 ms 5,248 KB
testcase_35	AC	1 ms 5,248 KB
testcase_36	AC	2 ms 5,248 KB
testcase_37	AC	2 ms 5,248 KB
testcase_38	AC	2 ms 5,248 KB
testcase_39	AC	1 ms 5,248 KB
testcase_40	AC	2 ms 5,248 KB
testcase_41	AC	2 ms 5,248 KB
testcase_42	AC	1 ms 5,248 KB
testcase_43	AC	1 ms 5,248 KB
testcase_44	AC	1 ms 5,248 KB
testcase_45	AC	2 ms 5,248 KB
testcase_46	AC	1 ms 5,248 KB
testcase_47	AC	1 ms 5,248 KB
testcase_48	AC	1 ms 5,248 KB
testcase_49	AC	2 ms 5,248 KB
testcase_50	AC	2 ms 5,248 KB
testcase_51	AC	2 ms 5,248 KB
testcase_52	AC	2 ms 5,248 KB
testcase_53	AC	1 ms 5,248 KB
testcase_54	AC	1 ms 5,248 KB
testcase_55	AC	1 ms 5,248 KB
testcase_56	AC	2 ms 5,248 KB
testcase_57	AC	1 ms 5,248 KB
testcase_58	AC	1 ms 5,248 KB
testcase_59	AC	2 ms 5,248 KB
testcase_60	AC	2 ms 5,248 KB
testcase_61	AC	2 ms 5,248 KB
testcase_62	AC	2 ms 5,248 KB
testcase_63	AC	2 ms 5,248 KB
testcase_64	AC	2 ms 5,248 KB
testcase_65	AC	1 ms 5,248 KB
testcase_66	AC	2 ms 5,248 KB
testcase_67	AC	2 ms 5,248 KB
testcase_68	AC	2 ms 5,248 KB
testcase_69	AC	2 ms 5,248 KB
testcase_70	AC	1 ms 5,248 KB
testcase_71	AC	1 ms 5,248 KB
testcase_72	AC	1 ms 5,248 KB
testcase_73	AC	1 ms 5,248 KB
testcase_74	AC	1 ms 5,248 KB
testcase_75	AC	2 ms 5,248 KB
testcase_76	AC	2 ms 5,248 KB
testcase_77	AC	2 ms 5,248 KB
testcase_78	AC	2 ms 5,248 KB
testcase_79	AC	2 ms 5,248 KB
testcase_80	AC	2 ms 5,248 KB
testcase_81	AC	2 ms 5,248 KB
testcase_82	AC	2 ms 5,248 KB
testcase_83	AC	2 ms 5,248 KB
testcase_84	AC	2 ms 5,248 KB
testcase_85	AC	2 ms 5,248 KB
testcase_86	AC	2 ms 5,248 KB
testcase_87	AC	2 ms 5,248 KB
testcase_88	AC	2 ms 5,248 KB
testcase_89	AC	2 ms 5,248 KB
testcase_90	AC	1 ms 5,248 KB
testcase_91	AC	2 ms 5,248 KB
testcase_92	AC	2 ms 5,248 KB
testcase_93	AC	2 ms 5,248 KB
testcase_94	AC	2 ms 5,248 KB
testcase_95	AC	2 ms 5,248 KB
testcase_96	AC	2 ms 5,248 KB
testcase_97	AC	1 ms 5,248 KB
testcase_98	AC	2 ms 5,248 KB
testcase_99	AC	1 ms 5,248 KB
testcase_100	AC	1 ms 5,248 KB
testcase_101	AC	1 ms 5,248 KB
testcase_102	AC	1 ms 5,248 KB
testcase_103	AC	1 ms 5,248 KB
testcase_104	AC	2 ms 5,248 KB
testcase_105	AC	2 ms 5,248 KB
testcase_106	AC	2 ms 5,248 KB
testcase_107	AC	2 ms 5,248 KB
testcase_108	AC	2 ms 5,248 KB
testcase_109	AC	2 ms 5,248 KB
testcase_110	AC	2 ms 5,248 KB
testcase_111	AC	2 ms 5,248 KB
testcase_112	AC	1 ms 5,248 KB
testcase_113	AC	1 ms 5,248 KB
testcase_114	AC	2 ms 5,248 KB
testcase_115	AC	2 ms 5,248 KB
testcase_116	AC	2 ms 5,248 KB
testcase_117	AC	1 ms 5,248 KB
testcase_118	AC	2 ms 5,248 KB

権限があれば一括ダウンロードができます

ソースコード

raw source code

// {{{ Templates
#include <bits/stdc++.h>

#define show(x) cerr << #x << " = " << x << endl

using namespace std;
using ll = long long;
using pii = pair<int, int>;
using vi = vector<int>;

template <typename T>
ostream& operator<<(ostream& os, const vector<T>& v)
{
    os << "sz:" << v.size() << "\n[";
    for (const auto& p : v) {
        os << p << ",";
    }
    os << "]\n";
    return os;
}

template <typename S, typename T>
ostream& operator<<(ostream& os, const pair<S, T>& p)
{
    os << "(" << p.first << "," << p.second
       << ")";
    return os;
}


constexpr ll MOD = (ll)1e9 + 7LL;

template <typename T>
constexpr T INF = numeric_limits<T>::max() / 100;

// }}}

using pll = pair<ll, ll>;

inline ll convert(const char c)
{
    return c == '0' ? 10LL : (ll)(c - '0');
}

inline pll calcst(const string& s, const ll l, const ll r)
{
    assert(0 <= l);
    assert(l < r);
    assert(r <= s.size());
    ll sum = 0;
    ll prod = 1;
    for (ll i = l; i < r; i++) {
        const char c = s[i];
        sum = (sum + convert(c));
        prod = (prod * convert(c)) % MOD;
    }
    return make_pair(sum, prod);
}

pll rec(const ll n, const ll l, const ll r, const vector<ll>& length, const vector<pll>& calced)  // f(n)の[L,R)の(和,積)
{
    assert(l >= 0);
    assert(l < r);
    if (n == 0) {
        return make_pair(0, 1);
    }
    if (l == 0 and r == length[n]) {
        return calced[n];
    } else {
        const string intst = to_string(n * n);
        const ll interval = intst.size();
        const ll half = (length[n] - interval) / 2;
        pll ret = make_pair(0, 1);
        if (l < half) {
            const pll p = rec(n - 1, l, min(r, half), length, calced);
            ret.first = (ret.first + p.first);
            ret.second = (ret.second * p.second) % MOD;
        }
        if (r > half + interval) {
            const pll p = rec(n - 1, max(l - half - interval, 0LL), r - half - interval, length, calced);
            ret.first = (ret.first + p.first);
            ret.second = (ret.second * p.second) % MOD;
        }
        if (r > half and l < half + interval) {
            const pll p = calcst(intst, max(l - half, 0LL), min(interval, r - half));
            ret.first = (ret.first + p.first);
            ret.second = (ret.second * p.second) % MOD;
        }
        return ret;
    }
}

inline pll cal(const string& s)
{
    ll sum = 0;
    ll prod = 1;
    for (const char c : s) {
        sum += convert(c);
        sum = sum;
        prod *= convert(c);
        prod = prod % MOD;
    }
    return make_pair(sum, prod);
}

int main()
{
    cin.tie(0);
    ios::sync_with_stdio(false);

    ll k, l, r;
    cin >> k >> l >> r;
    l--, r--;

    constexpr ll STOP = (ll)2e18;
    vector<ll> length(1, 0);
    ll len = 0;
    for (ll i = 1;; i++) {
        len = len * 2 + (ll)log10(i * i) + 1;
        length.push_back(len);
        if (len > STOP) {
            break;
        }
    }
    const ll size = length.size();
    if (k >= size) {
        k = size - 1;
    }
    vector<pll> calced(size, make_pair(0, 1));
    for (ll i = 1; i < size; i++) {
        const string s = to_string(i * i);
        const pll p = cal(s);
        calced[i].first = (calced[i - 1].first * 2) + p.first;
        calced[i].second = ((calced[i - 1].second * calced[i - 1].second) % MOD) * p.second;
        calced[i].second = calced[i].second % MOD;
    }

    if (length[k] <= r) {
        cout << -1 << endl;
        return 0;
    }
    const pll ans = rec(k, l, r + 1, length, calced);
    cout << ans.first << " " << ans.second << endl;
    return 0;
}

yukicoder

結果

テストケース

ソースコード