给定长度为n的数组a，定义f(l, r) = (a[l]+...+a[r])-(a[l]^...^a[r]),有q次询问：L，R，求对于L＜=l＜=r＜=R, f(l, r)取最大值时l，r的值

若f(l, r) 取最大值时l，r有多种可能，求r - l + 1最小的l，r。

题目

思路：

#include <bits/stdc++.h>
using namespace std;
#define int long long
#define pb push_back
#define fi first
#define se second
#define lson p << 1
#define rson p << 1 | 1
const int maxn = 1e6 + 5, inf = 1e9, maxm = 4e4 + 5, mod = 998244353, N = 1e6;
int a[maxn], b[maxn];
int n, m;
string s;
int sum[maxn];
int Xor[maxn];

void solve(){
    int res = 0;
    // int k;
    int x;
    int q;
    cin >> n;
    cin >> q;
    int L, R;
    // int Xor = 0, sum = 0;
    vector<int> p;
    p.pb(0);
    for(int i = 1; i <= n; i++){
        cin >> a[i];
        if(a[i]){
            p.pb(i);
        }
    }
    int m = p.size() - 1;
    for(int i = 1; i <= m; i++){
        sum[i] = sum[i - 1] + a[p[i]];
        Xor[i] = Xor[i - 1] ^ a[p[i]];
    }
    int num = 0;
    while(q--){
        num++;
        cin >> L >> R;
        if(!m){
            cout << L << ' ' << L << '\n';
            continue;
        }
        int x = lower_bound(p.begin() + 1, p.end(), L) - p.begin();
        int y = upper_bound(p.begin() + 1, p.end(), R) - p.begin();
        y--;
        if(x > y){
            cout << L << ' ' << L << '\n';
            continue;
        }
        L = p[x], R = p[y];
        // cout << x << ' ' << y << '\n';
        int l = p[x], r = p[y];
        int mx = (sum[y] - sum[x - 1]) - (Xor[y] ^ Xor[x - 1]);
        for(int i = x; i <= min(y, x + 32); i++){
            for(int j = max(i, y - 32); j <= y; j++){
                int tmp = (sum[j] - sum[i - 1]) - (Xor[j] ^ Xor[i - 1]);
                if(tmp == mx){
                    if(p[j] - p[i] + 1 < r - l + 1){
                        l = p[i];
                        r = p[j];
                    }
                }
            }
        }
        if(l > r){
            l = L, r = L;
        }
        cout << l << ' ' << r << '\n';
    }
}
    
signed main(){
    ios::sync_with_stdio(0);
    cin.tie(0);
    // fac[0] = 1;
    // for(int i = 1; i <= N; i++){
    //     fac[i] = fac[i - 1] * i % mod;
    // }
    int T = 1;
    cin >> T;
    while (T--)
    {
        solve();
    }
    return 0;
}