結果

問題 No.206 数の積集合を求めるクエリ
ユーザー cormorancormoran
提出日時 2016-08-16 23:14:30
言語 C++14
(gcc 13.2.0 + boost 1.83.0)
結果
AC  
実行時間 369 ms / 7,000 ms
コード長 2,534 bytes
コンパイル時間 1,284 ms
コンパイル使用メモリ 90,552 KB
実行使用メモリ 28,920 KB
最終ジャッジ日時 2024-04-25 08:10:20
合計ジャッジ時間 6,364 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 2 ms
5,376 KB
testcase_02 AC 2 ms
5,376 KB
testcase_03 AC 2 ms
5,376 KB
testcase_04 AC 2 ms
5,376 KB
testcase_05 AC 2 ms
5,376 KB
testcase_06 AC 8 ms
5,376 KB
testcase_07 AC 7 ms
5,376 KB
testcase_08 AC 8 ms
5,376 KB
testcase_09 AC 8 ms
5,376 KB
testcase_10 AC 1 ms
5,376 KB
testcase_11 AC 2 ms
5,376 KB
testcase_12 AC 11 ms
5,376 KB
testcase_13 AC 12 ms
5,376 KB
testcase_14 AC 13 ms
5,376 KB
testcase_15 AC 11 ms
5,376 KB
testcase_16 AC 11 ms
5,376 KB
testcase_17 AC 252 ms
28,916 KB
testcase_18 AC 244 ms
28,788 KB
testcase_19 AC 243 ms
28,788 KB
testcase_20 AC 235 ms
28,788 KB
testcase_21 AC 235 ms
28,788 KB
testcase_22 AC 237 ms
28,788 KB
testcase_23 AC 242 ms
28,916 KB
testcase_24 AC 364 ms
28,788 KB
testcase_25 AC 363 ms
28,920 KB
testcase_26 AC 345 ms
28,788 KB
testcase_27 AC 316 ms
28,896 KB
testcase_28 AC 350 ms
28,912 KB
testcase_29 AC 369 ms
28,656 KB
testcase_30 AC 341 ms
28,784 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<iostream>
#include<vector>
#include<cassert>
#include<algorithm>
#include<cstdlib>
#include<cmath>
#include<complex>
using namespace std;

using ll = long long;
#define rep(i, j) for(int i=0; i < (int)(j); i++)
#define repeat(i, j, k) for(int i = (j); i < (int)(k); i++)
#define all(v) v.begin(),v.end()
// vector
template<class T> istream& operator >> (istream &is , vector<T> &v) {
    for(T &a : v) is >> a; return is;
}

void fft_base(vector<complex<double>> &f, int sign) {
    int n = f.size();
    if (n == 1) return;
    assert(n == pow(2, log2(n)));  // f.size() must be 2^i

    vector<complex<double>> f0(n / 2), f1(n / 2);
    for (int i = 0; i < n / 2; i++) {
        f0[i] = f[2 * i + 0];
        f1[i] = f[2 * i + 1];
    }
    fft_base(f0, sign);
    fft_base(f1, sign);

    // zeta := exp{2pi * sqrt(-1) / n}
    const complex<double> zeta = polar(1.0, sign * 2 * M_PI / n);
    complex<double> pow_zeta   = 1;

    for (int i = 0; i < n; i++) {
        f[i] = f0[i % (n / 2)] + pow_zeta * f1[i % (n / 2)];
        pow_zeta *= zeta;
    }
    return;
}

vector<complex<double>> fft(vector<complex<double>> f) {
    int n = f.size();
    n = pow(2, ceil(log2(n)));
    f.resize(n, 0);
    fft_base(f, 1);
    return f;
}

vector<complex<double>> ifft(vector<complex<double>> f) {
    int n = f.size();
    n = pow(2, ceil(log2(n)));
    f.resize(n, 0);
    fft_base(f, -1);
    for (int i = 0; i < n; i++) f[i] /= n;
    return f;
}

//畳み込み
template <typename RET_T, typename T>
vector<RET_T> convolution(const vector<T> &A, const vector<T> &B) {
    int n = A.size() + B.size() + 1;
    n = pow(2, ceil(log2(n)));
    vector<complex<double>> f(n, 0), g(n, 0);
    for (int i = 0; i < (int)A.size(); i++) f[i] = A[i];
    for (int i = 0; i < (int)B.size(); i++) g[i] = B[i];
    f = fft(f);
    g = fft(g);
    vector<complex<double>> h(n);
    for (int i = 0; i < n; i++) h[i] = f[i] * g[i];
    h = ifft(h);
    vector<RET_T> ret(n);
    for (int i = 0; i < n; i++) ret[i] = h[i].real();
    return ret;
}

bool solve() {
    int L, M, N; cin >> L >> M >> N;
    vector<int> a(N + 1), b(N + 1);
    rep(i, L) {
        int A; cin >> A;
        a[A] = 1;
    }
    rep(i, M) {
        int B; cin >> B;
        b[N - B] = 1;
    }
    vector<double> c = convolution<double, int>(a, b);

    int Q; cin >> Q;
    rep(i, Q) {
        cout << (int)round(c[N + i]) << endl;
    }
    
    return false;
}

int main() {
    cin.tie(0);
    ios::sync_with_stdio(false);
    while(solve());
    return 0;
}
0