結果

問題 No.5001 排他的論理和でランニング
ユーザー treeone
提出日時 2018-03-17 02:14:45
言語 C++17(gcc12)
(gcc 12.3.0 + boost 1.87.0)
結果
AC  
実行時間 1,488 ms / 1,500 ms
コード長 2,861 bytes
コンパイル時間 2,149 ms
実行使用メモリ 7,740 KB
スコア 51,354,845
最終ジャッジ日時 2020-03-12 19:53:04
ジャッジサーバーID
（参考情報） 		judge7 /
純コード判定しない問題か言語
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 50
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#define rep(i, a, n) for(int i = a; i < n; i++)
#define REP(i, n) rep(i, 0, n)
#define repb(i, a, b) for(int i = a; i >= b; i--)
#define all(a) a.begin(), a.end()
#define int long long
#define chmax(x, y) x = max(x, y)
#define chmin(x, y) x = min(x, y)
using namespace std;
typedef pair<int, int> P;
const int mod = 1000000007;
const int INF = 1e12;
unsigned w;
unsigned xor128() {
  static unsigned x = 123456789, y = 362436069, z = 521288629/*, w = 88675123*/;
  unsigned t = (x ^ (x << 11));
  x = y; y = z; z = w;
  return (w = (w ^ (w >> 19)) ^ (t ^ (t >> 8)));
}
class Timer {
private:
  clock_t start;
public:
  Timer() {
    start = clock();
  }
  double sec() {
    clock_t end = clock();
    return (double) (end - start) / CLOCKS_PER_SEC;
  }
};
const double TimeLimit = 1.450;
Timer timer;
signed main(){
    ios::sync_with_stdio(false);
    cin.tie(0);
    clock_t start, end;
	start = clock();
	random_device rnd;
	mt19937 mt((int)time(0));
    w = mt();
    int n, m;
    cin >> n >> m;
    vector<int> a(n);
    vector<int> d[20];
    rep(i, 0, n){
        cin >> a[i];
    }
    sort(all(a));
    rep(i, 0, n){
        repb(j, 19, 0){
            if(a[i] & (1LL << j)){
                d[j].push_back(i);
                break;
            }
        }
    }
    // vector<int> ans;
    set<int> st;
    int now = 0;
    repb(i, 19, 0){
        if(st.size() == m) break;
        if(now & (1LL << i)) continue;
        int MAX = now, idx = -1;
        rep(j, 0, d[i].size()){
            if(now ^ a[d[i][j]] > MAX){
                MAX = now ^ a[d[i][j]];
                idx = d[i][j];
            }
        }
        if(idx != -1){
            now = MAX;
            st.insert(idx);
        }
    }
    rep(i, 0, 20){
        rep(j, 0, d[i].size()){
            if(st.size() == m) break;
            if(!st.count(d[i][j])){
                now ^= a[d[i][j]];
                // ans.push_back(i);
                st.insert(d[i][j]);
            }
        }
    }
    cerr << st.size() << " " << now << endl;
    int cnt = 0, grow = 0;
	while(1){
		double tmp = timer.sec();
		if(tmp > TimeLimit){
            break;
        }
        int id = xor128() % n;
        int id2 = xor128() % n;
        if(!(st.count(id) ^ st.count(id2))) continue;
        if(st.count(id2)) swap(id, id2); 
        cnt++;
        int tmp2 = now ^ a[id] ^ a[id2];
        if(tmp2 > now){
            now = tmp2;
            st.erase(id);
            st.insert(id2);
            grow++;
        }
		// if((int)(tmp * 1000) % 1000 <= 20){
		// 	cerr << tmp << " " << now << endl;
		// }
    }
    int sum = 0;
    for(auto it = st.begin(); it != st.end(); it++){
        if(it != st.begin()) cout << ' ';
        cout << a[*it];
        sum ^= a[*it];
    }
    cout << endl;
    cerr << now << " " << sum << " " << cnt << " " << grow << endl;
}
0