数学模板-CSDN博客

本文链接：https://blog.csdn.net/qq_37868325/article/details/81604078

~~1博弈~~

巴什博奕(Bash Game)

只有一堆n个物品，两个人轮流从这堆物品中取物，规定每次至少取一个，最多取m个。最后取光者得胜。要保持给对手留下（m+1）的倍数，就能最后获胜。

威佐夫博奕（Wythoff Game）

有两堆各若干个物品，两个人轮流从某一堆或同时从两堆中取同样多的物品，规定每次至少取一个，多者不限，最后取光者得胜。

前几个奇异局势是：(0,0),(1,2),(3,5),(4,7),(6,10),(8,13),(9,15),(11，18),(12,20).

面对非奇异局势,先拿者必胜;反之,则后拿者取胜。

Fibonacci博弈

有一堆个数为n的石子，游戏双方轮流取石子，满足：

（1）先手不能在第一次把所有的石子取完；

（2）之后每次可以取的石子数介于1到对手刚取的石子数的2倍之间（包含1和对手刚取的石子数的2倍）。约定取走最后一个石子的人为赢家。

当n为Fibonacci数时，先手必败。即存在先手的必败态当且仅当石头个数为Fibonacci数。

尼姆博弈（https://www.cnblogs.com/jiangjun/archive/2012/11/01/2749937.html）

① 有三堆各若干个物品（个数分别为a,b,c），两个人轮流从某一堆取任意多的物品，规定每次至少取一个，多者不限，最后取光者得胜。这里的三个一般可推广至n个

