나의 큰 O는 log x야

5차 정사각행렬 같은 경우는 어떻게 역행렬을 구할수 있을까??? 간단한 방법은 5차를 억지로 8차로 늘려버리는 것이다. 남는 자리에 0 (대각선은 1)을 채우면 된다. 그리고 얘의 역행렬을 구하고 왼쪽위에서 5x5개만 찝어내면 그게 원래의 역행렬이다. 그러면 이제 할일은 2^n꼴이 아닌 matrix를 2^n꼴로 확장시키는 것. template class roundbit { public: enum{value=roundbit::value*2}; }; template class roundbit { public: enum{value=2}; }; 2^m꼴이 아닌 수 n을 n이상의 2^m꼴로 바꿔주기 위해서 만든 클래스이다. roundbit::value하면 8이 되고, roundbit::value하면 64가 되겠다..

프로그래밍/테크닉 2009. 6. 28. 04:03

행렬을 작은 크기의 행렬 4개로 나눠어서 계산하는 것에 대해서는 이미 수학자들이 많이 연구해 놓았다. 그러나 글쓴이는 그걸 모르고 공책 한 바닥을 써가며 역행렬 구하는 공식을 구해는 삽질을 했다. (아이구 시간 아까워라.) 무궁한 지식의 원천 위키백과님에게 도움을 요청합니다. 그 공식은 위키백과에 상세히 나와있다. 이제 그럼 본격적으로 역행렬을 구하는 함수를 짜보자. bool has_inverse() { if(a.has_inverse() && (d-c*a.inverse()*b).has_inverse()) { return true; } if(d.has_inverse() && (a-b*d.inverse()*c).has_inverse()) { return true; } return false; } 역시나 te..

프로그래밍/테크닉 2009. 6. 28. 03:41

이번 떡밥은 템플릿 공부도 하고 수학 공부도 할겸 템플릿을 이용해서 n차행렬의 역행렬을 구해보는것이다. c++를 좀 써본 사람이라면 template어쩌구 하는 문법은 알거라고 생각하고 본론으로 들어간다. 이런 행렬을 이렇게 묶어주고 묶어서 생긴 새로운 행렬들을 A,B,C,D라 하면 이렇게 표현할수 있다. (행렬들의 행렬에 대해서도 덧셈, 뺄셈, 곱셈, 실수배가 똑같이 정의된다.) 이런 방법을 재귀적으로 사용한다면 2^n꼴의 행렬을 2*2행렬의 형태로 표현가능하다. 이 사실을 바탕으로 matrix 클래스를 짜본다. template class matrix { protected: enum{halfdim=dim/2}; matrix a,b,c,d; public: matrix() { } matrix(matrix _..

프로그래밍/테크닉 2009. 6. 28. 03:13