컴공 일기 246
… 예… 오랜만에 들어와서 요새 풀고 있는 알고리즘 코드 하나 올려놓고 갑니다…
여전히 꽤…는 아니어도 열심히는 사는 중입니다…
간단한 해설을 하자면 소수를 찾는 알고리즘입니다.
주로 Sieve of Eratosthenes, 에라토스테네스의 체라고 얘기합니다. 컴퓨터에서 소수를 찾는 여러가지 방법이 있습니다만, 학부 수준에서 가장 이해하기 쉽고 직관적인 알고리즘이라고 할 수 있겠네요.
소수가 아닌 게 확실한 수를 지워나가는 방식입니다. 그래서 ‘체’라는 말을 쓰죠. 걸러낸다는 거예요.
그럼 뭘 걸러내면 될까요? “배수”들을 걸러내는 겁니다.
2의 배수, 3의 배수, 4의 배수, 5의 배수, 6의 배수…. 등등 전부요.
예를 잠깐 들어볼까요? 만약에 1부터 100까지의 자연수 범위에서 소수를 찾고 싶으면, sqrt(100) 즉 10의 배수까지 다 지워내면 됩니다. 어? 100까지니까 100의 배수까지 지워야 되는 게 아니냐고요?
사실 맞습니다만, 10의 배수까지만 탐색해도 전부 탐색할 수 있습니다. 만약 N이 소수가 아니라서 a * b의 형태를 이룬다면 즉, N = a * b 라면, a와 b의 최댓값은 루트 N입니다. a와 b는 모두 자연수이기 때문이지요. 만약 둘 중 한 수가 루트 N을 초과하는 순간, a * b의 값은 N을 넘어서게 됩니다. 따라서, 루트N까지 탐색의 범위를 제한해도 무방합니다.
에라토스테네스의 체는 이중반복문으로 구현이 되어서 얼핏 O(N^2)의 부담스런 시간복잡도를 가지는 듯 하지만,
실상은 그렇지 않습니다. 왜냐하면 대다수의 경우가 if(prime[i] == 0) continue;를 충족시키기 때문이지요…
말하자면 그 전에 지워낸 게 많아서, 새로운 배수에서 소수가 아닌 수를 지울 때, 탐색해야 할 후보군이 많이 없다는 뜻입니다.
그 효율성 때문에 알고리즘 문제에서 자주 이용된다.. 뭐 그리 생각하면 됩니다.
오늘도 재미있는 공학 시간..
#include <iostream>
#include <cmath>
#define MAX 1000001
using namespace std;
int prime[MAX];
int n, a, b, result;
int main()
{
cin.tie(NULL);
ios::sync_with_stdio(false);
for(int i=2; i<MAX; i++)
{
prime[i] = i;
}
for(int i=2; i<sqrt(N); i++)
{
If(prime[i] == 0) continue;
for(int j=i * i; j < MAX; j+=i)
{
prime[j] = 0;
}
}
while(1)
{
cin >> n;
if(n == 0) break;
for(int i=3; i < MAX; i++)
{
if(prime[i] != 0)
{
if(n - prime[i] != 0)
{
a = i;
b = n - prime[i];
break;
}
}
}
cout << n << " = " << a << " + " << b << "\n";
}
}
0 XDK (+0)
유익한 글을 읽었다면 작성자에게 XDK를 선물하세요.
-
시불막지마
-
수1 질문 5
양변에 g(x) 나눠버리는거 왜 안되나여 ㅇㅅㅇ
-
현역 최저러고 시대 서바 모고 국어랑 수학 2개 다니고 있는데 둘 다 1?...은...
-
유당불내증 장점 7
변비 걸렸을때 라떼 350ml만 마셔주면 바로 쾌변
-
윤사 해봤는데 1
생윤하면 모래주머니 효과 남? 내신 윤사한다고 자이스토리를 벅벅 풀고 암기한 전적있음
-
그저 평범해지기 위해서 엄청난 에너지를 쓰는 사람들이 있다는 것을 아무도 알지...
-
불행한지도 모름 행복했을 때와 대조되어서 현실이 불행하다고 느끼는 건데 평생을...
-
100일의전사대기중.
-
영어 0
지금 2년째 유기중이고 작수(모의) 84 6모 77…입니다 수능까지 강사컨 딱...
-
낮에는 해커스 강남에서 토플이랑 hsk타고 밤에는 놀러다니고 재밌을듯 그리고 3월부터 다시 자취
-
헉 똥테 달성 6
감사합니다~
-
헉 0
반갑꼬리
-
안정 3이 목표고 7모 낮은 2등급입니다 이 세개 외에 더 좋은 엔제 있으면 추천 부탁드랴용
-
드릴 이미 끝냈는데 엄
-
N수행열차 미리 탑승합니다
-
집중 너무 못하네 요즘
-
대학 '전자책 구독' 추진…전공책 불법복제 막는다 23
[서울경제] 대학가에서 전공서적을 불법 복제하는 관행이 판을 치자 대학들이 전자책...
-
이미지쌤 미친기분 시작편 수원 수투를 플까요 아님 어삼쉬사를 풀까여? 뭐가 더 어렵나요
-
내 대가리가 너무 돌아가버린 탓이다 차라리 환상을 유지했어야 했다 진짜 왜 살지...
-
광역학 치료 아이디어 그대로 가져왔음 6번의 단,-- 조건을 왜 주었는지 연구 중...
-
첫 지문은 가볍게 그럼에도 불구하고 2번 문제는 미묘한 범주 차이를 파악해야 했음
-
확통 공부하면서도 미적 해보지도 않고 확통으로 바꾼게 참.. 바부같다 ㅠㅠ 미적하면 후회하려나~~~
-
정떨어져 1
이제 있던 정 없던 정 다 뗐으니 헤어질 일만 남았군 인간의 이중성이란 무엇인지...
-
김동하 듣는분들 3
이번 세미 잘보심? 난이도가 점점 올라가던데
-
수학 강k 받으려고 권현석T 대기 걸어놨는데 400번대라서.. 다른 수학 선생님들은...
-
수능 2점짜리 난이도?
-
그냥 버티면 심각해지려나? 왼손이랑 차이가 좀 많이 나긴함
-
ㄹㅇ크네
-
고밍이군
-
나만 존나 재미없냐..
-
국어가 2컷이라도 수학 1이면 가능한건가요
-
오부이들 머해 41
공부?
-
전부 알려주새여 ㅠ
-
두과목 번갈아 가면서 하시나요 아니면 같이 하시나요 탐구가 부족한편인데 하루 2시간은 해야겠죠
-
중앙대 영화과에 정시로 꼭 붙고 싶은데ㅠㅠ 다들 어느정도 등급(백분위)로...
-
30분정도를 개지랄을....
-
진심 맨날 띵가띵가 놀아도 수능1 찍는 씹재능러 존나 부럽다
-
진짜모름
-
미적선택인 재수생입니다 제 실력은 6모 59점 7모 61점으로 50-60점대의...
-
수리논술 준비하고 있는데 확통기하 거의 노베입니다 개념 대충 알아요 지금 다니는...
-
10만원대 우산을 본인꺼 마냥 자연스럽게 가져가고 사과 한마디 없이 돌려준...
-
수학 공부 0
제가 지금 미적분을 하는데 하루에 수학9시간투자해서 수1 수2 미적 다 3시간씩...
-
실제로 평가원이나 수능모고에서 p% 퍼센트포인트 단위로 문제나온적있음?? 한번도 없는거같은데
-
제발치려무나 ㅅㅂ
-
궁금한 것 남겨주세요. 세세한 거라도 좋습니다. 대학생활 / 컴공 관련이면...
-
저녁 ㅇㅈ 6
-
파이널때도 수능단단 같은 ㅂㄹㅈ 책 품???
-
트럼프 측근 “한동훈 브라보…美외교정책과 일치하는 훌륭한 방향 제시” 1
엘브리지 콜비 전 미국 국방부 부차관보 韓후보 전당대회 토론회 영상 SNS에 공유...
군대에서 코딩은 어떤 앱으로 하고 계신가요?