結果

問題 No.206 数の積集合を求めるクエリ
ユーザー koprickykopricky
提出日時 2017-07-23 07:11:38
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 166 ms / 7,000 ms
コード長 2,828 bytes
コンパイル時間 1,460 ms
コンパイル使用メモリ 152,120 KB
実行使用メモリ 17,260 KB
最終ジャッジ日時 2023-07-30 15:29:52
合計ジャッジ時間 5,101 ms
ジャッジサーバーID
(参考情報)
judge14 / judge11
このコードへのチャレンジ(β)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
4,376 KB
testcase_01 AC 1 ms
4,376 KB
testcase_02 AC 2 ms
4,380 KB
testcase_03 AC 2 ms
4,380 KB
testcase_04 AC 2 ms
4,376 KB
testcase_05 AC 1 ms
4,380 KB
testcase_06 AC 5 ms
4,380 KB
testcase_07 AC 5 ms
4,376 KB
testcase_08 AC 5 ms
4,376 KB
testcase_09 AC 5 ms
4,376 KB
testcase_10 AC 1 ms
4,376 KB
testcase_11 AC 2 ms
4,376 KB
testcase_12 AC 5 ms
4,380 KB
testcase_13 AC 5 ms
4,380 KB
testcase_14 AC 5 ms
4,380 KB
testcase_15 AC 5 ms
4,380 KB
testcase_16 AC 3 ms
4,380 KB
testcase_17 AC 159 ms
17,260 KB
testcase_18 AC 144 ms
16,924 KB
testcase_19 AC 154 ms
17,068 KB
testcase_20 AC 76 ms
10,380 KB
testcase_21 AC 151 ms
16,852 KB
testcase_22 AC 152 ms
16,704 KB
testcase_23 AC 158 ms
17,080 KB
testcase_24 AC 166 ms
17,144 KB
testcase_25 AC 160 ms
17,180 KB
testcase_26 AC 156 ms
16,812 KB
testcase_27 AC 150 ms
16,920 KB
testcase_28 AC 157 ms
16,832 KB
testcase_29 AC 155 ms
17,020 KB
testcase_30 AC 152 ms
16,756 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <bits/stdc++.h>
#define ll long long
#define INF 1000000005
#define MOD 1000000007
#define EPS 1e-10
#define rep(i,n) for(int i=0;i<(int)n;++i)
#define each(a, b) for(auto (a): (b))
#define all(v) (v).begin(),(v).end()
#define fi first
#define se second
#define pb push_back
#define show(x) cout <<#x<<" = "<<(x)<<endl
#define spair(p) cout <<#p<<": "<<p.fi<<" "<<p.se<<endl
#define svec(v) cout<<#v<<":";rep(kbrni,v.size())cout<<" "<<v[kbrni];cout<<endl
#define sset(s) cout<<#s<<":";each(kbrni,s)cout <<" "<<kbrni;cout<<endl

using namespace std;

typedef pair<int,int>P;
typedef complex<double>C;

const int MAX_N = 100005;
const double PI = 4*atan(1.0);

//aにはAiおよびBi(係数)が入っている
vector<C> FFT(double theta, const vector<C>& a){
    int n = (int)a.size();
    vector<C> ret = a;
    for(int m=n; m>=2; m>>=1){
        int mh = m>>1;
        rep(i,mh){
            C w = exp(i*theta*C(0,1));
            for(int j=i; j<n; j+=m){
                int k = j+mh;
                C x = ret[j] - ret[k];
                ret[j] += ret[k];
                ret[k] = w*x;
            }
        }
        theta *= 2;
    }
    int i=0;
    for(int j=1; j<n-1; j++){
        for(int k=n>>1; k>(i^=k); k>>= 1){

        }
        if(j < i) swap(ret[i], ret[j]);
    }
    return ret;
}

//畳み込み
template<class T>
vector<T> Convolution(const vector<T> &lhs, const vector<T> &rhs){
    int n = 1;
    while(n <= (int)(lhs.size() + rhs.size())){
        n <<= 1;
    }
    vector<C> temp1(n);
    vector<C> temp2(n);
    for(int i=0; i<n; i++){
        if(i < (int)lhs.size()){
            temp1[i] = C(lhs[i], 0);
            if(i < (int)rhs.size()){
                temp2[i] = C(rhs[i], 0);
            }
        }else if(i < (int)rhs.size()){
            temp2[i] = C(rhs[i], 0);
        }else{
            break;
        }
    }
    //FFTにかけて周波数分解
    temp1 = FFT(2.0*PI/n, temp1);
    temp2 = FFT(2.0*PI/n, temp2);
    //周波数同士で掛け算
    rep(i,n){
        temp1[i] *= temp2[i];
    }
    //逆FFTをかけて元に戻す
    temp1 = FFT(-2.0*PI/n, temp1);
    vector<T> ret(n);
    rep(i,n){
        ret[i] = (T)(temp1[i].real()/n + 0.5); //T=intの時用
    }
    return ret;
}

int main()
{
    int l,m,n;
    scanf("%d%d%d",&l,&m,&n);
    vector<int> a(l),b(m);
    rep(i,l){
        scanf("%d",&a[i]);
    }
    rep(i,m){
        scanf("%d",&b[i]);
        b[i] = n-b[i];
    }
    vector<int> u(*max_element(all(a))+1,0);
    vector<int> v(*max_element(all(b))+1,0);
    rep(i,l){
        u[a[i]]++;
    }
    rep(i,m){
        v[b[i]]++;
    }
    vector<int> res = Convolution(u,v);
    int q;
    scanf("%d",&q);
    rep(i,q){
        if(n+i < res.size()){
            printf("%d\n",res[n+i]);
        }else{
            printf("0\n");
        }
    }
}
0