結果

問題 No.206 数の積集合を求めるクエリ
コンテスト
ユーザー kopricky
提出日時 2017-07-23 07:11:38
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 142 ms / 7,000 ms
コード長 2,828 bytes
コンパイル時間 1,552 ms
コンパイル使用メモリ 165,972 KB
実行使用メモリ 17,292 KB
最終ジャッジ日時 2024-10-09 10:42:52
合計ジャッジ時間 4,662 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 28
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function ‘int main()’:
main.cpp:92:10: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   92 |     scanf("%d%d%d",&l,&m,&n);
      |     ~~~~~^~~~~~~~~~~~~~~~~~~
main.cpp:95:14: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   95 |         scanf("%d",&a[i]);
      |         ~~~~~^~~~~~~~~~~~
main.cpp:98:14: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   98 |         scanf("%d",&b[i]);
      |         ~~~~~^~~~~~~~~~~~
main.cpp:111:10: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
  111 |     scanf("%d",&q);
      |     ~~~~~^~~~~~~~~

ソースコード

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