先手必败态等价于a^b^c=0(^表示异或运算）。

② 推广至n个物品：

#include<iostream>
using namespacestd;
int main()
{
    int n,ans,temp;
    while(scanf(“%d”,&n)!=EOF)
    {
        temp=0;
        for(int i=0;i<n;i++)
        {
            scanf("%d",&ans);
            temp^=ans;
        }
        if(temp==0)  cout<<"后手必胜"<<endl;
        else cout<<"先手必胜"<<endl;
    }
    return 0;
}

③ Nim博弈中如果规定最后取光者输，情况是怎样的？

初看起来问题要复杂很多（因为不能主动拿了，而要“躲着”拿，以免拿到最后一个物品），但对于Nim游戏来说，几乎是一样的：

首先按照普通规则一样的策略进行，直到恰好有一个物品数大于1的堆x。在这样的情况下，只需要把堆x中的物品拿得只剩1个物品或者拿完，让对手面临奇数堆物品，这奇数堆物品每堆恰好1个物品。这样的状态显然是必败的。由于你每次操作后需要保证Nim和为0，因此不可能在你操作后首次出现“恰好有一个物品数大于1的堆”。新游戏得到了完美解决。

bool stringGame(int a[], int n) {
    int xor_sum = 0;
    bool flag = false; // 判断是否存在出现次数大于2的字符
    for(int i = 0; i < n; i++) {
        xor_sum ^= count[i];
        flag = flag || (count[i] >= 2);
    }
    return ((flag && xor_sum > 0) || (!flag && xor_sum == 0));
}

~~2.逆元~~

Mod是素数：
(a/b)%mod=(a*pow(b,(mod-2)))%mod;
ll pow_mod(ll n,ll k)
{
    ll res=1;
    while(k>0)
    {
        if(k&1)
            res=res*n%mod;
        n=n*n%mod;
        k>>=1;
    }
    return res;
}

void init_fac()
{
    fac[0]=fac[1]=inv[1]=1;
    for(int i=2;i<MAXN;++i)
    {
        fac[i]=i*fac[i-1]%mod;
        inv[i]=pow_mod(i,mod-2);
    }
}

~~3.筛素数~~

void AA()
{
    int i,j;
    num=0;
    for(i=2;i<M;i++)
    {
        if(!vis[i])
            prime[num++]=i;
        for(j=0;j<num;j++)
        {
            if(i*prime[j]>M)
                break;
            vis[i*prime[j]]=true;
            if(i%prime[j]==0)
                break;
        }
    }
}

~~4.矩阵快速幂~~

struct Matrix
{
    int m[3][3];
};
Matrix Mul(Matrix a,Matrix b)
{
    Matrix c;
    memset(c.m,0,sizeof(c.m));
    for(int i=0;i<3;i++)
        for(int j=0;j<3;j++)
            for(int k=0;k<3;k++)
                c.m[i][j] += ((a.m[i][k]*b.m[k][j])%SMod + SMod)%SMod;
    return c;
}
Matrix fastm(Matrix a,int n) 
{
    Matrix res;
    memset(res.m,0,sizeof(res.m));
    res.m[0][0] = res.m[1][1] = res.m[2][2] = 1;
    while(n)
    {
        if(n&1)
            res = Mul(res,a);
        n>>=1;
        a = Mul(a,a);
    }
    return res;
}

~~5.施罗德数字（+组合数模板）~~

/*
1
1 2
1 4 6
1 6 16 22
1 8 30 68 90

A[i][j]=a[i-1][j]+a[i-1][j-1]+a[i][j-1];

序列：1 2 6 22 90.。
*/

#if 01
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll MOD=998244353;
const ll MAXN=262144+5;
ll a[MAXN];
ll fact[MAXN],ifact[MAXN];//fact[i]是i的阶乘，ifact[i]是阶乘的除法逆元，两者用于求组合数
ll pow_mod(ll n,ll k,ll mod) //快速幂求n^k余m的结果
{
    ll res=1;
    n=n%mod;
    while(k>0)
    {
        if(k&1)
            res=res*n%mod;
        n=n*n%mod;
        k>>=1;
    }
    return res;
}
void init()//初始化
{
    fact[0]=ifact[0]=1;
	for(int i=1;i<MAXN;++i)
    {
        fact[i]=(fact[i-1]*i)%MOD;
		ifact[i]=pow_mod(fact[i],MOD-2,MOD);
	}
}
ll C(ll n,ll m)//求组合数
{
    if(n<m||m<0)
        return 0;
    return  fact[n]*ifact[m]%MOD*ifact[n-m]%MOD;
}
int main()
{
    init();
    ll n;
    while(scanf("%lld",&n)!=EOF)
    {
        if(n==1) {cout<<"1"<<endl;}
        else
        {
            n--;
        ll ans=0;
        for(int k=0;k<=n;k++)
        {
            ans=(ans+(pow_mod(2,k,MOD)*C(n,k)%MOD)*C(n,k-1)%MOD)%MOD;
        }
        cout<<ans*pow_mod(n,MOD-2,MOD)%MOD<<endl;
        }
    }
    return 0;
}
#endif // 1

~~6.判断两个正方形位置关系~~

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const double eps = 1e-6;
int sgn(double x)
{
    if(fabs(x) < eps)return 0;
    if(x < 0)return -1;
    else return 1;
}
struct Point
{
    double x,y;
    Point(){}
    Point(double _x,double _y)
    {
        x = _x;y = _y;
    }
    Point operator -(const Point &b)const
    {
        return Point(x - b.x,y - b.y);
    }
    //叉积
    double operator ^(const Point &b)const
    {
        return x*b.y - y*b.x;
    }
    //点积
    double operator *(const Point &b)const
    {
        return x*b.x + y*b.y;
    }
    //绕原点旋转角度B（弧度值），后x,y的变化
    void transXY(double B)
    {
        double tx = x,ty = y;
        x = tx*cos(B) - ty*sin(B);
        y = tx*sin(B) + ty*cos(B);
    }
};
struct Line
{
    Point s,e;
    Line(){}
    Line(Point _s,Point _e)
    {
        s = _s;e = _e;
    }
    //两直线相交求交点
    //第一个值为0表示直线重合，为1表示平行，为0表示相交,为2是相交
    //只有第一个值为2时，交点才有意义
    pair<int,Point> operator &(const Line &b)const
    {
        Point res = s;
        if(sgn((s-e)^(b.s-b.e)) == 0)
        {
            if(sgn((s-b.e)^(b.s-b.e)) == 0)
                return make_pair(0,res);//重合
            else return make_pair(1,res);//平行
        }
        double t = ((s-b.s)^(b.s-b.e))/((s-e)^(b.s-b.e));
        res.x += (e.x-s.x)*t;
        res.y += (e.y-s.y)*t;
        return make_pair(2,res);
    }
};

//判断线段相交
bool inter(Line l1,Line l2)
{
    return
    max(l1.s.x,l1.e.x) >= min(l2.s.x,l2.e.x) &&
    max(l2.s.x,l2.e.x) >= min(l1.s.x,l1.e.x) &&
    max(l1.s.y,l1.e.y) >= min(l2.s.y,l2.e.y) &&
    max(l2.s.y,l2.e.y) >= min(l1.s.y,l1.e.y) &&
    sgn((l2.s-l1.e)^(l1.s-l1.e))*sgn((l2.e-l1.e)^(l1.s-l1.e)) <= 0 &&
    sgn((l1.s-l2.e)^(l2.s-l2.e))*sgn((l1.e-l2.e)^(l2.s-l2.e)) <= 0;
}

//判断点在线段上
//判断点在线段上
bool OnSeg(Point P,Line L)
{
    return
    sgn((L.s-P)^(L.e-P)) == 0 &&
    sgn((P.x - L.s.x) * (P.x - L.e.x)) <= 0 &&
    sgn((P.y - L.s.y) * (P.y - L.e.y)) <= 0;
}
//判断点在凸多边形内
//点形成一个凸包，而且按逆时针排序（如果是顺时针把里面的<0改为>0）
//点的编号:0~n-1
//返回值：
//-1:点在凸多边形外
//0:点在凸多边形边界上
//1:点在凸多边形内
int inConvexPoly(Point a,Point p[],int n)
{
    for(int i = 0;i < n;i++)
    {
        if(sgn((p[i]-a)^(p[(i+1)%n]-a)) < 0)return -1;
        else if(OnSeg(a,Line(p[i],p[(i+1)%n])))return 0;
    }
    return 1;
}
//判断点在任意多边形内
//射线法，poly[]的顶点数要大于等于3,点的编号0~n-1
//返回值
//-1:点在凸多边形外
//0:点在凸多边形边界上
//1:点在凸多边形内
int inPoly(Point p,Point poly[],int n)
{
    int cnt;
    Line ray,side;
    cnt = 0;
    ray.s = p;
    ray.e.y = p.y;
    ray.e.x = -100000000000.0;//-INF,注意取值防止越界

    for(int i = 0;i < n;i++)
    {
        side.s = poly[i];
        side.e = poly[(i+1)%n];

        if(OnSeg(p,side))return 0;

        //如果平行轴则不考虑
        if(sgn(side.s.y - side.e.y) == 0)
            continue;

        if(OnSeg(side.s,ray))
        {
            if(sgn(side.s.y - side.e.y) > 0)cnt++;
        }
        else if(OnSeg(side.e,ray))
        {
            if(sgn(side.e.y - side.s.y) > 0)cnt++;
        }
        else if(inter(ray,side))
            cnt++;
    }
    if(cnt % 2 == 1)return 1;
    else return -1;
}
int main()
{
    int t;
    //cin >> t;
    int x1, x2, x3, x4, x5, x6, x7,x8, y1,y2,y3,y4,y5,y6,y7,y8;
    //while(t--)
    {
        int flag = 0;
        scanf("%d%d%d%d%d%d%d%d", &x1, &y1, &x2, &y2, &x3,&y3,&x4,&y4);
        scanf("%d%d%d%d%d%d%d%d", &x5, &y5, &x6, &y6,&x7,&y7,&x8,&y8);
        int xx[10],yy[10];
        xx[1]=x1;xx[2]=x2;xx[3]=x3;xx[4]=x4;xx[5]=x5;xx[6]=x6;xx[7]=x7;xx[8]=x8;
        yy[1]=y1;yy[2]=y2;yy[3]=y3;yy[4]=y4;yy[5]=y5;yy[6]=y6;yy[7]=y7;yy[8]=y8;
        for(int i=1;i<=8;i++)
        for(int j=1;j<=8;j++)
        if(i!=j){
            if(xx[i]==xx[j]&&yy[i]==yy[j]) flag=1;
        }
        Point p1[10], p2[10];
        p1[0] = Point(x1, y1);
        p1[1] = Point(x2, y2);
        p1[2] = Point(x3, y3);
        p1[3] = Point(x4, y4);
        p2[0] = Point(x5, y5);
        p2[1] = Point(x6, y6);
        p2[2] = Point(x7, y7);
        p2[3] = Point(x8, y8);
        if(inPoly(p1[0], p2, 4) == 1 &&  inPoly(p1[1], p2, 4) == 1 && inPoly(p1[2], p2, 4) == 1&& inPoly(p1[3], p2, 4) == 1)
            flag = 1;
        if(inPoly(p2[0], p1, 4) == 1 &&  inPoly(p2[1], p1, 4) == 1 && inPoly(p2[2], p1, 4) == 1&& inPoly(p2[3], p1, 4) == 1)
            flag = 1; // limian
        for(int i = 0; i < 4; i++)
        {
            for(int j = i; j < 4; j++)
            {
                Line l1 = Line(p1[i], p1[j]);
                for(int x = 0; x < 4; x++)
                {
                    for(int y = x; y < 4; y++)
                    {
                        Line l2 = Line(p2[x], p2[y]);
                        if(inter(l1, l2))
                            flag = 2;

                    }
                }
            }
        }

        if(flag == 1||flag==2)
            puts("YES");
        else
            puts("NO");

    }
    return 0;
}