https://ac.nowcoder.com/acm/contest/917/H
链接:https://ac.nowcoder.com/acm/contest/917/H
来源:牛客网
题目背景
把一张纸对折100次就和珠穆朗玛峰一样高了哦
——syh
题目描述
水宝宝某天突发奇想,在纸上写了这样一个问题:
给一个a数组,求∑∞i=1get(l,r,i)⋅∑∞i=1get(l1,r1,i)∑i=1∞get(l,r,i)⋅∑i=1∞get(l1,r1,i)
,其中get(l,r,x)表示求a[l]~a[r]中x出现了几次,他很快推出了规律,
但正当他把这道题录入电脑是发现作为一个蒟蒻的他不会打latex也没找到数学符号(主要是懒),
所以他省略了那个∑∑式子,于是,题面变为了求get(l,r,x)*get(l1,r1,x)的值了.
为了防数据过大,你要对20180623取模
不保证l<r,l1<r1,遇到这种情况,请先交换一下
注:本系列题不按难度排序哦
输入描述:
第一行一个n,m 接下来一行n个数表示a[i] 接下来m行,每行l,r,l1,r1,x,表示求get(l,r,x)*get(l1,r1,x)
输出描述:
3×m行,先输出get(l,r,x),再输出get(l1,r1,x),再输出get(l,r,x)*get(l1,r1,x)
示例1
输入
复制
5 1 2 2 2 2 2 1 5 1 3 2
输出
复制
5 3 15
说明
对于所有数据1<=n,m<=105;1<=l,r,l1,r1<=1051<=n,m<=105;1<=l,r,l1,r1<=105;a[i]≤10000,x≤100000;a[i]≤10000,x≤100000 提示:20180623不是质数
看了别人AC的代码,暴力过的!!!
我??????这n^2的复杂度就过了???搞不懂,实在搞不懂。
自己第一眼看去就就是莫队算法,奈何没注意到区间还能不规范,一直找了很久的错误,于是一两个小时就过去了。哎。。
这已经不知道是第几个低级错误了。
之所以发出来,也是因为这是自己第一次动手写莫队算法。
这题没什么说的,就是莫队板子:
#include<bits/stdc++.h>
#define forn(i,a,b) for(ll i=a;i<=b;i++)
#define INF 0x3f3f3f3f
using namespace std;
#define ll long long
const ll N =100*1000+10;
const ll mod = 20180623;
ll n,m;
ll a[N],ans[N*2],num[N],pos[N];
struct node{
ll L,R,x,id;
}t[N*2];
bool cmp(node a,node b){
if(pos[a.L]==pos[b.L])
return a.R<b.R;
else return pos[a.L]<pos[b.L];
}
ll L=1,R=0;
void solve(node p){
while(p.L>L){
num[a[L]]--;
L++;
}
while(p.R>R){
R++;
num[a[R]]++;
}
while(p.L<L){
L--;
num[a[L]]++;
}
while(p.R<R){
num[a[R]]--;
R--;
}
}
int main(){
scanf("%lld %lld",&n,&m);
memset(num,0,sizeof(num));
ll block=sqrt(n);
for(ll i=1;i<=n;i++){
scanf("%lld",a+i);
pos[i]=i/block;
}
for(ll i=1;i<=m;i++){
ll l1,r1,l2,r2,x;
scanf("%lld%lld%lld%lld%lld",&l1,&r1,&l2,&r2,&x);
if(l1>r1)swap(l1,r1);
if(l2>r2)swap(l2,r2);
t[i].L=l1,t[i].R=r1,t[i].id=i,t[i].x=x;
t[i+m].L=l2,t[i+m].R=r2,t[i+m].id=i+m,t[i+m].x=x;
}
sort(t+1,t+1+2*m,cmp);
memset(num,0,sizeof(num));
for(ll i=1;i<=2*m;i++){
solve(t[i]);
ans[t[i].id]=num[t[i].x];
}
for(ll i=1;i<=m;i++)
printf("%lld\n%lld\n%lld\n",ans[i]%mod,ans[i+m]%mod,(ans[i]%mod*ans[i+m])%mod);
return 0;
}