hdu 4630 多校第三场1010 树状数组

题意：给你N个数字，每次询问一个区间[ l , r ],输出区间内某两个数的最大公约数。

思路：记录下询问的区间，排序，这里我是按l从大到小排的，方法很多种。

从后开始往前扫，从n - 1，每次扫到一个a[i]，求出他所有的约数。

这里用一个数组vis记录这个约数上次出现的位置，这样就可以保证每次插入这个约数时，这个区间内至少是有两个这个约数的，也就是他至少是两个数的公约数。

然后查询这个区间内的最大值。

理解之后其实蛮水的。

#include <iostream>#include <cstdio>#include <algorithm>#include <string>#include <cmath>#include <cstring>#include <queue>#include <set>#include <vector>#include <stack>#include <map>#include <iomanip>#define PI acos(-1.0)#define Max 2505#define inf 1<<28#define LL(x) ( x << 1 )#define RR(x) ( x << 1 | 1 )#define REP(i,s,t) for( int i = ( s ) ; i <= ( t ) ; ++ i )#define ll long long#define mem(a,b) memset(a,b,sizeof(a))#define mp(a,b) make_pair(a,b)#define PII pair<int,int>using namespace std;inline void RD(int &ret) {    char c;    do {        c = getchar();    } while(c < '0' || c > '9') ;    ret = c - '0';    while((c=getchar()) >= '0' && c <= '9')        ret = ret * 10 + ( c - '0' );}inline void OT(int a) {    if(a >= 10)OT(a / 10) ;    putchar(a % 10 + '0') ;}#define N 50005int n , m ;int a[N] ;int c[N] ;int ans[N] ;struct QU {    int s , e , id ;} Q[N] ;bool cmp(const QU& aa ,const QU& bb) {    return aa.s > bb.s ;}int vis[N] ;void update(int pos ,int num) {    for (int i = pos ; i <= n ; i += i & (-i)){        if(c[i] < num)c[i] = num ;    }}int sum(int pos) {    int ans = 0 ;    for (int i = pos ; i >= 1 ; i -= i & (-i)){        if(ans < c[i])ans = c[i] ;    }    return ans ;}int main() {    int T ;    cin >> T ;    while(T -- ) {        cin >> n  ;        mem(vis, 0) ;        mem(c ,0) ;        REP(i , 1 , n ) {            RD(a[i]) ;        }        cin >> m ;        REP(i , 0 , m - 1) {            RD(Q[i].s) ;            RD(Q[i].e) ;            Q[i].id = i ;        }        sort(Q, Q + m , cmp) ;        int j = 0 ;        int i = n ;        while(1) {            while(i > 0 && Q[j].s <= i) {                for (int j = 1 ; j * j <= a[i] ; j ++ ) {                    if(a[i] % j)continue ;                    int k = a[i] / j ;                    if(vis[j])update(vis[j] , j) ;//vis数组表示这个约数j上次出现的位置，所以上上个位置插入j，因为上次出现过了，所以约数j的数量肯定大于等于2                    if(k != j && vis[k]) {                        update(vis[k] , k) ;                    }                    vis[j] = i ;                    vis[k] = i ;                }                i --;            }            while(Q[j].s > i && j < m) {                ans[Q[j].id] = sum(Q[j].e) ;                j ++ ;            }            if(j >= m)break ;        }        REP(i , 0 , m - 1 ) {            printf("%d\n",ans[i]) ;        }    }    return 0 ;}