hdu5726

題目大概說給一個包含n個數的序列,多次詢問有多少個區間GCD值等於某個區間的gcd值。

任何一個區間不同的GCD個數是log級別的,因爲隨着右端點向右延伸GCD是單調不增的,而每次遞減GCD至少除以2。

考慮固定左端點,最多就nlogn種GCD,可以直接把所有區間GCD值預處理出來,用map存儲各種GCD值的個數,查詢時直接輸出。

具體是這樣處理的:枚舉左端點,進行若干次二分查找,看當前GCD值最多能延伸到哪兒,進而統計當前GCD值的數量。

而求區間GCD,用ST表,預處理一下,就能在O(1)時間複雜度求出任意區間的gcd了。
注意第二個結果個數有可能是long long, 還有輸出別寫錯了,二分的時候應該是找最右邊的那個, 別加錯了

#include <bits/stdc++.h>
using namespace std;
const int maxn = 1e5 + 10;
int n; int a[maxn];
int dp[maxn][50];
void st(){
    for(int j = 1; (1 << j) <= n; j ++){
        for(int i = 1; i + (1 << j) - 1 <= n ; i ++){
            dp[i][j] = __gcd(dp[i][j - 1], dp[i + (1 << (j - 1))][j - 1]);
        }
    }
}
int query(int l , int r){
    int len = r - l + 1;
    int k = 0 ;
    while((1 << (k + 1)) <= len) k ++;
    return __gcd(dp[l][k], dp[r - (1 << k) + 1][k]);    
}
map<int, long long> mp;
int binary_search(int pos, int l, int r, int num){
    int s = r + 1;
    while(l <= r){
        int mid = (l + r) / 2;
        if(query(pos, mid) >= num){
            l = mid + 1; s = mid;
        }
        else r = mid - 1;
    }
    return s;
}
void solve(){
    for(int i = 1 ; i <= n ; i ++){
        int num = a[i];
        int l = i; int r = n; int pre = i;
        while(l <= n){
            l = binary_search(i, l ,r ,num);
        //  cout << "i = " << i << " l  =" << l << endl;
        //  system("pause");
            mp[num] += (1ll * (l - pre + 1));
            l ++; num = __gcd(num, a[l]);
            pre = l;
        }
    }
}
int main(){
//  freopen("data.out", "r", stdin);
//  freopen("data1.out", "w", stdout);
    int t ;  cin >> t; int cas = 1;
    while(t --){
        scanf("%d", &n); mp.clear();
        for(int i = 1; i <= n ; i++){
            scanf("%d", &a[i]); dp[i][0] = a[i];
        }
        st();
        solve();
        int q; scanf("%d", &q);printf("Case #%d:\n", cas ++);
        while(q --){
            int l, r; scanf("%d%d", &l, &r);
            int gcd = query(l, r);
            printf("%d %lld\n", gcd, mp[gcd]);
        }
    }   
    return 0;
}
發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章