結果
| 問題 | 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;
}
            
            
            
        