結果

問題 No.300 平方数
コンテスト
ユーザー newlife171128
提出日時 2018-01-30 00:12:25
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 11 ms / 1,000 ms
コード長 1,092 bytes
コンパイル時間 1,692 ms
コンパイル使用メモリ 160,544 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-12-29 18:19:33
合計ジャッジ時間 2,874 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 43
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <bits/stdc++.h>
#include "bits/stdc++.h"
#include <iostream>
#include <iomanip>
#include <vector>
#include <algorithm>
#include <string>
#include <sstream>
#include <cmath>
#include <stack>
#include <queue>
#include <cctype>
#include <stdio.h>
#include <map>

#include <string.h>
#include <utility>

typedef long long ll;
#define INF (1e9+1)
#define rep(i,n) for(ll i=0;i<(ll)(n);i++)
using namespace std;

typedef pair<int, int> P;

bool erasu[1000001];

vector<int> prim;


ll prime_factorization(ll n){

	ll ans =1;

	/*for(int i=2; i*i<=n; i++)*/
	for(int i=0; i<prim.size(); i++){
		int  cp=0;

		while(n%prim[i] ==0){
			cp++;

			n /=prim[i];
		}

		if(cp%2 !=0){
			ans *= prim[i];
		}
	}

	if( n !=1) ans *=n;

	return ans;

}


void shieve(int n){

	rep(i,n+1){
		erasu[i] =true;
	}

	for(int i=2; i<=n; i++){
		if(erasu[i]){
			prim.push_back(i);
			for(int j= 2*i; j<=n; j +=i){
				erasu[j] = false;
			}
		}
	}
}
int main(){


	ll x=0;
	cin>>x;

	shieve(sqrt(x));
/*	for(int i=0; i<prim.size(); i++){
		cout<<prim[i]<<endl;
	}*/

	cout<<prime_factorization(x)<<endl;
}
0