Problem p02574 (C++) — buggy submission, verdict: Time Limit Exceeded

#include <algorithm>
#include <cmath>
#include <fstream>
#include <iostream>
#include <map>
#include <numeric>
#include <queue>
#include <random>
#include <set>
#include <sstream>
#include <stack>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <string>
#include <vector>

using namespace std;

#define I_MAX 2147483647
#define LL_MAX 9223372036854775807
#define ll long long
#define ld long double

struct XX {
  ll a;
  int i;
};

class xxGreater {
public:
  bool operator()(const XX &riLeft, const XX &riRight) const {
    // 第2条件
    if ((riLeft.a) == (riRight.a)) {
      return riLeft.i <
             riRight.i; //<:昇順(小さいものから順番)、>:降順(大きいものから順番)
      // プライオリティキューの場合は > で、top()すると値の小さいものがとれる
    }
    // 第1条件
    return (riLeft.a) > (riRight.a);
  }
};

// map<long long,long long> prime_f(long long n){
//     map<long long,long long>res;
//     for(int i=2;i*i<=n;i++){
//         while(n%i==0){
//             ++res[i];
//             n/=i;
//         }
//     }
//     if(n!=1)res[n]=1;
//     return res;
// }

// int n;
////int dat[2*10000000];
////int dat2[2*10000000];
// int dat[10];
// int dat2[10];
//
// void init(int n_){
//     n=1;
//     while(n<n_)n*=2;
//     for(int i=0;i<2*n-1;i++){
//         dat[i]=0;
//         dat2[i]=0;
//     }
// }
//
// void initset(int k,int a){
//     k+=n-1;
//     dat[k]=a;
//     while(k>0){
//         k=(k-1)/2;
//         dat[k]=dat[k*2+1]+dat[k*2+2];
//     }
// }
//
////[a,b)の間を[l,r]区間で比較しアップデート
////引数のindexに注意
////nは固定。initで計算すみ
////update2(L[i],R[i]+1,0,0,n,D[i]);
// void update2(int a,int b,int k,int l,int r,int v){//v更新値、区間は0-index
//     if(r<=a || b<=l)return;
//     if(a<=l && r<=b){
//         dat[k]+=dat2[k];
//         if(r-l>1){
//             dat2[k*2+1]+=dat2[k]/2;
//             dat2[k*2+1]+=dat2[k]/2;
//         }
//         dat2[k]=v*(r-l);
//         return;
//     }else{
//         update2(a,b,k*2+1,l,(l+r)/2,v);
//         update2(a,b,k*2+2,(l+r)/2,r,v);
//         return;
//     }
// }
//
// int query(int a,int b,int k,int l,int r){
//     if(r<=a || b<=l)return 0;
//     if(a<=l && r<=b){
//         dat[k]+=dat2[k];
//         if(r-l>1){
//             dat2[k*2+1]+=dat2[k]/2;
//             dat2[k*2+1]+=dat2[k]/2;
//         }
//         dat2[k]=0;
//         return dat[k];
//     }
//     else{
//         int vl=query(a,b,k*2+1,l,(l+r)/2);
//         int vr=query(a,b,k*2+2,(l+r)/2,r);
//         return vl+vr;
//     }
// }

// void printb(unsigned int v) {
//     unsigned int mask = (int)1 << (sizeof(v) * CHAR_BIT - 1);
//     do putchar(mask & v ? '1' : '0');
//     while (mask >>= 1);
// }

#ifdef DEBUG
ll a[100];
ll is_prime[100 + 2];
#else
ll a[1000000];
ll is_prime[1000000 + 2];
#endif

ll gcd(ll a, ll b) {
  if (a < b) {
    swap(a, b);
  }
  while (b) {
    ll m = a % b;
    a = b;
    b = m;
  }
  return a;
}
ll prime[1000001];

