BZOJ4836(CDQ分治+FFT(NTT会T))

博客探讨了一种涉及特殊运算的题目,其中运算包括两种情况。当运算为加法时,问题转化为卷积,可用FFT解决。即使运算为减法,通过数组反转依然可以转换为卷积问题。作者对使用CDQ分治方法解决这个问题表示兴趣,通过二分查找和FFT进行区间卷积,实现O(N (log N)^2)的时间复杂度。然而,注意到在实际运行中,基于NTT的解决方案可能会在在线评测系统中超时。
摘要由CSDN通过智能技术生成

题面
题目定义了一种奇怪的运算,是这样的
这里写图片描述
还给出了两个数组A,B,数组长度和元素大小都<=5*10^4。
有多个询问,给出一个k,问有多少对(i,j)使得
Ai opt Bj = k。

如果 x opt y = x+y,那么显然是一个桶状数组的卷积,FFT可以解决
若 x opt y = x-y,根据FFT题的套路,把一个数组反过来,发现它依然是个卷积。而题目把运算按x与y的大小关系分成了两种,两种都很类似,可以独立地用相同的方法做。

这题是一个不熟FFT的大佬搜CDQ分治搜到的,让我非常感兴趣,所以就来考虑分治。
其实也很简单,不过思路我没听过而已。
对于权值区间[l,r],我们二分一个mid,然后用FFT将[l,mid]和[mid+1,r]做个卷积,翻转一边,再做个卷积,统计答案。正确地清零,保证复杂度。再递归处理两个区间,然后就搞定了。
时间复杂度 O(N (log N)^2)
N=5*10^4。

为什么NTT在电脑上跑得飞快,交上去老是T。

#include <iostream>
#include <fstream>
#include <algorithm>
#include <cmath>
#include <ctime>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <queue>

using namespace std;
#define mmst(a, b) memset(a, b, sizeof(a))
#define mmcp(a, b) memcpy(a, b, sizeof(b))

typedef long long LL;

void read(int &hy)
{
    hy=0;
    char cc=getchar();
    while(cc>'9'||cc<'0')
    cc=getchar();
    while(cc>='0'&&cc<='9')
    {
        hy=(hy<<3)+(hy<<1)+cc-
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值