数学知识3Accepted12

高斯消元

高斯消元解线性方程组

#include <iostream>
#include <cmath>

using namespace std;

const int N = 110;
const double eps = 1e-6;// 1乘以10的负6次方

int n;
double a[N][N];// 存放增广矩阵

int gauss()
{
    int c, r;// column列 row行
    for (c = 0, r = 0; c < n; c ++ )
    {
        int t = r;
        // 选择绝对值最大的行
        for (int i = r; i < n; i ++ )
        {
            if (fabs(a[i][c]) > fabs(a[t][c]))
            {
                t = i;
            }
        }
        
        if (fabs(a[t][c]) < eps) continue;// 当前行为0
        
        // 将该行换到最上面 
        for (int i = c; i <= n; i ++ ) swap(a[t][i], a[r][i]);// 交换第t列的第i行与第r列的第i行
        
        // 将该行第1个数变成1 
        for (int i = n; i >= c; i -- ) a[r][i] /= a[r][c];// 倒着
        
        // 将下面所有行的该列消成0
        for (int i = r + 1; i < n; i ++ )// 
        {
            if (fabs(a[i][c]) > eps) // 已经是0不用操作，不是0消成0
            {
                for (int j = n; j >= c; j -- )
                {
                    a[i][j] -= a[r][j] * a[i][c];
                }
            }
        }
        
        r ++;
    }
    
    if (r < n)
    {
        for (int i = r; i < n; i ++ )
        {
            if (fabs(a[i][n]) > eps)
            {
                return 2;// 无解
            }
        }
        
        return 1;// 无穷多解
    }
    
    for (int i = n - 1; i >= 0; i -- )
    {
        for (int j = i + 1; j < n; j ++ )
        {
            a[i][n] -= a[i][j] * a[j][n];// 递推出来每一个x 的值
        }
    }
    
    return 0;// 唯一解
}

int main()
{
    cin >> n;
    for (int i = 0; i < n; i ++ )
    {
        for (int j = 0; j <= n; j ++ )// n+1个数
        {
            cin >> a[i][j];
        }
    }
    
    int t = gauss();
    if (t == 0)
    {
        for (int i = 0; i < n; i ++ ) printf("%.2lf\n", a[i][n]);
    }
    else if (t == 1) puts("Infinite group solutions");
    else puts("No solution");
    
    return 0;
}

高斯消元解异或线性方程组

在这里插入代码片

组合计数

求组合数Ⅰ

#include <iostream>

using namespace std;

const int N = 2010, mod = 1e9 + 7;

int c[N][N];

void init()
{
    for (int i = 0; i < N; i ++ )
        for (int j = 0; j <= i; j ++ )
            if (!j) c[i][j] = 1;
            else c[i][j] = (c[i - 1][j] + c[i - 1][j - 1]) % mod;
}

int main()
{
    init();
    
    int n;
    scanf("%d", &n);
    while(n --)
    {
        int a, b;
        scanf("%d%d", &a, &b);
        printf("%d\n", c[a][b]);
    }
    
    return 0;
}