#200-[莫队]数列互质

这是一个关于数列中互质次数查询的问题,描述了一个输入n和m,n个正整数a[i]构成的数列,以及m组询问,每组包含li, ri, ki,目标是找出区间[li, ri]内与ki互质的数的个数。题目限制1≤n, m≤5×10^4,1≤a[i]≤n,1≤li≤ri≤n,1≤ki≤n,来源于美团CodeM比赛,属于莫队算法范畴。" 134059562,12425543,CentOS系统中部署Tomcat的详细步骤,"['Linux', 'CentOS', 'Java', 'Tomcat']

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

PS:刷了好久到200篇题解......OIer不容易啊

 

Description

给出一个长度为 n 的数列 a1,a2,a3,...,an,以及 m 组询问 (li,ri,ki),求区间 [li,ri] 中有多少数在该区间中的出现次数与 ki 互质。

Input

第一行,两个正整数 n , m。 

第二行,n 个正整数 ai 描述这个数列。 

接下来 m 行,每行三个正整数 li,ri,ki,描述一次询问。 

Output

输出 m 行,即每次询问的答案。

10 5
1 1 1 1 1 2 2 2 2 2
4 7 2
4 7 3
4 8 2
4 8 3
3 8 3
  • Sample Input

0
2
1
1
0
  • Sample Output

HINT

  • 1≤n,m≤5×104
  • 1≤ai≤n
  • 1≤li≤ri≤n
  • 1≤ki≤n

Source/Category

美团CodeM   

## **题目描述** 给出一个长度为 $n$ 的数列 $\{a_1, a_2, a_3, \dots, a_n\}$,以及 $m$ 组询问 $(l_i, r_i, k_i)$。对于每组询问,求在区间 $[l_i, r_i]$ 中有多少种数字,其在该区间内的出现次数与 $k_i$ 互质(即 $\gcd(\text{出现次数}, k_i) = 1$)。 ## **输入格式** 第一行包含个正整数 $n$ 和 $m$,分别表示数列长度和询问次数。 第二行包含 $n$ 个正整数 $a_i$,描述该数列。 接下来 $m$ 行,每行包含三个正整数 $l_i, r_i, k_i$,描述一次区间询问## **输出格式** 输出 $m$ 行,每行一个整数,表示对应询问的答案。 ## **样例** #### 样例输入 ```plain 10 5 1 1 1 1 1 2 2 2 2 2 4 7 2 4 7 3 4 8 2 4 8 3 3 8 3 ``` #### 样例输出 ```plain 0 2 1 1 0 ``` ## **数据范围与提示** #### **样例解释** - 第一组询问 $[4,7]$:元素为 $[1,1,2,2]$,数字 $1$ 出现 $2$ 次($\gcd(2,2)=2≠1$),数字 $2$ 出现 $2$ 次(不互质),故答案为 $0$。 - 第二组询问 $[4,7]$:数字 $1$ 出现 $2$ 次($\gcd(2,3)=1$),数字 $2$ 出现 $2$ 次($\gcd(2,3)=1$),故答案为 $2$。 #### **数据范围** - $1 \le n, m \le 5 \times 10^4$ - $1 \le a_i \le n$ - $1 \le l_i \le r_i \le n$ - $1 \le k_i \le n$ #include <bits/stdc++.h> #define ll long long using namespace std; const ll N=5e4; struct node{ ll l,r,k,f; }q[N+5]; ll n,m; ll a[N+5]; ll blk[N+5],lenb; bool jo[N+5]; ll pre[N+5],nxt[N+5]; ll head,tail; ll lp=1,rp,s; ll cnt[N+5],f[N+5]; ll ans[N+5]; ll gcd(ll a,ll b){ return b?gcd(b,a%b):a; } bool cmp(node l1,node l2){ if(blk[l1.l]==blk[l2.l]){ if(blk[l1.l]&1) return l1.r<l2.r; return l1.r>l2.r; } return blk[l1.l]<blk[l2.l]; } void addd(ll x){ if(x<1 || x>n || jo[x]) return; jo[x]=1; pre[x]=head; nxt[x]=nxt[head]; pre[nxt[head]]=x; nxt[head]=x; return; } void dell(ll x){ if(x<1 || x>n || !jo[x]) return; jo[x]=1; nxt[pre[x]]=nxt[x]; pre[nxt[x]]=pre[x]; return; } void add(ll x){ f[cnt[a[x]]]--; if(f[cnt[a[x]]]==0) dell(cnt[a[x]]); cnt[a[x]]++; f[cnt[a[x]]]++; if(f[cnt[a[x]]]==1) addd(cnt[a[x]]); return; } void del(ll x){ f[cnt[a[x]]]--; if(f[cnt[a[x]]]==0) dell(cnt[a[x]]); cnt[a[x]]--; f[cnt[a[x]]]++; if(f[cnt[a[x]]]==1) addd(cnt[a[x]]); return; } void solve(ll l,ll r){ while(l<lp){ lp--; add(lp); } while(lp<l){ del(lp); lp++; } while(rp<r){ rp++; add(rp); } while(r<rp){ del(rp); rp--; } return; } int main(){ scanf("%lld%lld",&n,&m); lenb=(ll)sqrt(n); ll sx=1; for(ll i=1;i<=n;i++){ blk[i]=sx; if(i%lenb==0) sx++; } for(ll i=1;i<=n;i++) scanf("%lld",&a[i]); for(ll i=1;i<=m;i++){ scanf("%lld%lld%lld",&q[i].l,&q[i].r,&q[i].k); q[i].f=i; } tail=n+1; nxt[head]=tail; pre[tail]=head; f[0]=n; sort(q+1,q+m+1,cmp); for(ll i=1;i<=m;i++){ solve(q[i].l,q[i].r); for(ll j=nxt[head];j!=tail;j=nxt[j]){ if(gcd(j,q[i].k)==1) ans[q[i].f]+=f[j]; // cout<<j<<" "; } } for(ll i=1;i<=m;i++) printf("%lld\n",ans[i]); return 0; } 过不了样例,修正代码错误并标出
最新发布
08-05
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值