본문 바로가기

[Project Euler 123] (p-1)^n + (p+1)^n % p^2 http://projecteuler.net/index.php?section=problems&id=123 pn 이 n번째 소수 ( 2, 3, 5, 7, 11, ..., )이고, r이 (pn1)n + (pn+1)n 를 pn2으로 나눈 나머지라고 하자. 예를 들어, n = 3, p3 = 5 이고, 43 + 63 = 280 5 mod 25 이다. 이 나머지가 109을 넘는 가장 작은 n은 7037이다. 나머지가 1010을 넘는 가장 작은 n을 구하라. #!/usr/bin/env python # Project Euler Problem 123 def remainder0(p, n): return ((p-1)**n + (p+1)**n)%(p*p) def remainder(p, n): ''' p is a prime. r..
[Project Euler 148] 파스칼 삼각형 mod 7 http://projecteuler.net/index.php?section=problems&id=148 파스칼 삼각형을 7번째 줄까지 써보면, 7로 나누어 떨어지는 것이 없다는 걸 알 수 있다.: 1 1 1 1 2 1 1 3 3 1 1 4 6 4 1 1 5 10 10 5 1 1 6 15 20 15 6 1 그러나, 100번째 줄까지를 보면, 5050개의 숫자 중에서 2361개만 7로 나누어 떨어지지 않는다. 10억(109)번째 줄까지 중에서 7로 나누어 떨어지지 않는 것의 갯수를 구하라. 우선, a+b mod n = (a mod n) + (b mod n) mod n 과 같다는 걸 이용해서 그림을 그려보면, 아주 예쁜 프랙탈모양을 이룸을 알 수 있다. 그리고, 숫자를 7진수로 변환한 후 생각하면 훨씬 생각이..
[MFC|CPP] 사구모양의 포텐셜 만들기 다음 신지식에 재미있어 보이는 게 있길래 답변달다가 만들어 본 것. 2차원 평면상에 주어진 점 P를 중심으로 둥그런 모양의 산을 만들어 보란다. 3차원 그래픽까지 구현할라면 죽을 것 같아서, 2차원 평면의 각 점에서의 함수값은 색으로 표현했다. 단색으로 표현할라니까 분해능이 256가지 뿐이다. 왼쪽 버튼을 클릭하면, 분지의 중앙점이 되는 P가 바뀌고, 오른쪽 버튼을 클릭하면 산이 바라보는 타겟 T가 바뀐다. 휠을 돌리면 분지의 반지름이 넓어진다. 소스코드다. VC6.0 에서 만들었다.
[ProjectEuler 206] 1_2_3_4_5_6_7_8_9_0 = n*n 인 유일한 n구하기. #!/usr/bin/env python ''' for x in xrange(101010101, 138902663): if xx % 10 != 7 or xx % 10 != 3: continue xx = 10*x s = str(xx*xx) # print s # print "1_2_3_4_5_6_7_8_9_0" MeetCondition=True for i in xrange(1, 10): if s[2*(i-1)] != str(i): MeetCondition=False break if MeetCondition: print xx*xx, xx print "1_2_3_4_5_6_7_8_9_0" ''' s = "1_2_3_4_5_6_7_8_9_0" def GetList(lst, length): if length == 9..
직육면체가 특정 면으로 착지할 확률 관련 포스팅 http://astralepic.egloos.com/4022771 http://zariski.egloos.com/2229177 http://mcfrog.org/tt/848 #!/usr/bin/env python import math def length(V): sq_sum = 0 for i in xrange(3): sq_sum += V[i]*V[i] return math.sqrt(sq_sum) def dot(U, V): d = 0 for i in xrange(3): d += U[i]*V[i] return d def omega(A, B, C): ''' computes the steradian defined by three vectors A, B, C. ref : http://en.wikipedi..
[복잡] For all e in G, and for all f not in G, e # G is significantly bigger than f # G.By above definition, we should be able to define a group G out of a linked structure.Remainig problems are how should we define a product of an element on a set of elements, denoted above with #.how should we mathematically define the language "significantly bigger".Idea for 1 e # G = ratio of links from e into GIde..
비트 인버트는 선형적인가? int main() { unsigned short i; unsigned short inv; unsigned short ineg; for(i = 0 ; i < 0xFFFF; i++) { inv = ~i; ineg = 0xFFFF - i; if(inv != ineg) { printf("For %d, ~i (%d) != 0xFFFF - i (%d)\n", i, inv, ineg); } } return 0; } 지금 생각하면 너무 당연하게 둘이 같아야 하는데, 잠결에 이게 아닐 것 같아서 일찍 출근했다. i가 선형적으로 변할 때, 그 인버트된 값도 선형적이다.
[C#] 페졸드책연습중에 한글로 장난 재밌다. ^^;