結果

問題 No.888 約数の総和
ユーザー d_sei
提出日時 2020-04-11 20:28:43
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 5 ms / 2,000 ms
コード長 1,752 bytes
コンパイル時間 852 ms
コンパイル使用メモリ 101,164 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-09-19 09:29:46
合計ジャッジ時間 1,784 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 30
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<iostream>
#include<string>
#include<algorithm>
#include<vector>
#include<queue>
#include<map>
#include<math.h>
#include<iomanip>
#include<set>
#include<numeric>
#include<cstring>
#include<cstdio>
#include<functional>
#include<bitset>
#include<limits.h>
#include<cassert>
#include<iterator>
#include<complex>
#include<stack>
#include<sstream>
#include<iterator>
using namespace std;
typedef  long long int lint;
typedef pair<int, int> IP;
typedef pair<lint, lint> LLP;
typedef pair<char, char>CP;
#define rep(i, n) for (int i = 0; i < n; i++)
#define repr(i, n) for (int i = n; i >= 0; i--)
#define sort(v) sort((v).begin(), (v).end())
#define reverse(v) reverse((v).begin(), (v).end())
#define upper(v,hoge) upper_bound(v.begin(),v.end(),hoge)
#define lower(v,hoge) lower_bound(v.begin(),v.end(),hoge)
#define llower(v,hoge) *lower_bound(v.begin(), v.end(), hoge)
#define lupper(v,hoge) *upper_bound(v.begin(), v.end(), hoge)

lint name(lint A, lint B) {
	vector<lint>V(B + 1);
	V[0] = 1;
	for (lint i = 1; i <= B; i++) {
		V[i] = V[i - 1] * A;
	}
	lint ans = 0;
	rep(i, V.size()) {
		ans += V[i];
	}
	return ans;
}

vector<lint> prime_factorization(lint N) {
	//素因数分解して配列で返す
	vector<lint> R;
	if (N < 2)return R;
	for (lint i = 2; i * i <= N; i++) {
		while (N % i == 0) {
			R.push_back(i);
			N /= i;
		}
	}
	if (N != 1) {
		R.push_back(N);
	}
	return R;
}

int main() {
	lint N;
	cin >> N;
	vector<lint>A = prime_factorization(N);
	if (A.size() == 0) {
		cout << N << endl;
		return 0;
	}
	lint cnt = 1;
	lint ans = 1;
	for (int i = 1; i < A.size(); i++) {
		if (A[i] == A[i - 1]) {
			cnt++;
		}
		else {
			ans *= name(A[i - 1], cnt);
			cnt = 1;
		}
	}
	ans *= name(A[A.size() - 1], cnt);
	cout << ans << endl;
}
0