正题

题目链接:https://www.luogu.com.cn/problem/U138580

题目大意

两个长度为nnn的序列，要求重排后同位置的相加后众数的个数最多。

解题思路

定义aia_iai为第一个序列中iii的出现次数，那么同理有bbb
不难对于每个众数iii的出现次数cic_ici，那么有cn=∑i=1n−1min{ai,bn−i}c_n=\sum_{i=1}^{n-1}min\{a_i,b_{n-i}\}cn=i=1∑n−1min{ai,bn−i}

发现这个东西很难处理，但是如果每个ai,bia_i,b_iai,bi都是111那么同化为卷积可以用NTTNTTNTT进行计算了。

所以我们可以对于每个zzz分开计算也就是cn=∑i=1n−1(ai≥z)∗(bn−i≥z)c_n=\sum_{i=1}^{n-1}(a_i\geq z)*(b_{n-i}\geq z)cn=i=1∑n−1(ai≥z)∗(bn−i≥z)

但是这样zzz的个数非常之多，但是当zzz很大时这个序列会变的很稀疏，那么我们定义一个值TTT，当z>Tz>Tz>T时我们就用暴力处理即可。

codecodecode

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
#define ll long long
using namespace std;
const ll N=3e5+10,T=10,XJQ=998244353;
ll n,m,lg,invn,maxs,ans[N],a[N],b[N];
ll g[N],gi[N],r[N],x[N],y[N];
vector<ll > pa,pb;
ll power(ll x,ll b){ll ans=1;while(b){if(b&1)ans=ans*x%XJQ;x=x*x%XJQ;b>>=1;}return ans;
}
void NTT(ll *f,ll op){for(ll i=1;i<n;i++)if(r[i]<i)swap(f[i],f[r[i]]);for(ll p=2;p<=n;p<<=1){ll len=p>>1,tmp=g[p];if(op==-1)tmp=gi[p];for(ll k=0;k<n;k+=p){ll buf=1;for(ll i=k;i<k+len;i++){ll tt=buf*f[i+len]%XJQ;f[i+len]=(f[i]-tt+XJQ)%XJQ;f[i]=(f[i]+tt)%XJQ;buf=buf*tmp%XJQ;}}}if(op==-1)for(ll i=1;i<=n;i++)f[i]=f[i]*invn%XJQ;return;
}
void solve(ll k){for(ll i=0;i<n;i++)x[i]=(a[i]>=k);for(ll i=0;i<n;i++)y[i]=(b[i]>=k);NTT(x,1);NTT(y,1);for(ll i=0;i<n;i++)x[i]=x[i]*y[i]%XJQ;NTT(x,-1);for(ll i=0;i<n;i++)ans[i]+=x[i];return;
}
int main()
{scanf("%lld",&m);for(n=1;n<=2e5;n<<=1)lg++;invn=power(n,XJQ-2);for(ll i=2;i<=n;i<<=1)g[i]=power(3,(XJQ-1)/i),gi[i]=power(g[i],XJQ-2);for(ll i=0;i<n;i++)r[i]=(r[i>>1]>>1)|((i&1)<<(lg-1));for(ll i=1;i<=m;i++){ll x;scanf("%lld",&x);a[x]++;}for(ll i=1;i<=m;i++){ll x;scanf("%lld",&x);b[x]++;}for(ll i=1;i<=T;i++)solve(i);for(ll i=1;i<=1e5;i++)if(a[i]>T)pa.push_back(i);for(ll i=1;i<=1e5;i++)if(b[i]>T)pb.push_back(i);for(ll i=0;i<pa.size();i++)for(ll j=0;j<pb.size();j++)ans[pa[i]+pb[j]]+=min(a[pa[i]],b[pb[j]])-T;for(ll i=0;i<=n;i++) maxs=max(maxs,ans[i]);printf("%lld",maxs);return 0;
}

[2020.11.4NOIP模拟赛]简单的打击【NTT】相关推荐

[2020.11.27NOIP模拟赛]中位数之中位数【二分,树状数组】
正题题目链接:https://www.luogu.com.cn/problem/U142585?contestId=37855 题目大意 nnn个数,求所有区间中位数的中位数. 解题思路二分一个答 ...
[2020.11.26NOIP模拟赛]勇者的后缀【SA,RMQ,主席树,二分】
正题题目链接:https://www.luogu.com.cn/problem/U142356?contestId=37784 题目大意一个字符串,询问给出(x,l,r)(x,l,r)(x,l,r ...
[2020.11.25NOIP模拟赛]出租车【dp】
正题题面链接:https://www.luogu.com.cn/problem/U142298?contestId=37766 题目大意 nnn个人有起点和终点,按顺序上车,但下车可以任意顺序,车最 ...
[2020.11.27NOIP模拟赛]拼图王【dp】
正题题面链接:https://www.luogu.com.cn/problem/U142584 题目大意 nnn个010101串,按顺序分成两个序列,然后拼接成一个序列(拼接串x,yx,yx,y的话 ...
[2020.11.26NOIP模拟赛]询问【字符串hash】
正题题目链接:https://www.luogu.com.cn/problem/U142342?contestId=37784 题目大意一个字符串,定义两个字符串相似为用一些字母代替相同的字母后可 ...
[2020.11.25NOIP模拟赛]下棋【dp】
正题题目链接:https://www.luogu.com.cn/problem/U142297?contestId=37766 题目大意 nnn个白棋mmm个黑棋排成一排,要求没有任何一段黑白棋的个 ...
[2020.11.3NOIP模拟赛]选数字【容斥】
正题题目链接:https://www.luogu.com.cn/problem/U138404?contestId=36493 题目大意 nnn个数字,每次询问一个区间有多少个三对数或为xxx. 解 ...
2017.6.11 校内模拟赛
题面及数据及std(有本人的也有原来的) :2017.6.11 校内模拟赛 T1 自己在纸上模拟一下后就会发现可以用栈来搞一搞事情受了上次zsq 讲的双栈排序的启发.. 具体就是将原盘子大小cop ...
9.11 myl模拟赛
9.11 myl 模拟赛 100 + 100 + 0 第一题耗费了太多的时间,导致最后一题没有时间想,直接去写了暴力,而且出题人没有给暴力分.... Problem 1. superman [题目描述 ...

[2020.11.4NOIP模拟赛]简单的打击【NTT】

正题

题目大意

解题思路

codecodecode

[2020.11.4NOIP模拟赛]简单的打击【NTT】相关推荐

最新文章

热门文章