結果
問題 | No.1711 Divide LCM |
ユーザー |
![]() |
提出日時 | 2021-08-14 20:56:31 |
言語 | C++17(gcc12) (gcc 12.3.0 + boost 1.87.0) |
結果 |
CE
(最新)
AC
(最初)
|
実行時間 | - |
コード長 | 3,783 bytes |
コンパイル時間 | 30,416 ms |
コンパイル使用メモリ | 5,632 KB |
最終ジャッジ日時 | 2025-01-23 22:14:03 |
ジャッジサーバーID (参考情報) |
judge2 / judge2 |
(要ログイン)
コンパイルエラー時のメッセージ・ソースコードは、提出者また管理者しか表示できないようにしております。(リジャッジ後のコンパイルエラーは公開されます)
ただし、clay言語の場合は開発者のデバッグのため、公開されます。
ただし、clay言語の場合は開発者のデバッグのため、公開されます。
コンパイルメッセージ
コンパイルが30秒の制限時間を超えました
ソースコード
#include <cstdio>#include <cstring>#include <iostream>#include <string>#include <cmath>#include <bitset>#include <vector>#include <map>#include <set>#include <queue>#include <deque>#include <algorithm>#include <complex>#include <unordered_map>#include <unordered_set>#include <random>#include <cassert>#include <fstream>#include <utility>#include <functional>#include <time.h>#include <stack>#include <array>#include <list>#include <atcoder/all>#include "testlib.h"#define popcount __builtin_popcountusing namespace std;using namespace atcoder;typedef long long ll;typedef pair<int, int> P;const int MAXN=100000;const int MAXA=1000000000;const int MAXP=1000000;const int MAXM=200000;bitset<MAXP+1> isprime;void sieve(){for(int i=3; i<MAXP+1; i++, i++) isprime[i]=1;isprime[2]=1;for(int i=3; i<MAXP+1; i++){if(isprime[i]){for(int j=(i<<1); j<MAXP+1; j+=i) isprime[j]=0;}}}ll a[MAXN];vector<P> f[MAXN];int mx[MAXP], p1[MAXP];int main(int argc, char* argv[]){registerValidation(argc, argv);sieve();int n=inf.readInt(1, MAXN);inf.readEoln();vector<int> w;unordered_map<int, int> mp;int ms=0;for(int i=0; i<n; i++){int m=inf.readInt(1, MAXM);inf.readEoln();ms+=m;a[i]=1;for(int j=0; j<m; j++){int p, e;p=inf.readInt(2, MAXP);ensure(isprime[p]);inf.readSpace();e=inf.readInt(1, 100);inf.readEoln();for(int j=0; j<e; j++){a[i]*=p;ensure(a[i]<=MAXA);}f[i].push_back({p, e});mx[p]=max(mx[p], e);}}ensure(ms<=MAXM);set<ll> sta;for(int i=0; i<n; i++) sta.insert(a[i]);ensure(sta.size()==n);inf.readEof();for(int i=2; i<MAXP; i++){if(mx[i]==0) continue;p1[i]=1;for(int j=0; j<mx[i]; j++) p1[i]*=i;w.push_back(p1[i]);}for(int i=0; i<n; i++){if(f[i].size()<w.size()) continue;bool ok=1;for(auto q:f[i]){if(mx[q.first]!=q.second){ok=0; break;}}if(ok){cout<<-1<<endl;return 0;}}for(int i=0; i<n; i++){vector<int> f1;for(auto q:f[i]){if(mx[q.first]==q.second) f1.push_back(p1[q.first]);}int m=f1.size();for(int j=0; j<(1<<m); j++){int c=0, x=1;for(int k=0; k<m; k++){if(j&(1<<k)){c++;x*=f1[k];}}if(c>1) mp[x]++;}}queue<P> que;for(int i=0; i<w.size(); i++){que.push({w[i], i});}auto ans=[&](ll x){int c=0;for(auto q:w){if(x%q==0) c++;}cout<<c<<endl;bool used[100010]={};for(auto q:w){if(x%q!=0) continue;vector<int> v;for(int i=0; i<n; i++){if(used[i]) continue;if(a[i]%q!=0){v.push_back(i);used[i]=1;}}cout<<v.size();for(auto i:v){cout<<" "<<a[i];}cout<<endl;}};while(!que.empty()){auto r=que.front(); que.pop();ll x=r.first;int i=r.second;for(int j=i+1; j<w.size(); j++){ll x1=x*w[j];if(mp.find(x1)==mp.end()){ans(x1);return 0;}que.push({x1, j});}}return 0;}