int main(int argc, const char *argv[]) {
  // scanf("%s",S);
  // scanf("%d",&N);
  // scanf("%lld %lld",&target1,&target2);
  // sscanf(tmp.c_str(),"%dd%d%d",&time[i], &dice[i], &z[i]);
  // getline(cin, target);
  // ifstream ifs("01");//テスト用
  // ifs >> a;
  // ここから

  // 入力高速化
  ios::sync_with_stdio(false);
  cin.tie(0);

  ll n;
  cin >> n;

  for (int i = 0; i < n; i++) {
    cin >> a[i];
  }

  int c[1000002];
  int p = 0;
  for (int i = 0; i <= 1000001; i++) {
    is_prime[i] = 1;
    c[i] = 0;
  }
  is_prime[0] = 1;
  is_prime[1] = 1;
  // この2重ループで十分高速
  for (ll i = 2; i <= 1000001; i++) {
    if (is_prime[i] == 1) {
      is_prime[i] = i;
      for (ll j = 2 * i; j <= 1000001; j += i) {
        if (is_prime[j] == 1) {
          is_prime[j] = i;
        }
      }
      prime[p++] = i;
    }
  }

  int f = 0;
  for (int i = 0; i < n; i++) {
    vector<ll> stock;
    ll tmp = a[i];
    while (tmp != 1) {
      if (c[is_prime[tmp]] > 0) {
        f = 1;
        break;
      } else {
        stock.push_back(is_prime[tmp]);
        tmp = tmp / is_prime[tmp];
      }
    }
    for (int j = 0; j < stock.size(); j++) {
      c[stock[j]]++;
    }
    if (f == 1) {
      break;
    }
  }
  if (f == 0) {
    cout << "pairwise coprime" << endl;
    return 0;
  }
  ll g = a[0];
  for (int i = 1; i < n; i++) {
    g = gcd(g, a[i]);
  }
  if (g == 1) {
    cout << "setwise coprime" << endl;
  } else {
    cout << "not coprime" << endl;
  }

  // ここまで
  // cout << "ans" << endl;
  // cout << " " << "ans" << endl;
  // printf("%.0f\\n",ans);//小数点以下表示なし
  // printf("%.7f\\n",p);

  return 0;
}


Error: TLE

Fixed submission:

#include <algorithm>
#include <cmath>
#include <fstream>
#include <iostream>
#include <map>
#include <numeric>
#include <queue>
#include <random>
#include <set>
#include <sstream>
#include <stack>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <string>
#include <vector>

using namespace std;

#define I_MAX 2147483647
#define LL_MAX 9223372036854775807
#define ll long long
#define ld long double

struct XX {
  ll a;
  int i;
};

class xxGreater {
public:
  bool operator()(const XX &riLeft, const XX &riRight) const {
    // 第2条件
    if ((riLeft.a) == (riRight.a)) {
      return riLeft.i <
             riRight.i; //<:昇順(小さいものから順番)、>:降順(大きいものから順番)
      // プライオリティキューの場合は > で、top()すると値の小さいものがとれる
    }
    // 第1条件
    return (riLeft.a) > (riRight.a);
  }
};

// map<long long,long long> prime_f(long long n){
//     map<long long,long long>res;
//     for(int i=2;i*i<=n;i++){
//         while(n%i==0){
//             ++res[i];
//             n/=i;
//         }
//     }
//     if(n!=1)res[n]=1;
//     return res;
// }

