推荐一个讲得不错的博客:
洛谷模板题:
模板题代码:
1 #include2 #include 3 #include 4 #include 5 #include 6 #define LL long long 7 #define RI register int 8 #define eps 1e-7 9 using namespace std;10 const int INF = 0x7ffffff ;11 const int N = 100 + 10 ;12 13 inline int read() { // 快读 14 int k = 0 , f = 1 ; char c = getchar() ;15 for( ; !isdigit(c) ; c = getchar())16 if(c == '-') f = -1 ;17 for( ; isdigit(c) ; c = getchar())18 k = k*10 + c-'0' ;19 return k*f ;20 }21 int n ; double hh[N][N] ;22 23 inline bool guass() {24 for(int i=0;i hh[r][i]) r = j ;28 }29 if(fabs(hh[r][i]) < eps) return 0 ; // 因为有精度误差,所以除了x1有无数解时是严格等于0的,其他时候都只是接近0而已 30 if(r != i) for(int j=i;j<=n;j++) swap(hh[i][j],hh[r][j]) ;31 for(int j=i+1;j =0;i--) { // 回带过程 39 for(int j=i+1;j
大家看着代码照着样例手动模拟一下过程就能理解了。 (实践出真知啊qwq)