結果

問題 No.990 N×Mマス計算(Kの倍数)
ユーザー NOSS
提出日時 2020-02-14 22:12:41
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 195 ms / 2,000 ms
コード長 1,790 bytes
コンパイル時間 1,933 ms
コンパイル使用メモリ 182,960 KB
実行使用メモリ 17,280 KB
最終ジャッジ日時 2024-11-16 00:47:48
合計ジャッジ時間 3,692 ms
ジャッジサーバーID
(参考情報)
judge5 / judge2
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 2
other AC * 19
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <bits/stdc++.h>
using namespace std;

using ll = long long int;
using P = pair<int, int>;
using P3 = pair<int,P>;
using PP = pair<P, P>;
constexpr int INF = 1<<30;
constexpr ll MOD = ll(1e9)+7;
constexpr int di[] = {0,1,0,-1};
constexpr int dj[] = {1,0,-1,0};

vector<ll> divisor(ll n) {
    vector<ll> res;
    for (ll i = 1; i * i <= n; i++) {
        if (n % i == 0) {
            res.push_back(i);
            if (i != n / i) res.push_back(n / i);
        }
    }
    return res;
}

ll gcd(ll p, ll q){
    if(q==0) return p;
    else return gcd(q, p%q);
}

int main(){
    ll n, m, k;
    char op;
    cin >> n >> m >> k >> op;
    vector<ll> a(n), b(m);
    for(int i=0;i<m;i++){
        cin >> b[i];
    }
    for(int i=0;i<n;i++){
        cin >> a[i];
    }
    if(op == '+'){
        map<ll,ll> mp;
        for(int i=0;i<m;i++){
            mp[b[i]%k]++;
        }
        ll cnt = 0;
        for(int i=0;i<n;i++){
            cnt += mp[(k-a[i]%k)%k];
        }
        cout << cnt << endl;
    }
    else{
        auto div = divisor(k);
        sort(div.begin(), div.end());
        map<ll,ll> mp;
        vector<ll> cum(div.size());
        for(int i=0;i<int(div.size());i++){
            mp[div[i]] = i;
        }
        for(int i=0;i<m;i++){
            ll g = gcd(b[i],k);
            cum[mp[g]]++;
        }
        for(int i=0;i<int(div.size());i++){
            ll d = div[i];
            ll tmp = 0;
            for(int j=i;j<int(div.size());j++){
                if(div[j]%d == 0){
                    tmp += cum[j];
                }
            }
            cum[i] = tmp;
        }
        ll cnt = 0;
        for(int i=0;i<n;i++){
            ll g = gcd(a[i],k);
            cnt += cum[mp[k/g]];
        }
        cout << cnt << endl;
    }
    return 0;
}
0