// int n;
////int dat[2*10000000];
////int dat2[2*10000000];
// int dat[10];
// int dat2[10];
//
// void init(int n_){
//     n=1;
//     while(n<n_)n*=2;
//     for(int i=0;i<2*n-1;i++){
//         dat[i]=0;
//         dat2[i]=0;
//     }
// }
//
// void initset(int k,int a){
//     k+=n-1;
//     dat[k]=a;
//     while(k>0){
//         k=(k-1)/2;
//         dat[k]=dat[k*2+1]+dat[k*2+2];
//     }
// }
//
////[a,b)の間を[l,r]区間で比較しアップデート
////引数のindexに注意
////nは固定。initで計算すみ
////update2(L[i],R[i]+1,0,0,n,D[i]);
// void update2(int a,int b,int k,int l,int r,int v){//v更新値、区間は0-index
//     if(r<=a || b<=l)return;
//     if(a<=l && r<=b){
//         dat[k]+=dat2[k];
//         if(r-l>1){
//             dat2[k*2+1]+=dat2[k]/2;
//             dat2[k*2+1]+=dat2[k]/2;
//         }
//         dat2[k]=v*(r-l);
//         return;
//     }else{
//         update2(a,b,k*2+1,l,(l+r)/2,v);
//         update2(a,b,k*2+2,(l+r)/2,r,v);
//         return;
//     }
// }
//
// int query(int a,int b,int k,int l,int r){
//     if(r<=a || b<=l)return 0;
//     if(a<=l && r<=b){
//         dat[k]+=dat2[k];
//         if(r-l>1){
//             dat2[k*2+1]+=dat2[k]/2;
//             dat2[k*2+1]+=dat2[k]/2;
//         }
//         dat2[k]=0;
//         return dat[k];
//     }
//     else{
//         int vl=query(a,b,k*2+1,l,(l+r)/2);
//         int vr=query(a,b,k*2+2,(l+r)/2,r);
//         return vl+vr;
//     }
// }

// void printb(unsigned int v) {
//     unsigned int mask = (int)1 << (sizeof(v) * CHAR_BIT - 1);
//     do putchar(mask & v ? '1' : '0');
//     while (mask >>= 1);
// }

#ifdef DEBUG

#else

#endif

ll gcd(ll a, ll b) {
  if (a < b) {
    swap(a, b);
  }
  while (b) {
    ll m = a % b;
    a = b;
    b = m;
  }
  return a;
}
int prime[1000001];
bool is_prime[1000000 + 2];
ll a[1000000];

int main(int argc, const char *argv[]) {
  // scanf("%s",S);
  // scanf("%d",&N);
  // scanf("%lld %lld",&target1,&target2);
  // sscanf(tmp.c_str(),"%dd%d%d",&time[i], &dice[i], &z[i]);
  // getline(cin, target);
  // ifstream ifs("01");//テスト用
  // ifs >> a;
  // ここから

  // 入力高速化
  ios::sync_with_stdio(false);
  cin.tie(0);

  ll n;
  cin >> n;

  for (int i = 0; i < n; i++) {
    cin >> a[i];
  }

  int p = 0;
  for (int i = 0; i <= 1000001; i++) {
    is_prime[i] = true;
  }
  is_prime[0] = false;
  is_prime[1] = false;
  for (int i = 2; i <= 1000001; i++) {
    if (is_prime[i]) {
      prime[p++] = i;
      for (int j = 2 * i; j <= 1000001; j += i) {
        is_prime[j] = false;
      }
    }
  }
  vector<int> pr;
  for (int i = 2; i <= 1000001; i++) {
    if (is_prime[i]) {
      pr.push_back(i);
    }
  }
  int f = 0;
  for (int i = 0; i < pr.size(); i++) {
    int cnt = 0;
    for (int j = 0; j < n; j++) {
      if (a[j] % pr[i] == 0) {
        cnt++;
      }
      if (cnt > 1) {
        break;
      }
    }
    if (cnt > 1) {
      f = 1;
      break;
    }
  }
  if (f == 0) {
    cout << "pairwise coprime" << endl;
    return 0;
  }
  ll g = a[0];
  for (int i = 1; i < n; i++) {
    g = gcd(g, a[i]);
  }
  if (g == 1) {
    cout << "setwise coprime" << endl;
  } else {
    cout << "not coprime" << endl;
  }

  // ここまで
  // cout << "ans" << endl;
  // cout << " " << "ans" << endl;
  // printf("%.0f\\n",ans);//小数点以下表示なし
  // printf("%.7f\\n",p);

  return 0;
}