結果

問題 No.990 N×Mマス計算（Kの倍数）
コンテスト
ユーザー kcvlex
提出日時 2020-02-14 22:40:41
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 142 ms / 2,000 ms
コード長 3,107 bytes
コンパイル時間 2,018 ms
コンパイル使用メモリ 183,740 KB
実行使用メモリ 17,280 KB
最終ジャッジ日時 2024-11-16 01:03:49
合計ジャッジ時間 3,542 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 19
権限があれば一括ダウンロードができます

ソースコード

diff # 

// #define DEBUGGING
#include <bits/stdc++.h>
#define endl '\n'
#define ALL(V) (V).begin(), (V).end()
#define ALLR(V) (V).rbegin(), (V).rend()
using ll = std::int64_t;
using ull = std::uint64_t;
using PLL = std::pair<ll, ll>;
using TLL = std::tuple<ll, ll, ll>;
template <typename T> using V = std::vector<T>;
template <typename T> using VV = V<V<T>>;
template <typename T> const T& var_min(const T &t) { return t; }
template <typename T> const T& var_max(const T &t) { return t; }
template <typename T, typename... Tail> const T& var_min(const T &t, const Tail&... tail) { return std::min(t, var_min(tail...)); }
template <typename T, typename... Tail> const T& var_max(const T &t, const Tail&... tail) { return std::max(t, var_max(tail...)); }
template <typename T, typename... Tail> void chmin(T &t, const Tail&... tail) { t = var_min(t, tail...); }
template <typename T, typename... Tail> void chmax(T &t, const Tail&... tail) { t = var_max(t, tail...); }
template <typename T> const T& clamp(const T &t, const T &low, const T &high) { return std::max(low, std::min(high, t)); }
template <typename T> void chclamp(T &t, const T &low, const T &high) { return t = clamp(t, low, high); }
namespace init__ { struct InitIO { InitIO() { std::cin.tie(nullptr); std::ios_base::sync_with_stdio(false); std::cout << std::fixed << std::setprecision(30); } } init_io; }
#define mv_rec make_v(init, tail...)
template <typename T> T make_v(T init) { return init; }
template <typename T, typename... Tail> auto make_v(T init, size_t s, Tail... tail) { return V<decltype(mv_rec)>(s, mv_rec); }
#undef mv_rec
using namespace std;

#ifdef DEBUGGING
#include "../../debug/debug.cpp"
#else
#define DEBUG(...) 0
#define DEBUG_SEPARATOR_LINE 0
#endif

ll gcd(ll a, ll b) { return b ? gcd(b, a % b) : a; }

ll solve_p(ll N, ll M, ll K, const V<ll> &A, const V<ll> &B) {
    map<ll, ll> cnt;
    for (ll e : A) cnt[e % K]++;
    ll ret = 0;
    for (ll e : B) {
        ll mod = e % K;
        ll key = (K - mod) % K;
        ret += cnt[key];
    }
    return ret;
}

ll solve_m(ll N, ll M, ll K, const V<ll> &A, const V<ll> &B) {
    map<ll, ll> gcd_v;
    for (ll e : B) gcd_v[gcd(e, K)]++;

    V<ll> divs;
    for (ll i = 1; i * i <= K; i++) {
        if (K % i) continue;
        divs.push_back(i);
        if (i * i < K) divs.push_back(K / i);
    }

    sort(ALL(divs));
    DEBUG(divs);
    map<ll, ll> cnt;
    for (ll i = 0; i < divs.size(); i++) {
        ll c = gcd_v[divs[i]];
        DEBUG(make_tuple(c, divs[i]));
        for (ll j = 0; j <= i; j++) if (divs[i] % divs[j] == 0) {
            DEBUG(make_tuple(divs[i], divs[j]));
            cnt[divs[j]] += c;
        }
    }
    
    for (auto &&e : cnt) DEBUG(e);
    ll ret = 0;
    for (ll e : A) {
        ll g = gcd(e, K);
        ret += cnt[K / g];
    }
    return ret;
}

int main() {
    ll N, M, K;
    cin >> N >> M >> K;
    char op;
    cin >> op;
    V<ll> B(M), A(N);
    for (ll &e : B) cin >> e;
    for (ll &e : A) cin >> e;
    ll ans = (op == '+' ?  solve_p(N, M, K, A, B) : solve_m(N, M, K, A, B));
    cout << ans << endl;
    return 0;
}
0