結果
| 問題 |
No.826 連絡網
|
| コンテスト | |
| ユーザー |
ahe100
|
| 提出日時 | 2019-05-03 23:00:54 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 25 ms / 2,000 ms |
| コード長 | 1,268 bytes |
| コンパイル時間 | 3,856 ms |
| コンパイル使用メモリ | 169,456 KB |
| 実行使用メモリ | 12,364 KB |
| 最終ジャッジ日時 | 2024-11-24 02:32:55 |
| 合計ジャッジ時間 | 3,323 ms |
|
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 2 |
| other | AC * 30 |
ソースコード
#include <bits/stdc++.h>
using namespace std;
#define rep(i,n) for(int i = 0;i<((int)(n));i++)
#define reg(i,a,b) for(int i = ((int)(a));i<=((int)(b));i++)
#define irep(i,n) for(int i = ((int)(n)-1);i>=0;i--)
#define ireg(i,a,b) for(int i = ((int)(b));i>=((int)(a));i--)
typedef long long ll;
/*
2に到達できれば勝ち
*/
ll n,p,ans=0;
vector<ll> prime,fact;
void prime_list(ll n,vector<ll>& p){
ll visited[1000010]={};
reg(i,2,n){
if(visited[i]==0){
p.push_back(i);
reg(j,2,n){
if(j*i>n)break;
visited[i*j]=1;
}
}
}
}
void factoring(ll n,vector<ll>& v,vector<ll>& p){
rep(i,p.size()){
while(n%p[i]==0){
n/=p[i];
v.push_back(p[i]);
}
}
if(n>1)v.push_back(n);
}
void init(){
cin>>n>>p;
prime_list(1000000,prime);
factoring(p,fact,prime);
}
int main(void){
init();
if(p==1){
ans=1;
}else if(fact.size()==1){
if(p*2>n){
ans=1;
}else{
ans=1;
//2がとれるので、n/2以下のすべての素数がOK
rep(i,prime.size()){
if(prime[i]>n)break;
if(prime[i]*2>n)ans++;
}
ans=n-ans;
}
}else{
ans=1;
//2がとれるので、n/2以下のすべての素数がOK
rep(i,prime.size()){
if(prime[i]>n)break;
if(prime[i]*2>n)ans++;
}
ans=n-ans;
}
cout<<ans<<endl;
return 0;
}
ahe100