結果

問題 No.2709 1975 Powers
ユーザー Basin-Bug
提出日時 2024-03-31 14:23:41
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
MLE  
実行時間 -
コード長 1,378 bytes
コンパイル時間 807 ms
コンパイル使用メモリ 79,616 KB
実行使用メモリ 718,976 KB
最終ジャッジ日時 2024-09-30 19:28:01
合計ジャッジ時間 18,338 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 5 WA * 12 MLE * 8
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <vector>
#include <cmath>
using namespace std;

long long pow(long long x, long long n, long long P) {
    long long ret = 1;
    while (n > 0) {
        if (n & 1) ret = ret * x % P;  // n の最下位bitが 1 ならば x^(2^i) をかける
        x = x * x % P;
        n >>= 1;  // n を1bit 左にずらす
    }
    return ret;
}

int main() {
    long long N, P, Q;
    cin >> N >> P >> Q;
    vector<long long> A(N);
    for (long long &a : A) cin >> a;

    vector<vector<vector<long long>>> dp(N + 1, vector<vector<long long>>(5, vector<long long>(P, 0)));
    dp[0][0][0] = 1;

    for (long long i = 1; i <= N; ++i) {
        long long a = pow(10, A[i - 1], P);
        long long b = pow(9, A[i - 1], P);
        long long c = pow(7, A[i - 1], P);
        long long d = pow(5, A[i - 1], P);

        for (long long j = 0; j < P; ++j) {
            dp[i][0][j] += dp[i - 1][0][j];
            dp[i][1][(j + a) % P] += dp[i - 1][0][j];
            dp[i][1][j] += dp[i - 1][1][j];
            dp[i][2][(j + b) % P] += dp[i - 1][1][j];
            dp[i][2][j] += dp[i - 1][2][j];
            dp[i][3][(j + c) % P] += dp[i - 1][2][j];
            dp[i][3][j] += dp[i - 1][3][j];
            dp[i][4][(j + d) % P] += dp[i - 1][3][j];
            dp[i][4][j] += dp[i - 1][4][j];
        }
    }

    cout << dp[N][4][Q] << endl;
    return 0;
}
0