結果

問題 No.300 平方数
コンテスト
ユーザー xxxasdfghjk
提出日時 2018-09-04 23:14:10
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
WA  
実行時間 -
コード長 2,027 bytes
コンパイル時間 2,245 ms
コンパイル使用メモリ 178,476 KB
実行使用メモリ 6,820 KB
最終ジャッジ日時 2024-11-06 13:05:59
合計ジャッジ時間 4,085 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 27 WA * 16
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
#define MAX_N 100001
#define INF_INT 2147483647
#define INF_LL 9223372036854775807
#define REP(i,n) for(int i=0;i<(int)(n);i++)
using namespace std;
typedef long long int ll;
typedef pair<ll,ll> P;
void init(int n);
int find(int n);
void unite(int x,int y);
bool same(int x, int y);
ll bpow(ll,ll,ll);
typedef vector<int> vec;
typedef vector<vec> mat;
mat mul(mat &A,mat &B);
mat pow(mat A,ll n);
int dx[4] = {1,0,0,-1};
int dy[4] = {0,1,-1,0};

ll gcd(ll a,ll b){
  return a==0 ? b : gcd(b%a,a);
}

ll lcm(ll a,ll b){
  return a*b/gcd(a,b);
}

int main()
{
  ll X,d,Y=1;
  cin >> X;
  d = X;
  set<int> sets;
  map<int,int> mp;
  ll sqr = sqrt(X)+1;
  for(int i=2;i<=sqr;i++){
    while(X%i == 0){
      X /= i;
      sets.insert(i);
      mp[i]++;
    }
  }
  if(X != 1){
    sets.insert(X);
    mp[X]++;
  }
  for(auto itr=sets.begin();itr!=sets.end();++itr){
    Y *= mp[*itr] % 2 ? *itr : 1;
  }
  cout << Y << endl;
  return 0;
}

int par[MAX_N];
int ranks[MAX_N];

//n要素で初期化
void init(int n){
  REP(i,n){
    par[i] = i;
    ranks[i] = 0;
  }

}

//木の根を求める
int find(int x){
  if(par[x] == x){
    return x;
  }else{
    return par[x] = find(par[x]);
  }
}

void unite(int x,int y){
  x = find(x);
  y = find(y);
  if(x == y) return ;
  if(ranks[x] < ranks[y]){
    par[x] = y;
  }else{
    par[y] = x;
    if(ranks[x] == ranks[y]) ranks[x]++;
  }
}

bool same(int x, int y){
  return find(x) == find(y);
}

ll bpow(ll a, ll n,ll mod){
  int i = 0;
  ll res=1;
  while(n){
    if(n & 1)
      res = (res*a) % mod;
    a = (a*a) % mod;
    n >>= 1;
  }
  return res;
}
const int MOD = 1000000007;
mat mul(mat &A, mat &B){
  mat C(A.size(),vec(B[0].size()));
    REP(i,A.size())REP(k,B.size())REP(j,B[0].size()){
    C[i][j] = (C[i][j] + A[i][k] * B[k][j]) % MOD;
  }
  return C;
}

mat pow(mat A,ll n)
{
  mat B(A.size(),vec(A.size()));
  REP(i,A.size()){
    B[i][i] = 1;
  }
  while(n > 0){
    if ( n & 1) B = mul(B,A);
    A = mul(A,A);
    n >>= 1;
  }
  return B;
}
0