結果

問題 No.937 Ultra Sword
ユーザー 👑 tute7627tute7627
提出日時 2019-11-29 22:59:59
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 4,377 bytes
コンパイル時間 1,815 ms
コンパイル使用メモリ 183,028 KB
実行使用メモリ 6,948 KB
最終ジャッジ日時 2024-05-01 02:44:54
合計ジャッジ時間 4,094 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 11 ms
6,816 KB
testcase_01 WA -
testcase_02 AC 12 ms
6,944 KB
testcase_03 AC 13 ms
6,944 KB
testcase_04 AC 6 ms
6,940 KB
testcase_05 AC 15 ms
6,940 KB
testcase_06 WA -
testcase_07 AC 31 ms
6,940 KB
testcase_08 WA -
testcase_09 WA -
testcase_10 AC 35 ms
6,940 KB
testcase_11 AC 22 ms
6,940 KB
testcase_12 AC 28 ms
6,944 KB
testcase_13 WA -
testcase_14 WA -
testcase_15 WA -
testcase_16 WA -
testcase_17 AC 16 ms
6,940 KB
testcase_18 AC 26 ms
6,940 KB
testcase_19 WA -
testcase_20 WA -
testcase_21 AC 39 ms
6,940 KB
testcase_22 WA -
testcase_23 AC 18 ms
6,944 KB
testcase_24 WA -
testcase_25 WA -
testcase_26 WA -
testcase_27 WA -
testcase_28 WA -
testcase_29 WA -
testcase_30 WA -
testcase_31 WA -
testcase_32 WA -
testcase_33 WA -
testcase_34 WA -
testcase_35 WA -
testcase_36 WA -
testcase_37 WA -
testcase_38 WA -
testcase_39 WA -
testcase_40 WA -
testcase_41 WA -
testcase_42 WA -
testcase_43 WA -
testcase_44 WA -
testcase_45 WA -
testcase_46 WA -
権限があれば一括ダウンロードができます

ソースコード

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=110000;
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);
  }
};

set<ll>st;
ll GaussJordan(BitMatrix &a, bool is_extended){
  ll h = a.height(), w = a.width();
  ll rank = 0;
  ll limit = w - is_extended;
  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];
    st.insert(pivot);
    swap(a[pivot],a[rank]);
    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,0,100001)if(t[i])cnt++;
  BitMatrix mat(20,cnt);
  rep(i,0,100001){
    if(t[i]){
      bitset<20>x=i;
      rep(j,0,20)if(x[j])mat[j][i]=1;
    }
  }
  GaussJordan(mat,false);
  ll maxbuf=0;
  bitset<20>y=a[0];
  rep(i,1,100001){
    bitset<20>x=i;
    ll minbuf=INF;
    judge=true;
    rep(j,0,20){
      if(!st.count(j)&&x[j]!=y[j])judge=false;
    }
    if(!judge)continue;
    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