結果

問題 No.28 末尾最適化
ユーザー nebukuro09nebukuro09
提出日時 2017-06-01 11:24:46
言語 D
(dmd 2.106.1)
結果
WA  
実行時間 -
コード長 928 bytes
コンパイル時間 923 ms
コンパイル使用メモリ 109,356 KB
実行使用メモリ 6,752 KB
最終ジャッジ日時 2023-09-03 13:48:25
合計ジャッジ時間 1,820 ms
ジャッジサーバーID
（参考情報） 		judge15 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,380 KB
testcase_01 WA -
権限があれば一括ダウンロードができます

ソースコード

diff # 

import std.stdio, std.array, std.string, std.conv, std.algorithm;
import std.typecons, std.range, std.random, std.math, std.container;
import std.numeric, std.bigint, core.bitop;

immutable long MOD = 100000009;
long[] primes = [2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31];

void solve() {
    auto s = readln.split.map!(to!long);
    auto seed = s[0];
    auto N = s[1].to!int;
    auto K = s[2].to!int;
    auto B = s[3];

    auto A = new long[](N+1);
    A[0] = seed;
    foreach (i; 0..N) {
        A[i+1] = 1 + ((A[i] * A[i] % MOD + A[i] * 12345 % MOD)) % MOD;
    }

    long bp;
    foreach (p; primes) if (B % p == 0) bp = p;

    auto cnt = new long[](N+1);
    fill(cnt, 0);
    foreach (i; 0..N+1) while (A[i] % bp == 0) cnt[i]++, A[i]/=bp;

    cnt.sort();
    long d = 0;
    while (B % bp == 0) d++, B/=bp;
    writeln( cnt[0..K].sum / d );
}

void main() {
    auto Q = readln.chomp.to!int;
    while (Q--) solve;
}
0