维护区间gcd(线段树 || ST表+二分)

来源:互联网 发布:ipad 壳淘宝店铺 编辑:程序博客网 时间:2024/06/04 19:05

题目:

因为hdu不知道为什么上不去了所以只有题意啦,不存在的 hdu5726

题意:

给你n个数,然后有m个询问,每个询问一个区间,问你这个区间的gcd是多少,并且问你从1到n有多少个区间的gcd和这个区间的gcd是相同的。

题解:

第一问

gcd都是很好求的,因为两个会重合的区间一样可以求出正确的gcd

第二问

只要把所有区间的gcd求出来丢进map里就好了
对于从l到r这个区间,如果固定l区间,r区间不断变化的时候,gcd的值在不断下降,这个下降速度是很快的。那么只需要枚举左端点l,然后二分l到n的区间,找到最大公约数为gcd的最大区间,然后记录这个gcd的数量,并且把gcd更新

for(int i = 1; i <= n; ++i){       int j = i, g = a[i];       while(j <= n){       int l=j, r = n;            while(l < r){            int mid = (l + r) >> 1;            if(Find(l, r) == g) l = mid + 1;             else r = mid - 1;              }       mp[g]+=l-j+1;       j=l+1;       g=gcd(g,a[j]);       }  } 

其实不用二分也可以,l[i]表示左端点固定,右端点为i,可以跳到多早的j使j~i的gcd都是gcd(l,i)

for(int i=1;i<=n;i++)     for(v[i]=a[i],j=l[i]=i ; j ; j=l[j]-1)      {                   v[j]=gcd(v[j],a[i]);                while(l[j]>1 && gcd(a[i],v[l[j]-1])==gcd(a[i],v[j]))                  l[j]=l[l[j]-1];                mp[v[j]]+=j-l[j]+1;        }    

其实本质上是一样的?

原创粉丝点击