题目

给出两个等长的序列\(a,b\),

重排序列\(b\),使得\(a+b\)众数出现的次数最多


分析

设\(f[i]\)表示众数为\(i\)的贡献,那么

\(f[i]=\sum_{j<i}min(A[j],B[i-j])\)

其中大写的\(a,b\)表示次数,但是这个东西很难做,

考虑把它变成正常的卷积的形式,那么就转换成判定,

判定\(f[i]\)是否能够达到阈值,那么显然就变成了只有\(0,1\)的卷积,

但是这样会超时,考虑阈值仅限于不超过一个常数就可以了


代码

#include <cstdio>
#include <cctype>
#include <queue>
#define rr register
using namespace std;
const int mod=998244353,N=100011,inv3=332748118; priority_queue<pair<int,int> >q;
int a[N],b[N],A[N],B[N],ff[N<<2],n,TOT[2],lim,gg[N<<2],ans[N<<2],Ans,Gmi[31],Imi[31];
inline signed iut(){
rr int ans=0; rr char c=getchar();
while (!isdigit(c)) c=getchar();
while (isdigit(c)) ans=(ans<<3)+(ans<<1)+(c^48),c=getchar();
return ans;
}
inline signed min(int a,int b){return a<b?a:b;}
inline signed mo1(int x,int y){return x+y>=mod?x+y-mod:x+y;}
inline signed mo2(int x,int y){return x<y?x-y+mod:x-y;}
inline signed ksm(int x,int y){
rr int ans=1;
for (;y;y>>=1,x=1ll*x*x%mod)
if (y&1) ans=1ll*ans*x%mod;
return ans;
}
namespace Theoretic{
int rev[N<<2],LAST,tt[N<<2];
inline void Pro(int n){
if (LAST==n) return; LAST=n;
for (rr int i=0;i<n;++i)
rev[i]=(rev[i>>1]>>1)|((i&1)?n>>1:0);
}
inline void NTT(int *f,int n,int op){
Pro(n);
for (rr int i=0;i<n;++i)
if (i<rev[i]) swap(f[i],f[rev[i]]);
rr int p=2,len=1;
for (rr int o=1;p<=n;++o){
rr int W=(op==1)?Gmi[o]:Imi[o];
for (rr int i=0;i<n;i+=p){
rr int t=1;
for (rr int j=i;j<i+len;++j){
rr int z=1ll*f[len+j]*t%mod;
f[len+j]=mo2(f[j],z),f[j]=mo1(f[j],z);
t=1ll*t*W%mod;
}
}
p<<=1,len<<=1;
}
}
inline void Cb(int *f,int *g,int n){
for (rr int i=0;i<n;++i) f[i]=1ll*f[i]*g[i]%mod;
}
inline void times(int *f,int *g,int len,int lim){
rr int n=1,invn; for (;n<lim;n<<=1);
for (rr int i=0;i<len;++i) tt[i]=g[i];
for (rr int i=len;i<n;++i) tt[i]=0;
NTT(f,n,1),NTT(tt,n,1),Cb(f,tt,n),NTT(f,n,-1);
for (rr int i=lim;i<n;++i) f[i]=0;
for (rr int i=0;i<n;++i) tt[i]=0;
invn=ksm(n,mod-2);
for (rr int i=0;i<lim;++i)
f[i]=1ll*f[i]*invn%mod;
}
}
inline void GmiImi(){
for (rr int i=0;i<31;++i) Gmi[i]=ksm(3,(mod-1)/(1<<i));
for (rr int i=0;i<31;++i) Imi[i]=ksm(inv3,(mod-1)/(1<<i));
}
inline void Prep(int lim){
for (rr int i=1;i<=n;++i) ff[i]=a[i]>=lim;
for (rr int i=1;i<=n;++i) gg[i]=b[i]>=lim;
Theoretic::times(ff,gg,n,n*2);
for (rr int i=1;i<=n;++i) ans[i]+=ff[i];
}
signed main(){
n=iut(),GmiImi();
for (rr int i=1;i<=n;++i) ++a[iut()];
for (rr int i=1;i<=n;++i) ++b[iut()];
for (rr int i=1;i<N;++i) if (a[i]) q.push(make_pair(a[i],0));
for (rr int i=1;i<N;++i) if (b[i]) q.push(make_pair(b[i],1));
while (!q.empty()){
rr int t=q.top().second; q.pop();
if ((++TOT[t])*TOT[t^1]>1e8) break;
}
if (!q.empty()) lim=q.top().first+1;
for (rr int i=1;i<=lim;++i) Prep(i);
for (rr int i=1;i<=n;++i) if (a[i]>lim) A[++A[0]]=i;
for (rr int i=1;i<=n;++i) if (b[i]>lim) B[++B[0]]=i;
for (rr int i=1;i<=A[0];++i)
for (rr int j=1;j<=B[0];++j)
ans[A[i]+B[j]]+=min(a[A[i]],b[B[j]])-lim;
for (rr int i=1;i<=2*n;++i)
if (Ans<ans[i]) Ans=ans[i];
return !printf("%d",Ans);
}

