#403125 (C++14) No.937 Ultra Sword

提出ソース

結果

問題	No.937 Ultra Sword
ユーザー	👑 tute7627
提出日時	2019-11-30 10:17:23
言語	C++14 (gcc 12.3.0 + boost 1.83.0)
結果	WA
実行時間	-
コード長	4,618 bytes
コンパイル時間	2,244 ms
コンパイル使用メモリ	181,908 KB
実行使用メモリ	6,824 KB
最終ジャッジ日時	2024-11-21 03:21:21
合計ジャッジ時間	4,401 ms
ジャッジサーバーID （参考情報）	judge1 / judge2

このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	9 ms 6,820 KB
testcase_01	WA	-
testcase_02	AC	6 ms 6,816 KB
testcase_03	AC	9 ms 6,816 KB
testcase_04	AC	4 ms 6,820 KB
testcase_05	AC	13 ms 6,820 KB
testcase_06	WA	-
testcase_07	WA	-
testcase_08	WA	-
testcase_09	WA	-
testcase_10	WA	-
testcase_11	WA	-
testcase_12	WA	-
testcase_13	AC	16 ms 6,820 KB
testcase_14	AC	22 ms 6,820 KB
testcase_15	WA	-
testcase_16	WA	-
testcase_17	WA	-
testcase_18	WA	-
testcase_19	WA	-
testcase_20	WA	-
testcase_21	AC	29 ms 6,816 KB
testcase_22	AC	9 ms 6,816 KB
testcase_23	AC	28 ms 6,816 KB
testcase_24	AC	15 ms 6,816 KB
testcase_25	AC	13 ms 6,816 KB
testcase_26	AC	10 ms 6,816 KB
testcase_27	AC	11 ms 6,816 KB
testcase_28	AC	11 ms 6,816 KB
testcase_29	AC	4 ms 6,816 KB
testcase_30	AC	9 ms 6,820 KB
testcase_31	AC	13 ms 6,820 KB
testcase_32	AC	9 ms 6,816 KB
testcase_33	AC	13 ms 6,816 KB
testcase_34	AC	5 ms 6,820 KB
testcase_35	AC	4 ms 6,816 KB
testcase_36	AC	13 ms 6,820 KB
testcase_37	AC	4 ms 6,816 KB
testcase_38	AC	10 ms 6,816 KB
testcase_39	AC	4 ms 6,820 KB
testcase_40	AC	11 ms 6,820 KB
testcase_41	AC	6 ms 6,816 KB
testcase_42	AC	10 ms 6,816 KB
testcase_43	WA	-
testcase_44	AC	6 ms 6,816 KB
testcase_45	WA	-
testcase_46	AC	8 ms 6,820 KB

権限があれば一括ダウンロードができます

ソースコード

raw source code

#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;
  ll lim=(1<<17);
  vector<ll>num(lim);
  vector<bool>t(lim);
  rep(i,0,n){
    bitset<20>x=a[i];
    res+=a[i];
    num[a[i]]++;
    for(ll j=a[i];j<lim;j*=2){
      t[j]=true;
    }
  }
  ll cnt=0;
  rep(i,1,lim)if(t[i])cnt++;
  BitMatrix mat(20,cnt);
  cnt=0;
  vector<ll>b;
  rep(i,1,lim){
    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(lim,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;
  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;
}

yukicoder

結果

テストケース

ソースコード