#include using namespace std; #if __has_include() #include using namespace atcoder; templateistream &operator>>(istream &is,static_modint &a){long long b;is>>b;a=b;return is;} istream &operator>>(istream &is,modint &a){long long b;cin>>b;a=b;return is;} #endif #ifdef LOCAL #include "debug.h" #else #define debug(...) static_cast(0) #define debugg(...) static_cast(0) templateostream &operator<<(ostream &os,const pair&p){os<; templateusing minque=priority_queue,greater>; templatebool chmax(T &a,const T &b){return (abool chmin(T &a,const T &b){return (a>b?(a=b,true):false);} templateistream &operator>>(istream &is,pair&p){is>>p.first>>p.second;return is;} templateistream &operator>>(istream &is,vector &a){for(auto &i:a)is>>i;return is;} templatevoid operator++(pair&a,int n){a.first++,a.second++;} templatevoid operator--(pair&a,int n){a.first--,a.second--;} templatevoid operator++(vector&a,int n){for(auto &i:a)i++;} templatevoid operator--(vector&a,int n){for(auto &i:a)i--;} #define reps(i,a,n) for(int i=(a);i<(n);i++) #define rep(i,n) reps(i,0,n) #define all(x) x.begin(),x.end() #define pcnt(x) __builtin_popcountll(x) ll myceil(ll a,ll b){return (a+b-1)/b;} template auto vec(const int (&d)[n],const T &init=T()){ if constexpr (id(d,init)); else return init; } template T bs(T ok,T ng,const auto &judge){ static_assert(is_integral::value||is_floating_point::value); static_assert(is_convertible_v>); if constexpr(is_integral::value){ while(abs(ok-ng)>1){ T mid=(ok>>1)+(ng>>1)+(ok&ng&1); (judge(mid)?ok:ng)=mid; } return ok; } else{ for(int i=0;i<120;i++){ T mid=(ok+ng)/2; (judge(mid)?ok:ng)=mid; } return ok; } } void SOLVE(); int main(){ ios::sync_with_stdio(false); cin.tie(nullptr); cout<>testcase; for(int i=0;ileft,right; ll sumleft,sumright; median_min():sumleft(0),sumright(0){} void eval(){ if(left.size()>=right.size()+2){ right.push_front(left.back()); sumright+=left.back(); sumleft-=left.back(); left.pop_back(); } } void add(ll x){ left.push_front(x); sumleft+=x; eval(); } ll getmin()const{ return (left.back()*left.size()-sumleft)+(sumright-left.back()*right.size()); } }; void SOLVE(){ int n; cin>>n; vectorx(n),a(n); cin>>x>>a; vectordp(n+1,1e18); dp[0]=0; reps(i,1,n+1){ ll xorsum=0; median_min mm; for(int j=i-1;j>=0;j--){ xorsum^=a[j]; mm.add(x[j]); chmin(dp[i],dp[j]+xorsum+mm.getmin()); } } cout<