結果

問題 No.990 N×Mマス計算（Kの倍数）
ユーザー SSRS
提出日時 2021-12-17 03:11:27
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 210 ms / 2,000 ms
コード長 1,500 bytes
コンパイル時間 2,032 ms
コンパイル使用メモリ 181,492 KB
実行使用メモリ 18,048 KB
最終ジャッジ日時 2024-09-14 03:19:06
合計ジャッジ時間 3,720 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 19
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
int gcd(int a, int b){
  if (b == 0){
    return a;
  } else {
    return gcd(b, a % b);
  }
}
int main(){
  int N, M, K;
  cin >> N >> M >> K;
  char op;
  cin >> op;
  vector<int> B(M);
  for (int i = 0; i < M; i++){
    cin >> B[i];
  }
  vector<int> A(N);
  for (int i = 0; i < N; i++){
    cin >> A[i];
  }
  if (op == '+'){
    map<int, int> cntA;
    for (int i = 0; i < N; i++){
      cntA[A[i] % K]++;
    }
    map<int, int> cntB;
    for (int i = 0; i < M; i++){
      cntB[B[i] % K]++;
    }
    long long ans = 0;
    for (auto P : cntA){
      ans += (long long) P.second * cntB[(K - P.first) % K];
    }
    cout << ans << endl;
  }
  if (op == '*'){
    vector<int> F;
    for (int i = 1; i * i <= K; i++){
      if (K % i == 0){
        F.push_back(i);
        if (i * i < K){
          F.push_back(K / i);
        }
      }
    }
    sort(F.begin(), F.end());
    int c = F.size();
    vector<int> cntA(c, 0);
    for (int i = 0; i < N; i++){
      int p = lower_bound(F.begin(), F.end(), gcd(K, A[i])) - F.begin();
      cntA[p]++;
    }
    vector<int> cntB(c, 0);
    for (int i = 0; i < M; i++){
      int p = lower_bound(F.begin(), F.end(), gcd(K, B[i])) - F.begin();
      cntB[p]++;
    }
    long long ans = 0;
    for (int i = 0; i < c; i++){
      for (int j = 0; j < c; j++){
        if ((long long) F[i] * F[j] % K == 0){
          ans += (long long) cntA[i] * cntB[j];
        }
      }
    }
    cout << ans << endl;
  }
}
0