結果

問題 No.2709 1975 Powers
コンテスト
ユーザー Basin-Bug
提出日時 2024-03-31 14:21:56
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 1,306 bytes
コンパイル時間 924 ms
コンパイル使用メモリ 90,164 KB
最終ジャッジ日時 2025-02-20 17:11:10
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 9 WA * 16
権限があれば一括ダウンロードができます

ソースコード

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() {
    int N, P, Q;
    cin >> N >> P >> Q;
    vector<int> A(N);
    for (int &a : A) cin >> a;

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

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

        for (int 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