容斥原理以及Nim基础(异或,SG函数)

容斥原理:

容斥的复杂度为O(2^m),所以可以通过,对于实现,一共2^n-1种,我们可以用二进制来实现

下面是AC代码:

cpp 复制代码
#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int N=20;
int n,m;
int p[N];
int main()
{
    cin>>n>>m;
    for(int i=0;i<m;i++) cin>>p[i];
    int res=0;
    for(int i=1;i<1<<m;i++)
    {
        int t=1,cnt=0;
        for(int j=0;j<m;j++)
        {
            if(i>>j&1)
            {
                cnt++;
                if((LL)t*p[j]>n)
                {
                    t=-1;
                    break;
                }
                t*=p[j];
            }
        }
        if(t!=-1)
        {
            if(cnt%2) res+=n/t;
            else res-=n/t;
        }
    }
    cout<<res;
}

博弈论:

先手必胜状态:可以走到某一个必败状态。

先手必败状态:可以走到任何一个必胜状态。

异或是0就是先手必败,反之必胜。

等价于证a1^a2^...=x!=0,走一步使其变成0

找到x的最高k位1,那么ai中第k位1,那么ai^x<ai,我们从ai拿走ai-(a1^x)个即可。

这样子,假如第一个人!=0,那么下一个人一定是0,那么它又不是0,而一个人败即0^0..=0,因此第一个人一定不会败,即必胜。

sg函数:

首先定义mex(S)表示不属于S的最小非负整数即mex(S)=min(x),x属于自然数并且不属于S

SG函数:

一个状态可以走到不同的状态,这样就形成了有向图,我们令终点为0,SG(终点(即没有出边))=0,假如当前X局面可以通过k个操作变成y1,...yk,那么SG(x)=mex{SG(y1),...SG(yk)}

SG(x)=0:必败,!=0:必胜

因为!=0,那么X可以到0,=0的就到不了0,这样子先手!=0就可以保证自己一定不是0,而终点一定是0,那么先手必胜。

实际情况中可能有不止1个图,每一次可以从其中一个走一步,因此答案就是SG(X1)^SG(X2)...

Xi表示每一个图的起始状态

证明类似与Nim游戏

终点:所有终点SG=0,为0

若其不是0,那么一定存在使其变成0的方案(由mex性质)

假如只有一堆10并只可以选2和5我们可以画出下面的状态图:

下面是AC代码:

cpp 复制代码
#include<bits/stdc++.h>
using namespace std;
const int N=110,M=10010;
int n,m;
int s[N],f[M];
int res;
int sg(int x)//记忆化搜素
{
    if(f[x]!=-1) return f[x];
    unordered_set<int> S;
    for(int i=0;i<m;i++)
    {
        int sum=s[i];
        if (x >= sum) S.insert(sg(x - sum));
    }
    for (int i = 0; ; i ++ )
        if (!S.count(i))
            return f[x] = i;
}
int main()
{
    cin>>m;
    for(int i=0;i<m;i++) cin>>s[i];
    cin>>n;
    memset(f,-1,sizeof(f));//sg是唯一的,因为就算不同图,当前值一样并且可以的操作也一样。
    for(int i=0;i<n;i++)
    {
        int x;
        cin>>x;
        res^=sg(x);
    }
    if(res) cout<<"Yes";
    else cout<<"No";
}

Nim的几个变种:

对于样例,先手把3拿一个给2,变成2,2,2

若对手拿的是第二台阶的,那么它就把相应的放到地面,这样子1和3个数就不变。

若从3拿1给2,我们就从1拿1个给地面,因此对于我们来说,1,3的个数都不一致,而对手都一致,最终都为0,因此先手必胜。

因此,对于一般情况,我们看奇数台阶:

a1^a3^...=x!=0那么先手必胜,反之必败

类似上面的SG函数,不过这里注意SG(b1,b2)=SG(b1)^SG(b2)

cpp 复制代码
#include<bits/stdc++.h>
using namespace std;
const int N=110;
int f[N];
int sg(int x)
{
    if(f[x]!=-1) return f[x];
    unordered_set<int> S;
    for(int i=0;i<x;i++){
        for(int j=0;j<=i;j++)
        {
            S.insert(sg(i)^sg(j));
        }
    }
    for(int i=0;;i++)
    {
        if(!S.count(i)) return f[x]=i;
    }
}
int main()
{
    int n;
    cin>>n;
    memset(f,-1,sizeof(f));
    int res=0;
    for(int i=0;i<n;i++)
    {
        int x;
        cin>>x;
        res^=sg(x);
    }
    if(res) cout<<"Yes";
    else cout<<"No";
}
相关推荐
sc写算法几秒前
基于nlohmann/json 实现 从C++对象转换成JSON数据格式
开发语言·c++·json
SunkingYang4 分钟前
C++中如何遍历map?
c++·stl·map·遍历·方法
Andrew_Xzw6 分钟前
数据结构与算法(快速基础C++版)
开发语言·数据结构·c++·python·深度学习·算法
库库的里昂7 分钟前
【C++从练气到飞升】03---构造函数和析构函数
开发语言·c++
momo卡7 分钟前
MinGW-w64的安装详细步骤(c_c++的编译器gcc、g++的windows版,win10、win11真实可用)
c语言·c++·windows
超的小宝贝1 小时前
数据结构算法(C语言)
c语言·数据结构·算法
凤年徐3 小时前
【数据结构初阶】单链表
c语言·开发语言·数据结构·c++·经验分享·笔记·链表
oioihoii3 小时前
C++11 右值引用:从入门到精通
开发语言·c++
木子.李3477 小时前
排序算法总结(C++)
c++·算法·排序算法
闪电麦坤958 小时前
数据结构:递归的种类(Types of Recursion)
数据结构·算法