結果

問題 No.937 Ultra Sword
コンテスト
ユーザー tute7627
提出日時 2019-11-30 10:10:14
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 4,639 bytes
コンパイル時間 1,972 ms
コンパイル使用メモリ 181,248 KB
実行使用メモリ 6,824 KB
最終ジャッジ日時 2024-11-21 03:20:29
合計ジャッジ時間 3,706 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 33 WA * 14
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
using namespace std;

#define endl '\n'
#define ALL(a)  (a).begin(),(a).end()
#define ALLR(a)  (a).rbegin(),(a).rend()
#define spa << " " <<
#define test cout<<"test"<<endl;
#define fi first
#define se second
#define MP make_pair
#define PB push_back
#define EB emplace_back
#define rep(i,n,m) for(ll i = (n); i < (ll)(m); i++)
#define rrep(i,n,m) for(ll i = (m) - 1; i >= (ll)(n); i--)
using ll = long long;
using ld = long double;
const ll MOD = 1e9+7;
//const ll MOD = 998244353;
const ll INF = 1e18;
using P = pair<ll, ll>;
template<typename T>
void chmin(T &a,T b){if(a>b)a=b;}
template<typename T>
void chmax(T &a,T b){if(a<b)a=b;}
void pmod(ll &a,ll b){a=(a+b)%MOD;}
void pmod(ll &a,ll b,ll c){a=(b+c)%MOD;}
void qmod(ll &a,ll b){a=(a*b)%MOD;}
void qmod(ll &a,ll b,ll c){a=(b*c)%MOD;}
ll median(ll a,ll b, ll c){return a+b+c-max({a,b,c})-min({a,b,c});}
void ans1(bool x){if(x) cout<<"Yes"<<endl;else cout<<"No"<<endl;}
void ans2(bool x){if(x) cout<<"YES"<<endl;else cout<<"NO"<<endl;}
void ans3(bool x){if(x) cout<<"Yay!"<<endl;else cout<<":("<<endl;}
template<typename T1,typename T2>
void ans(bool x,T1 y,T2 z){if(x)cout<<y<<endl;else cout<<z<<endl;}  
template<typename T>
void debug(vector<vector<T>>&v,ll h,ll w){for(ll i=0;i<h;i++)
{cout<<v[i][0];for(ll j=1;j<w;j++)cout spa v[i][j];cout<<endl;}};
void debug(vector<string>&v,ll h,ll w){for(ll i=0;i<h;i++)
{for(ll j=0;j<w;j++)cout<<v[i][j];cout<<endl;}};
template<typename T>
void debug(vector<T>&v,ll n){if(n!=0)cout<<v[0];
for(ll i=1;i<n;i++)cout spa v[i];cout<<endl;};
template<typename T>
vector<vector<T>>vec(ll x, ll y, T w){
  vector<vector<T>>v(x,vector<T>(y,w));return v;}
ll gcd(ll x,ll y){ll r;while(y!=0&&(r=x%y)!=0){x=y;y=r;}return y==0?x:y;}
vector<ll>dx={1,0,-1,0,1,1,-1,-1};
vector<ll>dy={0,1,0,-1,1,-1,1,-1};
template<typename T>
vector<T> make_v(size_t a,T b){return vector<T>(a,b);}
template<typename... Ts>
auto make_v(size_t a,Ts... ts){
  return vector<decltype(make_v(ts...))>(a,make_v(ts...));
}
ostream &operator<<(ostream &os, pair<ll, ll>&p){
  return os << p.first << " " << p.second;
}  
const ll MAX_W=30;
struct BitMatrix {
  ll h,w;
  vector<bitset<MAX_W>>A;
  BitMatrix() {};
  BitMatrix(ll n, ll m) :h(n),w(m){
    A.assign(h,bitset<MAX_W>(0));
  };
  BitMatrix(vector<vector<ll>>v):h(v.size()),w(v[0].size()){
    A.assign(h,bitset<MAX_W>(0));
    for(ll i=0;i<h;i++)for(ll j=0;j<w;j++)A[i][j]=(v[i][j]&1);  
  }
  ll height() {
    return h;
  }

  ll width() {
    return w;
  }

  inline bitset<MAX_W> &operator[](ll k) {
    return A[k];
  }
  friend ostream &operator<<(ostream &os, BitMatrix &p) {
    ll n = p.height(), m = p.width();
    for(ll i = 0; i < n; i++) {
      os << "[";
      for(ll j = 0; j < m; j++) {
        os << p[i][j] << (j + 1 == m ? "]\n" : ",");
      }
    }
    return (os);
  }
};
vector<ll>st;
ll GaussJordan(BitMatrix a, bool is_extended){
  ll h = a.height(), w = a.width();
  ll rank = 0;
  ll limit = w - is_extended;
  vector<ll>pre(20);
  rep(i,0,20)pre[i]=i;
  for(ll col=0;col<limit;col++){
    ll pivot = -1;
    for(ll row=rank;row<h;row++){
      if(a[row][col]){
        pivot=row;
        break;
      }
    }
    if(pivot==-1)continue;
    auto tmp=a[pivot][col];
    swap(a[pivot],a[rank]);
    st.PB(col);
    for(ll row=0;row<h;row++){
      if(row!=rank&&a[row][col]!=0){
        a[row]^=a[rank];
      }  
    }
    rank++;
  }
  return rank;
}
int main(){
  cin.tie(nullptr);
  ios_base::sync_with_stdio(false);
  //cout<<fixed<<setprecision(10);
  ll res=0,buf=0;
  bool judge = true;
  ll n;cin>>n;
  vector<ll>a(n);
  rep(i,0,n)cin>>a[i];
  ll minbit=INF;
  vector<ll>num(100001);
  vector<bool>t(100001);
  rep(i,0,n){
    bitset<20>x=a[i];
    res+=a[i];
    num[a[i]]++;
    for(ll j=a[i];j<=100000;j*=2){
      t[j]=true;
    }
  }
  ll cnt=0;
  rep(i,1,100001)if(t[i])cnt++;
  BitMatrix mat(20,cnt);
  cnt=0;
  vector<ll>b;
  rep(i,1,100001){
    if(t[i]){
      //cout<<i<<endl;
      bitset<20>x=i;
      rep(j,0,20)if(x[j])mat[j][cnt]=1;
      b.PB(i);
      cnt++;
    }
  }
  //cout<<mat<<endl;
  ll rank=GaussJordan(mat,false);
  //cout<<mat<<endl;
  vector<bool>t2(150001,false);
  //cout<<rank<<endl;
  //debug(st,st.size());
  rep(i,0,1<<rank){
    bitset<20>x=i;
    ll ret=0;
    rep(j,0,rank)if(x[j])ret^=b[st[j]];
    t2[ret]=true;
  }
  ll maxbuf=0;
  bitset<20>y=a[0];
  rep(i,1,100001){
    if(!t2[i])continue;
    //if(i<=20)cout<<i<<endl;
    ll sumbuf=0;
    for(ll j=i;j<=100000;j+=i){
      sumbuf+=(j-j/i)*num[j];
    }
    //if(i<=30)cout<<i spa sumbuf<<endl;
    chmax(maxbuf,sumbuf);
  }
  cout<<res-maxbuf<<endl;
  return 0;
}
0