高斯消元+叉积求面积+记忆化搜索

Posted 钟钟终

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了高斯消元+叉积求面积+记忆化搜索相关的知识,希望对你有一定的参考价值。

P3389 【模板】高斯消元法

模板采用高斯-约旦消元法
步骤如下:
1.选择一个尚未被选择过的,且为该列的最大系数的未知数作为主元,进行整行交换。
2.再将这一列的其他系数通过行初等变换都变为0。具体做法:该系数除以主元,再用每一行中每一列的值减去这个值乘上求得每一行算出的系数。
3.每一列都这个重复步骤,最后每一行都只存在一个系数和值。做出运算得出值。

#include <bits/stdc++.h>
#define endl '\\n'
#define int long long
#define ios (ios::sync_with_stdio(false),cin.tie(0),cout.tie(0))
#define PII pair<int,int>

using namespace std;
const int N =1e2+5;
const int inf=1e18;
const int mod=1e9+7;
int n;
double d[N][N];
void guess()

    for(int i=1;i<=n;i++)   //枚举列项
    
        int mx=i;
        for(int j=i+1;j<=n;j++)   //选出该列的最大系数
        
            if(abs(d[j][i])>abs(d[mx][i]))
                mx=j;
        
        for(int j=1;j<=n+1;j++)
            swap(d[i][j],d[mx][j]);
        if(!d[i][i])    //最大值为0说明该列都为0,无解
        
            cout<<"No Solution"<<endl;return;
        
        for(int j=1;j<=n;j++)
        
            if(j!=i)
            
                double tmp=d[j][i]/d[i][i];
                for(int k=i;k<=n+1;k++)
                
                    d[j][k]-=d[i][k]*tmp;
                
            
        
    
    for(int i=1;i<=n;i++)
        cout<<fixed<<setprecision(2)<<d[i][n+1]/d[i][i]<<endl;

void solve()

    cin>>n;
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n+1;j++)
            cin>>d[i][j];
    guess();

signed main()

    ios;
    //int T;cin>>T;
    //while(T--)
    solve();
    return 0;



Matrix Equation

抄了大佬们的一个高斯消元求自由元的板子。

#include <bits/stdc++.h>
#define endl '\\n'
#define int long long
#define ios (ios::sync_with_stdio(false),cin.tie(0),cout.tie(0))
#define PII pair<int,int>

using namespace std;
const int N =2e2+5;
const int inf=1e18;
const int mod=998244353;
int n,id;
int A[N][N],B[N][N];
int a[N][N];

int qpow(int x,int y)

    int res=1;
    while(y)
    
        if(y&1) res=res*x%mod;
        x=x*x%mod;
        y>>=1;
    
    return res;

int x[N];       //解集
int freeX[N];   //自由变元
//equ:方程个数 var:变量个数
int Gauss(int equ, int var)

    for(int i = 1; i <= var; i++)
    
        x[i] = 0;
        freeX[i] = 0;
    
    int col = 1;//当前处理的列
    int num = 0;//自由变元的序号
    int row;//当前处理的行
    for(row = 1; row <= equ && col <= var; row++, col++) //枚举当前处理的行
    
        int maxRow = row;//当前列绝对值最大的行
        for(int i = row + 1; i <= equ; i++) //寻找当前列绝对值最大的行
        
            if(abs(a[i][col]) > abs(a[maxRow][col]))
                maxRow = i;
        
        if(maxRow != row) //与第row行交换
        
            for(int j = row; j <= var + 1; j++)
                swap(a[row][j], a[maxRow][j]);
        
        if(a[row][col] == 0) //col列第row行以下全是0,处理当前行的下一列
        
            freeX[num++] = col;//记录自由变元
            row--;
            continue;
        
        for(int i = row + 1; i <= equ; i++)
        
            if(a[i][col] != 0)
            
                for(int j = col; j <= var + 1; j++) //对于下面出现该列中有1的行,需要把1消掉
                
                    a[i][j] ^= a[row][j];
                
            
        
    
    /*求解*/
    //无解:化简的增广阵中存在(0,0,...,a)这样的行,且a!=0
    for (int i = row; i <= equ; i++)
        if (a[i][col] != 0)
            return -1;

    //无穷解: 在var*(var+1)的增广阵中出现(0,0,...,0)这样的行
    if (row <= var)//返回自由变元数
        return var - row+1;//自由变元有var-k个

    //唯一解: 在var*(var+1)的增广阵中形成严格的上三角阵
    for (int i = var ; i >= 1; i--) //计算解集
        x[i] = a[i][var];
        for (int j = i + 1; j <=var; j++)
            x[i] ^= (a[i][j] && x[j]);
    
    return 0;


