結果

問題 No.1006 Share an Integer
ユーザー KKT89KKT89
提出日時 2020-06-09 02:28:39
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 565 ms / 2,000 ms
コード長 1,085 bytes
コンパイル時間 786 ms
コンパイル使用メモリ 80,544 KB
最終ジャッジ日時 2025-01-11 00:17:56
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 19
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <vector>
#include <algorithm>
using namespace std;
typedef long long int ll;

int main(){
	cin.tie(nullptr);
	ios::sync_with_stdio(false);
	int n; cin >> n;
	vector<int> p,lp(2e6+1);
	// lp[i]: iの最小素因数
	lp[0]=lp[1]=-1;
	for(int i=2;i<lp.size();i++){
		if(lp[i]==0){
			lp[i]=i;
			p.push_back(i);
		}
		for(int j=0;j<p.size()&&p[j]<=lp[i]&&i*p[j]<=lp.size();j++){
			lp[i*p[j]]=p[j];
		}
	}
	// xを素因数分解して返す
	auto pfact=[&](int x){
		vector<pair<int,int>> vec;
		while(x>1){
			if(vec.empty()||vec.back().first!=lp[x]){
				vec.push_back({lp[x],1});
			}
			else{
				vec.back().second++;
			}
			x/=lp[x];
		}
		return vec;
	};
	vector<pair<int,int>> res;
	int ans=1e9;
	for(int i=1;i<n;i++){
		int j=n-i;
		auto p=pfact(i),q=pfact(j);
		int x=1,y=1;
		for(auto pp:p){
			x*=(pp.second+1);
		}
		for(auto qq:q){
			y*=(qq.second+1);
		}
		int cost=abs((i-x)-(j-y));
		if(ans>cost){
			ans=cost;
			res.clear();
		}
		if(ans==cost){
			res.push_back({i,j});
		}
	}
	for(auto v:res){
		printf("%d %d\n",v.first,v.second);
	}
}
0