結果

問題 No.917 Make One With GCD
ユーザー Jiro_tech15
提出日時 2019-10-25 22:45:05
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 4 ms / 2,000 ms
コード長 2,815 bytes
コンパイル時間 1,060 ms
コンパイル使用メモリ 93,860 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-06-24 11:05:25
合計ジャッジ時間 1,943 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 32
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <cstring>
#include <algorithm>
#include <vector>
#include <string>
#include <math.h>
#include <iomanip>
#include <limits>
#include <list>
#include <queue>
#include <tuple>
#include <map>
#include <stack>
#include <set>
using namespace std;
#define MOD (long long int)(1e9+7)
#define ll long long int
#define rep(i,n) for(int i=0; i<(int)(n); i++)
#define reps(i,n) for(int i=1; i<=(int)(n); i++)
#define REP(i,n) for(int i=n-1; i>=0; i--)
#define REPS(i,n) for(int i=n; i>0; i--)
#define INF (int)(1123456789)
#define LINF (long long int)(1e18+1)
#define chmax(a, b) a = (((a)<(b)) ? (b) : (a))
#define chmin(a, b) a = (((a)>(b)) ? (b) : (a))
#define all(v) v.begin(), v.end()

const int N = (int)3e5;

ll mpow(ll a, ll b){
  if(b==0) return 1;
  else if(b%2==0){ll memo = mpow(a,b/2); return memo*memo%MOD;}
  else return mpow(a,b-1) * a % MOD;
}
ll lpow(ll a, ll b){
  if(b==0) return 1;
  else if(b%2==0){ll memo = lpow(a,b/2); return memo*memo;}
  else return lpow(a,b-1) * a;
}
ll gcd(ll a, ll b){
  if(b==0) return a;
  else return gcd(b, a%b);
}
vector<ll> kaijo_memo;
ll kaijo(ll n){
  if(kaijo_memo.size() > n) return kaijo_memo[n];
  if(kaijo_memo.size() == 0) kaijo_memo.push_back(1);
  while(kaijo_memo.size() <= n) kaijo_memo.push_back(kaijo_memo[kaijo_memo.size()-1] * kaijo_memo.size() % MOD);
  return kaijo_memo[n];
}
vector<ll> gyaku_kaijo_memo;
ll gyaku_kaijo(ll n){
  if(gyaku_kaijo_memo.size() > n) return gyaku_kaijo_memo[n];
  if(gyaku_kaijo_memo.size() == 0) gyaku_kaijo_memo.push_back(1);
  while(gyaku_kaijo_memo.size() <= n) gyaku_kaijo_memo.push_back(gyaku_kaijo_memo[gyaku_kaijo_memo.size()-1] * mpow(gyaku_kaijo_memo.size(), MOD-2) % MOD);
  return gyaku_kaijo_memo[n];
}

ll nCr(ll n, ll r){
  if(n == r) return 1;//0個の丸と-1個の棒みたいな時に時に効く？不安.
  if(n < r || r < 0) return 0;
  ll ret = 1;
  ret *= kaijo(n); ret %= MOD;
  ret *= gyaku_kaijo(r); ret %= MOD;
  ret *= gyaku_kaijo(n-r); ret %= MOD;
  return ret;
}

multiset<ll> s,kind;

int main(void){
  ll n;cin>>n;
  vector<ll> A;
  rep(i,n){
    ll a;
    cin>>a;
    A.push_back(a);
  }
  ll ans = 0;
  rep(i,n){
    ll a = A[i];
    vector<ll> memo;
    reps(p,(ll)sqrt(a+0.1)+10){
      if(p==1) continue;
      if(a%p == 0){
        while(a%p == 0){
          a /= p;
        }
        memo.push_back(p);
      }
    }
    if(a != 1){
      memo.push_back(a);
    }

    rep(b, lpow(2,memo.size())){
      ll mul = 1;
      ll pc = 0;
      rep(j, memo.size()){
        if((b>>j)&1){
          mul *= memo[j];
          pc++;
        }
      }

      ll count = 0;
      rep(j,n){
        if(j <= i) continue;
        if(A[j]%mul == 0){
          count++;
        }
      }
      ans += lpow(2,count) * (pc%2==0?1:-1);
    }
  }
  cout<<ans<<endl;
  return 0;
}
0