BZOJ4836[Lydsy2017年4月月赛]二元运算 分治+FFT
Posted CQzhangyu
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了BZOJ4836[Lydsy2017年4月月赛]二元运算 分治+FFT相关的知识,希望对你有一定的参考价值。
【BZOJ4836】[Lydsy2017年4月月赛]二元运算
Description
定义二元运算 opt 满足
现在给定一个长为 n 的数列 a 和一个长为 m 的数列 b ,接下来有 q 次询问。每次询问给定一个数字 c
你需要求出有多少对 (i, j) 使得 a_i opt b_j=c 。
Input
第一行是一个整数 T (1≤T≤10) ,表示测试数据的组数。
对于每组测试数据:
第一行是三个整数 n,m,q (1≤n,m,q≤50000) 。
第二行是 n 个整数,表示 a_1,a_2,?,a_n (0≤a_1,a_2,?,a_n≤50000) 。
第三行是 m 个整数,表示 b_1,b_2,?,b_m (0≤b_1,b_2,?,b_m≤50000) 。
第四行是 q 个整数,第 i 个整数 c_i (0≤c_i≤100000) 表示第 i 次查询的数。
Output
对于每次查询,输出一行,包含一个整数,表示满足条件的 (i, j) 对的个数。
Sample Input
2
2 1 5
1 3
2
1 2 3 4 5
2 2 5
1 3
2 4
1 2 3 4 5
2 1 5
1 3
2
1 2 3 4 5
2 2 5
1 3
2 4
1 2 3 4 5
Sample Output
1
0
1
0
0
1
0
1
0 1
0
1
0
0
1
0
1
0 1
题解:本题应该算是分治fft的比较好的入门题了。
按数的大小分治,将a中>mid的数与b中<=mid的数进行fft,将a中<=mid的数与b中>mid的数进行fft,然后递归处理下去即可。
#include <cstdio> #include <iostream> #include <cstring> #include <cmath> #define pi acos(-1.0) using namespace std; const int maxn=200010; typedef long long ll; struct cp { double x,y; cp () {} cp (double a,double b){x=a,y=b;} cp operator + (const cp &a) const {return cp(x+a.x,y+a.y);} cp operator - (const cp &a) const {return cp(x-a.x,y-a.y);} cp operator * (const cp &a) const {return cp(x*a.x-y*a.y,x*a.y+y*a.x);} }C[maxn],D[maxn]; int A[maxn],B[maxn]; ll ans[maxn]; int n,m,q,mx; void FFT(cp *a,int len,int f) { int i,j,k,h; cp t; for(i=k=0;i<len;i++) { if(i>k) swap(a[i],a[k]); for(j=len>>1;(k^=j)<j;j>>=1); } for(h=2;h<=len;h<<=1) { cp wn(cos(2*pi*f/h),sin(2*pi*f/h)); for(j=0;j<len;j+=h) { cp w(1,0); for(k=j;k<j+h/2;k++) t=a[k+h/2]*w,a[k+h/2]=a[k]-t,a[k]=a[k]+t,w=w*wn; } } } void solve(int l,int r) { if(l==r) { ans[0]+=(ll)A[l]*B[l]; return ; } int mid=(l+r)>>1,len,i; for(len=1;len<(r-l+1);len<<=1); memset(C,0,sizeof(C[0])*len),memset(D,0,sizeof(D[0])*len); for(i=l;i<=mid;i++) C[i-l].x=A[i]; for(i=mid+1;i<=r;i++) D[i-mid-1].x=B[i]; FFT(C,len,1),FFT(D,len,1); for(i=0;i<len;i++) C[i]=C[i]*D[i]; FFT(C,len,-1); for(i=0;i<len;i++) ans[i+l+mid+1]+=ll(C[i].x/len+0.1); memset(C,0,sizeof(C[0])*len),memset(D,0,sizeof(D[0])*len); for(i=mid+1;i<=r;i++) C[i-mid-1].x=A[i]; for(i=l;i<=mid;i++) D[mid-i].x=B[i]; FFT(C,len,1),FFT(D,len,1); for(i=0;i<len;i++) C[i]=C[i]*D[i]; FFT(C,len,-1); for(i=0;i<len;i++) ans[i+1]+=ll(C[i].x/len+0.1); solve(l,mid),solve(mid+1,r); } inline int rd() { int ret=0,f=1; char gc=getchar(); while(gc<‘0‘||gc>‘9‘) {if(gc==‘-‘)f=-f; gc=getchar();} while(gc>=‘0‘&&gc<=‘9‘) ret=ret*10+gc-‘0‘,gc=getchar(); return ret*f; } void work() { n=rd(),m=rd(),q=rd(); int i,a; mx=0,memset(A,0,sizeof(A)),memset(B,0,sizeof(B)),memset(ans,0,sizeof(ans)); for(i=1;i<=n;i++) a=rd(),A[a]++,mx=max(mx,a); for(i=1;i<=m;i++) a=rd(),B[a]++,mx=max(mx,a); solve(0,mx); for(i=1;i<=q;i++) a=rd(),printf("%lld\n",ans[a]); } int main() { int T=rd(); while(T--) work(); return 0; }//2 2 1 5 1 3 2 1 2 3 4 5 2 2 5 1 3 2 4 1 2 3 4 5
以上是关于BZOJ4836[Lydsy2017年4月月赛]二元运算 分治+FFT的主要内容,如果未能解决你的问题,请参考以下文章
bzoj4836[Lydsy2017年4月月赛]二元运算 分治+FFT
[Bzoj4832][Lydsy2017年4月月赛]抵制克苏恩 (期望dp)
BZOJ 4832 4832: [Lydsy2017年4月月赛]抵制克苏恩 (期望DP)
BZOJ4832[Lydsy2017年4月月赛]抵制克苏恩 概率与期望