21牛客多校h（fft/ntt）

2023-11-22 12:12:51

题目大意：找一个最小的m，使得n个数对m求模后无哈希冲突

思路：当 |ai - aj | ==k*m时存在哈希冲突，问题转化为找最小m，使得任意两个数的差不为m的倍数

大方向是从小到大枚举m，看它的倍数是否存在即可，枚举的复杂度是调和级数nlongn

判存在的话需要预处理

转化成多项式乘法，若d存在则x^d的系数为1（不为0就行），反之为0

然后fft/ntt算下卷出来的系数，因为这里是ai-aj，我们给后者加个5e5，避免出现负数

#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
using namespace std;
namespace NTT
{
const long long g = 3;
const long long p = 998244353;
long long wn[35];
long long pow2(long long a, long long b)
{long long res = 1;while (b){if (b & 1) res = res * a % p;a = a * a % p;b >>= 1;}return res;
}
void getwn()
{for (int i = 0; i < 25; i++) wn[i] = pow2(g, (p - 1) / (1LL << i));
}
void ntt(long long *a, int len, int f)
{long long i, j = 0, t, k, w, id;for (i = 1; i < len - 1; i++){for (t = len; j ^= t >>= 1, ~j & t;);if (i < j) swap(a[i], a[j]);}for (i = 1, id = 1; i < len; i <<= 1, id++){t = i << 1;for (j = 0; j < len; j += t){for (k = 0, w = 1; k < i; k++, w = w * wn[id] % p){long long x = a[j + k], y = w * a[j + k + i] % p;a[j + k] = (x + y) % p;a[j + k + i] = (x - y + p) % p;}}}if (f){for (i = 1, j = len - 1; i < j; i++, j--) swap(a[i], a[j]);long long inv = pow2(len, p - 2);for (i = 0; i < len; i++) a[i] = a[i] * inv % p;}
}
void mul(long long *a, long long *b, int l1, int l2)
{int len, i;for (len = 1; len <= l1 + l2; len <<= 1);for (i = l1 + 1; i <= len; i++) a[i] = 0;for (i = l2 + 1; i <= len; i++) b[i] = 0;ntt(a, len, 0); ntt(b, len, 0);for (i = 0; i < len; i++) a[i] = a[i] * b[i] % p;ntt(a, len, 1);
}
};
typedef long long ll;
#define _for(i,a,b) for(int i=(a) ;i<=(b) ;i++)
#define _rep(i,a,b) for(int i=(a) ;i>=(b) ;i--)
#define scd(v) scanf("%d",&v)
#define scdd(a,b) scanf("%d %d",&a,&b)
#define endl "\n"
#define IOS ios::sync_with_stdio(false)
#define pb push_back
#define all(v) v.begin(),v.end()
#define int long long
#define odd(x) x&1
#define mst(v,a) memset(v,a,sizeof(v))
#define lson p<<1 ,l,mid
#define rson p<<1|1,mid+1,r
#define ls p<<1
#define rs p<<1|1
#define fi first
#define se second
#define pii pair
#define inf 0x7f7f7f7f
const int N=5e5+10;
//const int mod=998244353;
int n,m,k;
int ki_1[N*3],ki_2[N*3];
void solve();
signed main()
{//!!//
//    freopen("data.txt","r",stdin);//!!IOS;NTT::getwn();cin>>n;_for(i,1,n){int x;cin>>x;ki_1[x]++;//ai系数ki_2[500000-x]++;//-aj系数}NTT::mul(ki_1,ki_2,5e5,5e5);//ntt加速求ki_1与系数ki_2的卷积bool ok=1;for(int i=1;;i++)//枚举m{ok=1;for(int j=i ;j<=5e5 ;j+=i)//枚举m的倍数{if( ki_1[j+500000] ){ok=0;break;}}if( ok ){cout<

 
                        
                        
本文来自互联网用户投稿，文章观点仅代表作者本人，不代表本站立场，不承担相关法律责任。如若转载，请注明出处。 如若内容造成侵权/违法违规/事实不符，请点击【内容举报】进行投诉反馈！



                    



                    

    收藏
    



                    
    
        
        标签：技术
        
    

    
        
                
            上一篇 >
            FastAPI 1 - 基本使用
        
                
            下一篇 >
            【JS协议UI源码】autoJS协议源码，JavaScript界面源码，AJ脚本利用网盘就能远程关软
        
                
    



                    
    
        
        
            相关文章
        
                
            Duilib中list控件支持ctrl和shif多行选中的实现
        
                
            [ICML2015]Batch Normalization:Accelerating Deep Network Training by Reducing Internal Covariate Shif
        
                
            win10系统 微软输入法 于eclipse ctrl+shif+f冲突间接处理办法
        
                
            Codeforces Round #259 (Div. 2) B. Little Pony and Sort by Shif
        
                
            读LDD3，内存映射与DMA--PAGE_SHIF…
        
                
            VMware虚拟机安装XP【要先分区，再设置BOOT 启动CD，shif+上移】
        
                
            更换iBus五笔的左与右Shif
        
                
            sublime ctrl+shif+f 没用解决办法
        
                
            idea 对 ctrl + z 的撤销 是 ctrl + shif + z
        
                
            计算机最早的设计师应用于,计算机应用基础选择题doc.doc
        
                
            win10自带截图神器：Win+Shift+S
        
                
            Python基础之文件目录操作
        
                
            python简述目录_Python基础之文件目录操作(示例代码)
        
                
            tp5 如何做数据采集
        
                
            任务2-7(服务器字体+阿里巴巴矢量库)
        
                
            html标签（1)：h1~h6,p,br,pre,hr
        
                
            TI 电量计介绍与芯片选型指南
        
                
            几款TI电源芯片简介
        
                
            TI DSP芯片C2000系列读取FLASH数据
        
                
            德州仪器(Ti)平台嵌入式开发基础
        
                
            TI三相电机智能栅极驱动芯片特点分类
        
                
            省选模拟（12.08） T3 圈圈圈圈圈圈圈圈
        
                
            Hadoop生态圈技术栈（上）
        
                
            大数据开发基础入门与项目实战（三）Hadoop核心及生态圈技术栈之6.Impala交互式查询
        
                
            小猿圈之Linux下Mysql 操作命令
        
                
            大数据Hadoop生态圈常用面试题
        
                
            大数据开发基础入门与项目实战（三）Hadoop核心及生态圈技术栈之4.Hive DDL、DQL和数据操作
        
                
            备战Noip2018模拟赛11（B组）T3 Monogatari 物语
        
                
            【智能优化算法-圆圈搜索算法】基于圆圈搜索算法Circle Search Algorithm求解单目标优化问题附matlab代码
        
                
            NYOJ 78 圈水池
        
                
            递归问题 跑道 汽车 绕圈问题 Python实现
        
                
            Hadoop生态圈（三）：MapReduce