※ 글쓴이는 취미로 코딩을 익혀보는 사람이라 정확하지 않은 내용을 담고 있을 수 있다 ※

 

이번에 볼 문제는 백준 16532번 문제인 Looking for the Risk Factor이다.
문제는 아래 링크를 확인하자.

https://www.acmicpc.net/problem/16532 

 

16532번: Looking for the Risk Factor

For testing a new cryptographic algorithm, engineers working for a large investment bank need to compute a value they named the Risk Factor of the algorithm. Informally, the Risk Factor is the amount of numbers less than or equal to a certain value N, that

www.acmicpc.net

먼저 2 이상 10만 이하의 수에 대해서 해당 수의 factor 중 가장 큰 factor를 구하자.

이는 int 배열에서 에라토스테네스의 체를 만들듯 2, 3, 5의 배수를 둘러보면서, 매 새로운 factor로 순회할 때 해당 수로 나눠지는 숫자들에 그 숫자를 덮어씌워 기록하는 작업을 하는 것으로 빠르게 구할 수 있다.

 

남은 일은 주어진 쿼리들을 K순으로 정렬한 뒤(offline query), K 이하가 최대 factor인 수들에 가중치 1을 모두 업데이트 한 뒤 구간합을 구하는 세그먼트트리를 관리하는 것으로 문제를 해결하는 것이다.

 

아래는 제출한 소스코드이다.

#include <iostream>
#include <vector>
#include <algorithm>
using namespace std;

int seg[262145];

int update(int L, int R, int qI, int sI) {
	if (qI < L || R < qI) return seg[sI];
	if (L == R) return seg[sI] = 1;
	return seg[sI] = update(L, (L + R) / 2, qI, sI * 2) + update((L + R) / 2 + 1, R, qI, sI * 2 + 1);
}

int rangesum(int L, int R, int qL, int qR, int sI) {
	if (R < qL || qR < L) return 0;
	if (qL <= L && R <= qR) return seg[sI];
	return rangesum(L, (L + R) / 2, qL, qR, sI * 2) + rangesum((L + R) / 2 + 1, R, qL, qR, sI * 2 + 1);
}

int sieve[100001];
vector<pair<int, int>> num; //가장큰primefactor, 원래수
vector<pair<pair<int, int>, int>> query; // {기준primefactor(K), 숫자}, 쿼리idx

int ans[100000];

int main() {
	ios::sync_with_stdio(0);
	cin.tie(0);

	for (int i = 2; i < 100000; i++) {
		if (sieve[i]) continue;
		for (int j = i * 2; j < 100001; j += i) {
			sieve[j] = i;
		}
	}
	for (int i = 2; i < 100001; i++) {
		if (sieve[i] == 0) num.emplace_back(make_pair(i, i));
		else num.emplace_back(make_pair(sieve[i], i));
	}

	sort(num.begin(), num.end());

	int Q; cin >> Q;
	for (int i = 0; i < Q; i++) {
		int N, K; cin >> N >> K;
		query.emplace_back(make_pair(make_pair(K, N), i));
	}
	
	sort(query.begin(), query.end());
	auto iter = num.begin();
	for (int i = 0; i < Q; i++) {
		while (iter != num.end()) {
			if ((*iter).first <= query[i].first.first) {
				update(1, 100000, (*iter).second, 1);
				iter++;
			}
			else break;
		}
		ans[query[i].second] = rangesum(1, 100000, 1, query[i].first.second, 1);
	}

	for (int i = 0; i < Q; i++) cout << ans[i] << '\n';
}

 

728x90

'BOJ' 카테고리의 다른 글

[BOJ 15353 // C++] 큰 수 A+B (2)  (0) 2021.11.22
[BOJ 10986 // C++] 나머지 합  (0) 2021.11.21
[BOJ 17425 // C++] 약수의 합  (0) 2021.11.19
[BOJ 9359 // C++] 서로소  (0) 2021.11.18
[BOJ 11689 // C++] GCD(n, k) = 1  (0) 2021.11.17

+ Recent posts