結果
問題 | No.382 シャイな人たち (2) |
ユーザー |
|
提出日時 | 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 |
ソースコード
#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;}