GCD
题目大意:给出了n个数,m次查询,每次查询gcd(l,r)并且输出有多少个区间的gcd==gcd(l,r);
解题思路:做的时候误以为要求查询的区间在(l,r)内,并且也没有意识到固定左端点移动右端点gcd是单调不增的。
我们可以固定左端点,然后根据上边那个性质,二分右端点。然后统计贡献。
#include<bits/stdc++.h>
using namespace std;
#define LL long long
const int N = 1e5+5;
map<int,LL>mmp;
int a[N];
int f[N][20];
void build(int n)
{
for(int i=1;i<=n;i++)f[i][0]=a[i];
for(int i=1;i<=19;i++)
{
for(int j=1;j+(1<<i)-1<=n;j++)
{
f[j][i]=__gcd(f[j][i-1],f[j+(1<<(i-1))][i-1]);
}
}
}
int ask(int l,int r)
{
int k=log2(r-l+1);
return __gcd(f[l][k],f[r-(1<<k)+1][k]);
}
void solve(int n)
{
mmp.clear();
for(int i=1;i<=n;i++)
{
int now=f[i][0];
int L=i,R=i;
while(L<=n)
{
int l=L,r=n;
int mid;
while(l<=r)
{
int m=(l+r)>>1;
if(ask(i,m)==now)l=m+1,mid=m;
else r=m-1;
}
mmp[now]+=mid-L+1;
L=mid+1;
now=ask(i,L);
}
}
}
int main()
{
int t;
cin>>t;
int cas=1;
while(t--)
{
int n;
scanf("%d",&n);
for(int i=1;i<=n;i++)scanf("%d",&a[i]);
build(n);
solve(n);
int m;
scanf("%d",&m);
printf("Case #%d:\n",cas++);
for(int i=1;i<=m;i++)
{
int x,y;
scanf("%d%d",&x,&y);
int tmp=ask(x,y);
printf("%d %lld\n",tmp,mmp[tmp]);
}
}
}