結果

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

ソースコード

diff #
raw source code 

#include <bits/stdc++.h>
using namespace std;
template <typename T> using min_priority_queue = priority_queue<T,vector<T>,greater<T>>;

template<typename T>
void printv(vector<T> &v){
    bool b = false;
    for(auto i : v){
        if(b) cout << " ";
        else b = true;
        cout << i;
    }
    cout << endl;
}

template <typename T>
bool chmin(T &a, const T& b) {
    if (a > b) {
        a = b;  // aをbで更新
        return true;
    }
    return false;
}

template <typename T>
bool chmax(T &a, const T& b) {
    if (a < b) {
        a = b;  // aをbで更新
        return true;
    }
    return false;
}

int64_t pow_mod(int64_t a, int64_t b, int64_t n){
    a %= n;
    if(a == 0) return 0;
    if(a == 1 || b == 0) return 1;
    if(b == 1) return a;
    int64_t c = b/2;
    int64_t d = b%2;
    int64_t ac = pow_mod(a,c,n);
    int64_t ad = pow_mod(a,d,n);
    return ac*ac%n*ad%n;
}

void yn(bool b){
    if(b) cout << "Yes" << endl;
    else cout << "No" << endl;
}

bool debug;
bool randomInput;
bool debugOutput;
int numOfTestCase;

using ans_type = int;

void input(){
    if(numOfTestCase > 1){
        ;
    }

    if(randomInput){
        
    }
    else{
        
    }
    return;
}

void output_input(){
    ;
}

ans_type calc(){
    int N, P, Q;
    cin >> N >> P >> Q;
    vector<int64_t> A(N);
    for(auto &a : A) cin >> a;
    sort(A.begin(), A.end());

    int M = 4;
    vector<int> B = {10, 9, 7, 5};

    vector<vector<int64_t>> V(N, vector<int64_t>(M));
    for(int i = 0; i < N; i++){
        for(int j = 0; j < M; j++){
            V[i][j] = pow_mod(B[j], A[i], P);
        }
    }

    vector<vector<int>> W(N+1, vector<int>(P, 0));//[i, N)
    for(int i = N-1; i >= 0; i--){
        W[i] = W[i+1];
        W[i][V[i][M-1]]++;
    }

    int64_t ans = 0;
    for(int i = 0; i < N; i++){
        for(int j = i+1; j < N; j++){
            for(int k = j+1; k < N; k++){
                int64_t r = Q - V[i][0] - V[j][1] - V[k][2];
                r = (r%P+P)%P;
                ans += W[k+1][r];
            }
        }
    }
    cout << ans << endl;
    return ans_type();
}

ans_type calc_simple(){
    return ans_type();
}

void output(ans_type ans){
    return;
}

int main(){
    debug = 0;
    randomInput = 0;
    debugOutput = 0;
    numOfTestCase = 1;
    srand(time(NULL));
    cout << fixed << setprecision(12);
    
    if(numOfTestCase == 0) cin >> numOfTestCase;

    if(debug){
        for(int i = 0; i < numOfTestCase; i++){
            input();
            ans_type ans = calc();
            ans_type ansSimple = calc_simple();
            if(ans != ansSimple){
                output_input();
                output(ans);
                output(ansSimple);
            }
        }
    }
    else{
        for(int i = 0; i < numOfTestCase; i++){
            input();
            output(calc());
        }
    }

    return 0;
}
0