본문 바로가기

알고리즘 관련/BOJ

BOJ)11689 GCD(n,k) = 1

문제:icpc.me/11689


1~n까지의 수중에서 n과 서로소인 수를 구하는 문제이다.


소인수분해를 sqrt(n) 시간에 해결해준 뒤 포함-배제의 원리를 이용하여 서로소의 개수를 세주면된다.


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
#include <cstdio>
#include <algorithm>
#include <vector>
#include <queue>
using namespace std;
typedef long long ll;
ll n;
struct ele{
    ll num, pos, state;
    ele(ll num, ll pos, ll state) :num(num), pos(pos), state(state) {}
};
vector<ll> vt;
int main() {
    scanf("%lld"&n);
    if (n == 1) {
        puts("1");
        return 0;
    }
    ll idx = sqrt(n);
    ll tn = n;
    for (int i = 2; i <= idx; i++) {
        while (!(tn%i)) {
            tn /= i;
            vt.push_back(i);
        }
    }
    if (tn != 1)vt.push_back(tn);
    sort(vt.begin(), vt.end());
    vt.erase(unique(vt.begin(), vt.end()), vt.end());
    queue<ele> qu;
    idx = 0;
    for (auto next : vt) 
        qu.push(ele(next, idx++1));
    ll res = n;
    while (qu.size()) {
        ll here = qu.front().num;
        ll pos = qu.front().pos;
        ll state = qu.front().state;
        qu.pop();
        if (state)
            res -= n / here;
        else
            res += n / here;
        for (int i = pos + 1; i < vt.size(); i++) {
            ll next = here*vt[i];
            if (next <= n) {
                qu.push(ele(next, i, state ^ 1));
            }
            else
                break;
        }
    }
    printf("%lld\n", res);
    return 0;
}
cs


'알고리즘 관련 > BOJ' 카테고리의 다른 글

BOJ)1626 두 번째로 작은 스패닝 트리  (5) 2017.06.12
BOJ)10881 프로도의 선물 포장  (0) 2017.06.11
BOJ)2503 숫자 야구  (1) 2017.06.08
BOJ)1799 비숍  (0) 2017.06.08
BOJ)7154 Job Postings  (0) 2017.06.08