結果
| 問題 |
No.1514 Squared Matching
|
| コンテスト | |
| ユーザー |
umezo
|
| 提出日時 | 2021-05-22 00:54:03 |
| 言語 | C++17(gcc12) (gcc 12.3.0 + boost 1.87.0) |
| 結果 |
TLE
(最新)
AC
(最初)
|
| 実行時間 | - |
| コード長 | 725 bytes |
| コンパイル時間 | 11,584 ms |
| コンパイル使用メモリ | 256,280 KB |
| 最終ジャッジ日時 | 2025-01-21 17:06:06 |
|
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| other | AC * 24 TLE * 2 |
ソースコード
#define rep(i, n) for (int i = 0; i < (int)(n); i++)
#define ALL(v) v.begin(), v.end()
typedef long long ll;
#include <bits/stdc++.h>
using namespace std;
int B[50000010];
int main(){
int n;
cin>>n;
int MAX=50000010;
vector<int> prime(MAX,1);
prime[0]=0,prime[1]=0;
for(int i=2;i*i<=MAX;i++){
if(!prime[i]) continue;
for(int j=2*i;j<=MAX;j+=i) prime[j]=0;
}
vector<int> A;
for(int i=1;i*i<=n;i++){
if(prime[i]) A.push_back(i);
}
rep(i,n+1) B[i]=i;
for(auto p:A){
int q=p*p;
int r=q;
while(1){
for(int i=q;i<=n;i+=q) B[i]/=r;
if(q>n/r) break;
q*=r;
}
}
int ans=0;
for(int i=1;i<=n;i++) ans+=sqrt(n/B[i]);
cout<<ans<<endl;
return 0;
}
umezo