jzoj3895-数字对【RMQ,GCD,二分答案,单调队列】
生活随笔
收集整理的這篇文章主要介紹了
jzoj3895-数字对【RMQ,GCD,二分答案,单调队列】
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
正題
題目大意
一個序列aaa
對于[L..R][L..R][L..R]若ak∈[L..R]∣ai∈[L..R]a_k \in [L..R]|a_i\in[L..R]ak?∈[L..R]∣ai?∈[L..R]則這個一個特殊區間。
求最長特殊區間。
解題思路
先RMQRMQRMQ求區間GCDGCDGCD,然后二分答案。
之后aka_kak?肯定是這個區間最小的,單調隊列維護一下(或者RMQRMQRMQ維護),然后判斷區間GCDGCDGCD是否等于aka_kak?
codecodecode
#include<cstdio> #include<queue> #include<cmath> #include<cctype> #include<algorithm> #define N 500010 using namespace std; deque<int> q; int n,a[N],ans[N],tot,f[N][23],lg[N]; int read() {int x=0,f=1; char c=getchar();while(!isdigit(c)) {if(c=='-')f=-f;c=getchar();}while(isdigit(c)) x=(x<<1)+(x<<3)+c-48,c=getchar();return x*f; } void print(int x){if (x>9) print(x/10); putchar(x%10+48); return; } void Keep(int x){while(!q.empty()&&q.front()<x)q.pop_front(); } void Push(int x){while(!q.empty()&&a[q.back()]>a[x])q.pop_back();q.push_back(x); } int Gcd(int l,int r){int z=lg[r-l+1];return __gcd(f[l][z],f[r+1-(1<<z)][z]); } int check(int x){int ans=0;while(!q.empty())q.pop_front();for(int i=1;i<=x;i++)Push(i);for(int i=x+1;i<=n;i++){Keep(i-x);Push(i);int w=a[q.front()];if(Gcd(i-x,i)==w)ans++;}return ans; } void write(int x){while(!q.empty())q.pop_front();for(int i=1;i<=x;i++)Push(i);for(int i=x+1;i<=n;i++){Keep(i-x);Push(i);int w=a[q.front()];if(Gcd(i-x,i)==w)print(i-x),putchar(' ');} } int main() {n=read();lg[0]=-1;for(int i=1;i<=n;i++)a[i]=read(),f[i][0]=a[i],lg[i]=lg[i/2]+1;for(int j=1;(1<<j)<=n;j++)for(int i=1;i+(1<<j)-1<=n;i++)f[i][j]=__gcd(f[i][j-1],f[i+(1<<j-1)][j-1]);int l=0,r=n;while(l<=r){int mid=(l+r)/2;if(check(mid)) l=mid+1;else r=mid-1;}print(check(r));putchar(' ');print(r);putchar('\n');write(r); }總結
以上是生活随笔為你收集整理的jzoj3895-数字对【RMQ,GCD,二分答案,单调队列】的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: P3365,jzoj3894-改造二叉树
- 下一篇: jzoj3896-战争游戏【tarjan