結果

問題 No.382 シャイな人たち (2)
コンテスト
ユーザー fumofumofuni
提出日時 2022-09-10 17:52:04
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 2,174 bytes
コンパイル時間 2,059 ms
コンパイル使用メモリ 203,072 KB
最終ジャッジ日時 2025-02-07 03:53:55
ジャッジサーバーID
（参考情報） 		judge2 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 18 WA * 3
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
using namespace std;
#define rep(i,n) for(ll i=0;i<n;i++)
#define repl(i,l,r) for(ll i=(l);i<(r);i++)
#define per(i,n) for(ll i=(n)-1;i>=0;i--)
#define perl(i,r,l) for(ll i=r-1;i>=l;i--)
#define fi first
#define se second
#define pb push_back
#define ins insert
#define pqueue(x) priority_queue<x,vector<x>,greater<x>>
#define all(x) (x).begin(),(x).end()
#define CST(x) cout<<fixed<<setprecision(x)
#define vtpl(x,y,z) vector<tuple<x,y,z>>
#define rev(x) reverse(x);
using ll=long long;
using vl=vector<ll>;
using vvl=vector<vector<ll>>;
using pl=pair<ll,ll>;
using vpl=vector<pl>;
using vvpl=vector<vpl>;
const ll MOD=1000000007;
const ll MOD9=998244353;
const int inf=1e9+10;
const ll INF=4e18;
const ll dy[9]={0,1,0,-1,1,1,-1,-1,0};
const ll dx[9]={1,0,-1,0,1,-1,1,-1,0};
template<class T> inline bool chmin(T& a, T b) {
    if (a > b) {
        a = b;
        return true;
    }
    return false;
}
template<class T> inline bool chmax(T& a, T b) {
    if (a < b) {
        a = b;
        return true;
    }
    return false;
}
using bs=bitset<128>;
vector<int> maxindependent(vector<bs> mat,int trial=1000000){
    int n=mat.size();
    int ans=0;
    bs cand;
    vector<int> order(n);rep(i,n)order[i]=i;
    mt19937 mt(chrono::steady_clock::now().time_since_epoch().count());
    while(trial--){
        shuffle(all(order),mt);
        bs used;
        int ret=0;
        for(auto i:order){
            if((used&mat[i]).any())continue;
            used[i]=1;ret++;
        }
        if(ans<ret){
            ans=ret;
            cand=used;
        }
    }
    vector<int> lc;
    for(int i=0;i<n;i++){
        if(cand[i])lc.emplace_back(i);
    }
    return lc;
}
int main(){
    ll s;cin >> s;
    ll mod=1000003;
    s=s*12345%mod;
    ll n=s%120+2;
    s=s*12345%mod;
    ll p=s;
    vector<bs> mat(n);
    rep(i,n){
        repl(j,i+1,n){
            s=(s*12345)%mod;
            if(s>=p){
                mat[i][j]=1;
                mat[j][i]=1;
            }
        }
    }
    auto f=maxindependent(mat);
    cout << f.size()+1 << endl;
    rep(i,f.size()){
        if(i)cout << " ";
        cout << f[i]+1;
    }
    cout << endl;
}
0