結果

問題 No.382 シャイな人たち (2)
コンテスト
ユーザー yosupot
提出日時 2016-06-18 00:30:18
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 7,263 ms / 8,000 ms
コード長 2,968 bytes
コンパイル時間 858 ms
コンパイル使用メモリ 105,404 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-10-09 20:08:48
合計ジャッジ時間 161,494 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 21
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <cstdio>
#include <cassert>
#include <cstring>
#include <vector>
#include <valarray>
#include <array>
#include <queue>
#include <set>
#include <unordered_set>
#include <map>
#include <unordered_map>
#include <algorithm>
#include <cmath>
#include <complex>
#include <random>

using namespace std;
using ll = long long;
using ull = unsigned long long;
constexpr int TEN(int n) {return (n==0)?1:10*TEN(n-1);}

ll nex(ll S) {
    return (S * 12345) % 1000003;
}

const int MN = 130;
int N;
vector<int> g[MN];

int idx[MN];

namespace StopWatch {
    clock_t st;
    bool f = false;
    void start() {
        f = true;
        st = clock();
    }
    int msecs() {
        assert(f);
        return (clock()-st)*1000 / CLOCKS_PER_SEC;
    }
}


uint32_t xor96(void) {
  static uint32_t x = 123456789;
  static uint32_t y = 362436069;
  static uint32_t z = 521288629;
  uint32_t t;

  t = (x ^ (x << 3)) ^ (y ^ (y >> 19)) ^ (z ^ (z << 6));
  x = y; y = z;
  return z = t;
}


void solve() {
/*    printf("N: %d\n", N);
    for (int i = 0; i < N; i++) {
        printf("vec %d:", i);
        for (int d: g[i]) {
            printf("%d ", d);
        }
        printf("\n");
    }*/
    int ans = 0;
    vector<int> ans_p;
    iota(idx, idx+N, 0);
    StopWatch::start();
    while (StopWatch::msecs() < 7200) {
        for (int fe = 0; fe < 1000; fe++) {
            for (int i = 0; i < N; i++) {
                swap(idx[i], idx[xor96()%(i+1)]);
            }
            bool used[N];
            fill_n(used, N, false);
            int co = 0;
            for (int i = 0; i < N; i++) {
                int p = idx[i];
                if (!used[p]) {
                    co++;
                    for (int d: g[p]) {
                        used[d] = true;
                    }
                }
            }
            if (ans < co) {
                ans = co;
                ans_p.clear();
                fill_n(used, N, false);
                for (int i = 0; i < N; i++) {
                    int p = idx[i];
                    if (!used[p]) {
                        ans_p.push_back(p+1);
                        for (int d: g[p]) {
                            used[d] = true;
                        }
                    }
                }
            }
        }
    }
    if (ans == N) {
        printf("-1\n");
        return;
    }
    printf("%d\n", ans+1);
    int m = (int)ans_p.size();
    for (int i = 0; i < m; i++) {
        printf("%d", ans_p[i]);
        if (i < m-1) {
            printf(" ");
        } else {
            printf("\n");
        }
    }
}

int main() {
    ll S;
    cin >> S;
    S = nex(S);
    N = (S % 120) + 2;
    S = nex(S);
    ll P = S;
    for (int i = 0; i < N; i++) {
        for (int j = i+1; j < N; j++) {
            S = nex(S);
            if (S >= P) {
                g[i].push_back(j);
                g[j].push_back(i);
            }
        }
    }

    solve();
    return 0;
}
0