結果

問題 No.1514 Squared Matching
ユーザー inksamurai
提出日時 2021-07-03 20:59:51
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 1,404 ms / 4,000 ms
コード長 1,759 bytes
コンパイル時間 1,049 ms
コンパイル使用メモリ 104,900 KB
実行使用メモリ 394,176 KB
最終ジャッジ日時 2024-06-30 13:11:28
合計ジャッジ時間 39,001 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 26
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <cmath>
#include <deque>
#include <algorithm>
#include <iterator>
#include <list>
#include <tuple>
#include <map>
#include <unordered_map>
#include <queue>
#include <set>
#include <unordered_set>
#include <stack>
#include <string>
#include <vector>
#include <fstream>
#include <iostream>
#include <functional>
#include <numeric>
#include <iomanip> 
#include <stdio.h>
#include <assert.h>
//eolibraries
#define lnf 3999999999999999999
#define inf 999999999
#define fi first
#define se second
#define pb push_back
#define ll long long
#define ld long double
#define all(c) (c).begin(),(c).end()
#define sz(c) (int)(c).size()
#define make_unique(a) sort(all(a)),a.erase(unique(all(a)),a.end())
#define pii pair <int,int>
#define ftpii pair<pii,pii>
#define rep(i,n) for(int i = 0 ; i < n ; i++) 
#define drep(i,n) for(int i = n-1 ; i >= 0 ; i--)
#define crep(i,x,n) for(int i = x ; i < n ; i++)
#define vi vector <int> 
#define vec(...) vector<__VA_ARGS__>
#define fcin ios_base::sync_with_stdio(false),cin.tie(0),cout.tie(0);
//eodefine
using namespace std;

const int _maxn = 5e7+3;

int d[_maxn+100],prm[_maxn+10];

void fdprm(int nupto) {
	crep(i,2,sqrt(nupto)+1){
		if(prm[i] == 0) {
			for(int j = i*i ; j <= nupto ; j += i) prm[j] = i;
		}
	}
}

int main(){
fcin;
// 	fdprm(_maxn);

	rep(i,_maxn) d[i]=i;

	crep(p,2,sqrt(_maxn)){
	    if(prm[p]) continue;
	    if(p*p>_maxn) break;
		for(int j=p*p;j<=_maxn;j+=p*p){
			while(d[j]>1 and d[j]%(p*p)==0) d[j]/=(p*p);
		}
		for(int j = p*p ; j <= _maxn ; j += p) prm[j] = p;
	}
	
	rep(i,_maxn) prm[i]=0;
    
    int n;
    cin>>n;
    
    rep(i,n+1){
        prm[d[i]]++;
    }
    
    ll ans=0;
    crep(i,1,n+1){
        ans+=(ll)prm[i]*prm[i];
    }
    
    cout<<ans<<"\n";
/*
*/
	return 0;
}
0