bzoj 4836: [Lydsy2017年4月月赛]二元运算 -- 分治+FFT
4836: [Lydsy2017年4月月赛]二元运算
Time Limit: 8 Sec Memory Limit: 128 MB
Description
Input
Output
对于每次查询,输出一行,包含一个整数,表示满足条件的 (i, j) 对的个数。
Sample Input
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
0
1
0
0
1
0
1
0
1
HINT
Source
可能我的常数真的巨大。。。
转载一个姜大爷的题解http://blog.csdn.net/neither_nor/article/details/70889246
#include<map>
#include<cmath>
#include<queue>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
#define inf 1000000007
#define ll long long
#define PI acos(-1)
#define M 400010
inline int rd()
{
int x=,f=;char ch=getchar();
while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}
while(ch>=''&&ch<=''){x=x*+ch-'';ch=getchar();}
return x*f;
}
#define db double
struct cl{
db x,y;
cl(){}
cl(db _x,db _y){x=_x;y=_y;}
friend cl operator +(const cl &x,const cl &y){return cl(x.x+y.x,x.y+y.y);}
friend cl operator -(const cl &x,const cl &y){return cl(x.x-y.x,x.y-y.y);}
friend cl operator *(const cl &x,const cl &y){return cl(x.x*y.x-x.y*y.y,x.x*y.y+x.y*y.x);}
friend cl operator /(const cl &x,const db &y){return cl(x.x/y,x.y/y);}
cl con(){return cl(x,-y);}
};
cl a[M],b[M];
int n,L,R[M];
int A[M],B[M],mx,T;
ll ans[M];
void fft(cl *x,int f,int n)
{
for(int i=;i<n;i++) if(i<R[i]) swap(x[i],x[R[i]]);
for(int i=;i<n;i<<=)
{
cl wn(cos(PI/i),f*sin(PI/i));
for(int j=;j<n;j+=(i<<))
{
cl w(,),X,Y;
for(int k=;k<i;k++,w=w*wn)
{
X=x[j+k];Y=w*x[j+k+i];
x[j+k]=X+Y;x[j+k+i]=X-Y;
}
}
}
if(f==-) for(int i=;i<n;++i) x[i]=x[i]/n;
}
void cal(int l,int r,int n,int L)
{
if(l==r)
{
ans[]+=(ll)A[l]*B[l];
return ;
}
register int i,j,mid=l+r>>;
for(i=;i<n;++i) R[i]=(R[i>>]>>)|((i&)<<L);
for(i=;i<n;++i) a[i]=b[i]=cl(,);
for(i=l;i<=mid;++i) a[i-l].x=A[i];
for(i=mid+;i<=r;++i) a[i-mid-].y=B[i];
fft(a,,n);
for(i=;i<n;i++)
{
j=(n-i)&(n-);
b[i]=(a[i]*a[i]-(a[j]*a[j]).con())*cl(,-0.25);
}
fft(b,-,n);
for(i=;i<n;++i) ans[l+mid++i]+=(ll)(b[i].x+0.1); for(i=;i<n;++i) a[i]=b[i]=cl(,);
for(i=mid+;i<=r;i++) a[i-mid-].x=A[i];
for(i=l;i<=mid;i++) a[mid-i].y=B[i];
fft(a,,n);
for(i=;i<n;i++)
{
j=(n-i)&(n-);
b[i]=(a[i]*a[i]-(a[j]*a[j]).con())*cl(,-0.25);
}
fft(b,-,n);
for(i=;i<n;++i) ans[+i]+=(ll)(b[i].x+0.1); cal(l,mid,n>>,L-);
cal(mid+,r,n>>,L-);
}
int main()
{
register int n1,n2,q,i,x;
T=rd();
while(T--)
{
n1=rd();n2=rd();q=rd();mx=;
memset(A,,sizeof(A));
memset(B,,sizeof(B));
memset(ans,,sizeof(ans));
for(i=;i^n1;++i)
{
x=rd();A[x]++;
mx=max(mx,x);
}
for(i=;i^n2;i++)
{
x=rd();B[x]++;
mx=max(mx,x);
}
for(n=,L=-;n<(mx<<);n<<=) L++;
cal(,n-,n,L);
while(q--)
{
x=rd();
printf("%lld\n",ans[x]);
}
}
return ;
}
bzoj 4836: [Lydsy2017年4月月赛]二元运算 -- 分治+FFT的更多相关文章
- 【bzoj4836】[Lydsy2017年4月月赛]二元运算 分治+FFT
题目描述 定义二元运算 opt 满足 现在给定一个长为 n 的数列 a 和一个长为 m 的数列 b ,接下来有 q 次询问.每次询问给定一个数字 c 你需要求出有多少对 (i, j) 使得 a_ ...
- BZOJ 4881: [Lydsy2017年5月月赛]线段游戏
4881: [Lydsy2017年5月月赛]线段游戏 Time Limit: 3 Sec Memory Limit: 256 MBSubmit: 164 Solved: 81[Submit][St ...
- bzoj 4836 [Lydsy1704月赛]二元运算 分治FFT+生成函数
[Lydsy1704月赛]二元运算 Time Limit: 8 Sec Memory Limit: 128 MBSubmit: 577 Solved: 201[Submit][Status][Di ...
- bzoj4836 [Lydsy2017年4月月赛]二元运算
Description 定义二元运算 opt 满足 现在给定一个长为 n 的数列 a 和一个长为 m 的数列 b ,接下来有 q 次询问.每次询问给定一个数字 c 你需要求出有多少对 (i, j) ...
- bzoj 4885: [Lydsy2017年5月月赛]长方体
Description 给定一个a*b*c的长方体,定义其表面上两个点的距离为沿着长方体的表面走的最短路径的长度,请找到距离最远的点对,你需要保证找到的两个点里至少有一个是长方体顶点. Input 第 ...
- BZOJ 4883 [Lydsy2017年5月月赛]棋盘上的守卫(最小生成环套树森林)
[题目链接] http://www.lydsy.com/JudgeOnline/problem.php?id=4883 [题目大意] 在一个n*m的棋盘上要放置若干个守卫. 对于n行来说,每行必须恰好 ...
- BZOJ 4884 [Lydsy2017年5月月赛]太空猫(单调DP)
[题目链接] http://www.lydsy.com/JudgeOnline/problem.php?id=4884 [题目大意] 太空猫(SpaceCat)是一款画面精致.玩法有趣的休闲游戏, 你 ...
- BZOJ 4836: [Lydsy1704月赛]二元运算 分治FFT
Code: #include<bits/stdc++.h> #define ll long long #define maxn 500000 #define setIO(s) freope ...
- 【BZOJ 4832 】 4832: [Lydsy2017年4月月赛]抵制克苏恩 (期望DP)
4832: [Lydsy2017年4月月赛]抵制克苏恩 Time Limit: 1 Sec Memory Limit: 128 MBSubmit: 275 Solved: 87 Descripti ...
随机推荐
- 微信小程序迁移到头条小程序工具
最近公司需要将微信小程序迁移到头条小程序,比较得知微信和头条小程序的写法类似,只有文件名称不同,相关的指令不同,以及头条在ttml绑定的数据不可使用function,于是就写了node脚本来实现这些重 ...
- Django之组合搜索组件(一)
什么是组合搜索呢? 比如你想买车,但手里只有10万块!所以你只能在10万块的车里挑选,但你喜欢黑色,因为觉得很高端大气上档次,说白了就是装逼杠杠的!之后售车姐给你拿了个表表,你看到了低于10万块且颜色 ...
- 在Oracle中查询表的大小
SELECT segment_name AS TABLENAME,round(BYTES/1024/1024,2) FROM user_segments WHERE segment_name='表名 ...
- 42、和为S的两个数字
一.题目 输入一个递增排序的数组和一个数字S,在数组中查找两个数,是的他们的和正好是S,如果有多对数字的和等于S,输出两个数的乘积最小的. 二.解法 import java.util.ArrayLis ...
- Linux内核中的Cache段
Linux内核中的Cache段 原文地址:http://blogold.chinaunix.net/u2/85263/showart_1743693.html 最近移植LEON3的内核时,了解了一些简 ...
- MySQL之——如何添加新数据库到MySQL主从复制列表 【转】
转自 转载请注明出处:http://blog.csdn.net/l1028386804/article/details/54653691 MySQL主从复制一般情况下我们会设置需要同步的数据库,使用参 ...
- 设计模式之笔记--职责链模式(Chain of Responsibility)
职责链模式(Chain of Responsibility) 定义 职责链模式(Chain of Responsibility),使多个对象都有机会处理请求,从而避免请求的发送者和接收者之间的耦合关系 ...
- 003iptables 命令介绍
http://www.cnblogs.com/wangkangluo1/archive/2012/04/19/2457072.html iptables 防火墙可以用于创建过滤(filter)与NAT ...
- int(long) 类型转换为char
char类型占一个字节,8位 int类型四个字节32位 (long类型的转换跟int类型相同) #include <stdio.h> ]) { buffer[] = (char)tmp; ...
- Photon3Unity3D.dll 解析四——LitePeer
LitePeer 玩家 Connect 连接服务器 Disconnect 断开与服务器的连接 OpJoin 进入游戏 OpLeave 离开游戏,但仍与服务器保持连接 ...