void solve()

    cin>>n;
    for(int i=1; i<=n; i++) for(int j=1; j<=n; j++)
            cin>>A[i][j];
    for(int i=1; i<=n; i++) for(int j=1; j<=n; j++)
            cin>>B[i][j];
    int ans=0;
    for(int k=1; k<=n; k++)
    
        memset(a,0,sizeof a);
        for(int i=1; i<=n; i++)
        
            for(int j=1; j<=n; j++)
                a[i][j]=A[i][j];
            if(B[i][k]) a[i][i]^=1;
        
        ans+=max(0ll,Gauss(n,n));
    
    cout<<qpow(2,ans)%mod<<endl;

signed main()

    //ios;
    //int T;cin>>T;
    //while(T--)
    solve();
    return 0;





F. Icebergs

逆时针对两个点进行叉乘累加,最后取绝对值再除以2为多边形的面积。
此外,题目中没有要求,先不要考虑给出点的顺时针、逆时针顺序。

#include <bits/stdc++.h>
#define endl '\\n'
#define int long long
#define ios (ios::sync_with_stdio(false),cin.tie(0),cout.tie(0))
#define PII pair<int,int>

using namespace std;
const int N =1e6+5;
const int inf=1e18;
const int mod=1e9+7;
int n,m;
struct Node
	double x, y;
	double operator ^ (const Node &b) const
		return x * b.y - y * b.x;
	
;

vector<Node> e;
void solve()

    cin>>n;
    double ans=0;
    for(int i=1; i<=n; i++)
    
        e.clear();
        cin>>m;
        int lx=inf,rx=0,by=inf,ty=0;
        for(int j=1; j<=m; j++)
        
            int x,y;
            cin>>x>>y;
            e.push_back(x,y);
        
        e.push_back(e[0].x,e[0].y);
        double area=0;
        for(int i=1; i<e.size(); i++)
        
            area+=(e[i-1]^e[i]);
        
        ans+=abs(area)/2;
    
    cout<<(int)(ans)<<endl;

signed main()

    //ios;
    //int T;cin>>T;
    //while(T--)
    solve();
    return 0;



D1. Zero-One (Easy Version)

如果对应位置不一样的数目为奇数,是不存在变换方案的;此外,若只好有一对相邻位置不同的,在x和2*y中取最小值即可;若有偶数前提下有多对,肯定选y值不相邻交换才是最优的

#include<bits/stdc++.h>
#define int  long long
#define endl '\\n'
#define For(i,a,b) for(i=(a);i<=(b);++i)
#define ios (ios::sync_with_stdio(false),cin.tie(0),cout.tie(0))

using namespace std;
const int N=1e6+5;
const int mod=1e9+7;
int n,x,y;
string s1,s2;
void solve()

    cin>>n>>x>>y>>s1>>s2;
    s1=" "+s1;s2=" "+s2;
    int g=0,flag=0;
    for(int i=1;i<=n;i++)
        if(s1[i]!=s2[i]) g++;
    if(g%2)
    
        cout<<-1<<endl;return;
    
    for(int i=1;i<n;i++)
    
        if(s1[i]!=s2[i]&&s1[i+1]!=s2[i+1])
        
            if(g==2)
            
                flag=1;
                break;
            
        
    
    if(flag)
        cout<<min(x,2*y)<<endl;
    else
        cout<<(g/2)*y<<endl;

signed main()

    //ios;
    int T;cin>>T;
    while(T--)
        solve();
    return 0;



D2. Zero-One (Hard Version)

在D1的条件下,加上几组特判和记忆化搜索。
题解:非常棒的大佬题解

#include<bits/stdc++.h>
#define int long long
#define endl '\\n'
#define For(i,a,b) for(i=(a);i<=(b);++i)
#define ios (ios::sync_with_stdio(false),cin.tie(0),cout.tie(0))

using namespace std;
const int inf=1e18;
const int N=5e3+5;
const int mod=1e9+7

以上是关于高斯消元+叉积求面积+记忆化搜索的主要内容,如果未能解决你的问题,请参考以下文章

POJ 2954 /// 皮克定理+叉积求三角形面积

poj 1269 Intersecting Lines——叉积求直线交点坐标

高斯消元法讲解

高斯消元

改革春风吹满地 HDU - 2036 计算几何-----利用叉积计算多边形的面积

POJ 3348 Cows 凸包 求面积