#NTT,DP#U138580 简单的打击的更多相关文章

  1. ZOJ 4257 MostPowerful(状压DP,简单)

    题目大意:不超过10种气体,两两之间相互碰撞可以产生一定的能量,如a碰b,那么b气体就消失,自身不能碰自身,问最后所能得到的最大能量. 原代码链接:http://blog.csdn.net/accry ...

  2. 数位dp 的简单入门

    时间紧张,就不讲那么详细了. 之前一直被深搜代码误解,以为数位dp 其实就是记忆化深搜...(虽说爆搜确实很舒服而且还好想) 但是后来发现数位dp 的标准格式其实是 预处理 + dp ...... 数 ...

  3. 斜率优化dp 的简单入门

    不想写什么详细的讲解了...而且也觉得自己很难写过某大佬(大米饼),于是建议把他的 blog 先看一遍,然后自己加了几道题目以及解析...顺便建议看看算法竞赛(蓝皮书)的 0x5A 斜率优化(P294 ...

  4. PKU 1458 Common Subsequence(最长公共子序列,dp,简单)

    题目 同:ZJU 1733,HDU 1159 #include <stdio.h> #include <string.h> #include <algorithm> ...

  5. DP的简单应用

    Problem A:简单的图形覆盖 Time Limit:1000MS  Memory Limit:65536KTotal Submit:201 Accepted:104 Description 有一 ...

  6. dp优化简单总结

    1.二分优化 (使用二分查找优化查找效率) 典型例题:LIS dp[i]保存长度为 i 的上升子序列中最小的结尾,可以用二分查找优化到nlogn 2.数学优化 (通过数学结论减少状态数) 例题1:hd ...

  7. HDU 1024 Max Sum Plus Plus(DP的简单优化)

    Problem Description Now I think you have got an AC in Ignatius.L's "Max Sum" problem. To b ...

  8. dp的简单递推笔记1

    (1)转自rockZ的博文 UVa 10328 - Coin Toss (递推) 题意:给你一个硬币,抛掷n次,问出现连续至少k个正面向上的情况有多少种. 原题中问出现连续至少k个H的情况,很难下手. ...

  9. 2017-5-14 湘潭市赛 Similar Subsequence 分析+四维dp+一些简单优化

    Similar Subsequence Accepted : Submit : Time Limit : MS Memory Limit : KB Similar Subsequence For gi ...

  10. 树形DP(简单题)(Y HDU4705)

    题意:给出一个n个节点的树形图,统计{A,B,C}的数量,其中ABC分别是树上三个不同的节点,并且这三个节点不能被一条路径覆盖 分析:对于下图 进行dfs深搜统计,num[u]统计回溯到当前节点u,并 ...

随机推荐

  1. zynq7000 I2C RTC 与 串口使用

    RS485 串口 测试 硬件上2路串口,其中UART 1对应PS STD IN/OUT,UART 0对应RS485: 图 ‑1 RS485电路,自动转换输入.输出方向 可参考 https://blog ...

  2. Spring + JAX-WS : ‘xxx’ is an interface, and JAXB can’t handle interfaces 错误解决方法

    错误栈 Caused by: com.sun.xml.bind.v2.runtime.IllegalAnnotationsException: 2 counts of IllegalAnnotatio ...

  3. Apache Hudi 在 vivo 湖仓一体的落地实践

    作者:vivo 互联网大数据团队 - Xu Yu 在增效降本的大背景下,vivo大数据基础团队引入Hudi组件为公司业务部门湖仓加速的场景进行赋能.主要应用在流批同源.实时链路优化及宽表拼接等业务场景 ...

  4. Java 常用类 于 StringBuffer 和 StringBuilder的使用 + String三者的异同

    1 package com.bytezero.stringclass; 2 3 import org.junit.Test; 4 5 /** 6 * 关于 StringBuffer 和 StringB ...

  5. 使用内网nginx代理rancher

    需求: rancher 部署在 192.168.188.167服务器上 控制台访问地址 https://192.168.188.167:8443 在本地只有192.168.80.111有权限访问, 需 ...

  6. SPFA最短路

    目录 从Bellman-Ford开始 核心思想 模拟算法执行过程 时间复杂度 模板 spfa spfa优化的思想 模板 从Bellman-Ford开始 对于所有边权都大于等于0的图,任意两个顶点之间的 ...

  7. 小工具 --- 百度翻译API翻译工具

    引言 最近想把一些英文官方文档的资料翻译成中文,然后转化为Markdown文档,然后发现百度通用翻译的API有不错的免费额度,个人申请也能申请到高级版.这个额度足够个人的日常使用了. 如何使用 如何使 ...

  8. 没有有线网卡的笔记本如何在PVE下All in one?—NAS + Linux +win下载机

    没有有线网卡的笔记本在PVE下All in one | NAS + Linux + Win下载机 (保姆级未完成版) 目录: 1.前言 2.PVE的安装 3.PVE联网前的准备工作 4.PVE使用无线 ...

  9. 35_音视频播放器_seek&暂停

    目录 一.实现seek功能 二.解决点击seek操作时会出现画面快速闪过 三.实现暂停功能 3.1.音频暂停 3.2.视频暂停 一.实现seek功能 我们主要是使用ffmpeg的av_seek_fra ...

  10. PAT 甲级【1015 Reversible Primes】

    考察素数判断 考察进制转换 import java.io.IOException; import java.io.InputStreamReader; import java.io.StreamTok ...