因为每个操作取负的长度都是小于等于整个区间长度一半,根据字符串周期的dilemma,可以知道所有操作等价于每次能长度g的区间,其中g为所有bi的gcd。
发现这个操作下的一个性质,不管这个操作选在哪个开始位置,选到的区间中都恰好是1个下标mod g = 0, 1个下标mod g = 1,...,1个下标mod g = g - 1,然后考虑把取负的位置i记ci = 1否则ci = 0,能得到不管操作多少次,c[0] ^ c[g] ^ c[2g] ^... = c_1 ^ c[g + 1] ^ c[2g + 1] ^ ... = ...,反过来只要c满足这个性质也就是可以得到的(考虑从左到右不匹配的依次翻转即可)
因此就可以对mod g为某个值的所有下标位置做个dp,最后合并即可。
#include<bits/stdc++.h>
#define pii pair<int,int>
#define fi first
#define sc second
#define pb push_back
#define ll long long
#define trav(v,x) for(auto v:x)
#define all(x) (x).begin(), (x).end()
#define VI vector<int>
#define VLL vector<ll>
#define pll pair<ll, ll>
#define double long double
//#define int long long
using namespace std;
const int N = 1e6 + 100;
const int inf = 1e9;
//const ll inf = 1e18;
const ll mod = 998244353;//1e9 + 7;
#ifdef LOCAL
void debug_out(){cerr << endl;}
template<typename Head, typename... Tail>
void debug_out(Head H, Tail... T)
{
cerr << " " << to_string(H);
debug_out(T...);
}
#define debug(...) cerr << "[" << #__VA_ARGS__ << "]:", debug_out(__VA_ARGS__)
#else
#define debug(...) 42
#endif
void sol()
{
int n, m;
cin >> n >> m;
VLL a(n);
for(int i = 0; i < n; i++)
cin >> a[i];
ll d;
for(int i = 0; i < m; i++)
{
if(i == 0)
cin >> d;
else
{
ll x;
cin >> x;
d = __gcd(d, x);
}
}
VLL ans(2, 0);
for(int t = 0; t < d; t++)
{
VLL dp(2, 0);
dp[1] = -1e18;
for(int i = t; i < n; i += d)
{
VLL dq(2);
dq[0] = max(dp[0] + a[i], dp[1] - a[i]);
dq[1] = max(dp[0] - a[i], dp[1] + a[i]);
swap(dp, dq);
}
ans[0] += dp[0];
ans[1] += dp[1];
}
cout << max(ans[0], ans[1]) << '\n';
}
signed main()
{
ios::sync_with_stdio(0);
cin.tie(0);
int tt;
cin >> tt;
while(tt--)
sol();
}