BZOJ4589. Hard Nim

BZOJ4589. Hard Nim

题意:

Claris和NanoApe在玩石子游戏,他们有n堆石子,规则如下:

  1. Claris和NanoApe两个人轮流拿石子,Claris先拿。
  2. 每次只能从一堆中取若干个,可将一堆全取走,但不可不取,拿到最后1颗石子的人获胜。
    不同的初始局面,决定了最终的获胜者,有些局面下先拿的Claris会赢,其余的局面Claris会负。
    Claris很好奇,如果这n堆石子满足每堆石子的初始数量是不超过m的质数,而且他们都会按照最优策略玩游戏,那么NanoApe能获胜的局面有多少种。
    由于答案可能很大,你只需要给出答案对10^9+7取模的值。

题解:

首先要知道nim先手获胜条件是所有堆的数量异或为0
那么问题就抽象为:n个数,每个数取值范围是[2,m]中的质数,可以取重,问一共有多少种方案,使得这n个数异或为0
对于每一个2~m的质数p,我们都可以取,用数组b来存就是b[p]=1,质数都被标记为1
CkC_kCk表示异或和为k的方案数有多少种,b[i]=1说明第i位是质数,否则不是
有:Ck=∑i1⊕i2⊕.....⊕in=kbi1×bi2×....×binC_k=\sum_{i_1⊕i_2⊕.....⊕i_n=k}b_{i_1}×b_{i_2}×....×b_{i_n}Ck=i1i2.....in=kbi1×bi2×....×bin
因为我们要求异或和为0,所以k等于0
式子就是:
求n个序列,其中i1⊕i2⊕.....⊕in=0i_1⊕i_2⊕.....⊕i_n=0i1i2.....in=0
C0=∑i1⊕i2⊕.....⊕in=0bi1×bi2×....×binC_0=\sum_{i_1⊕i_2⊕.....⊕i_n=0}b_{i_1}×b_{i_2}×....×b_{i_n}C0=i1i2.....in=0bi1×bi2×....×bin
很明显这就是FWT能解决的问题
但是本题中n很大,n<=1e9,如果直接乘一定会超时,我们仔细观察,这个式子就是n个b序列,而所有b序列都是一样的,相当于其实是n个b相乘,也就是bnb^nbn,所以可以通过快速相乘,结合FWT模板理解一下就懂了

代码:

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <map>
#include <queue>
using namespace std;
typedef long long ll;
typedef int itn;
typedef pair<int, int>PII;
const int N = 5e5 + 7, mod = 1e9 + 7;
const double PI = acos(-1.0);int n, m, limit, t;
int a[N], b[N];
bool vis[N];
int primes[N], cnt;int qpow(int a, int b)
{int res = 1;while(b) {if(b & 1) res = 1ll * res * a % mod;a = 1ll * a * a % mod;b >>= 1;}return res;
}int inv2 = qpow(2, mod - 2);void init(int n)
{for(int i = 2; i <= n; ++ i) {if(vis[i] == 0) primes[ ++ cnt] = i;for(int j = 1; j <= cnt && i * primes[j] <= n; ++ j) {vis[i * primes[j]] = true;if(i % primes[j] == 0) break;}}
}void XOR(int *a, int n, int type = 1)
{for(int o = 2; o <= n; o <<= 1) {for(int i = 0, k = o >> 1; i < n; i += o) {for(int j = 0; j < k; ++j) {int X = a[i + j];int Y = a[i + j + k];a[i + j] = (1ll * X + Y) % mod;a[i + j + k] = ((1ll * X - Y) % mod + mod) % mod;if(type == -1) {a[i + j] = (1ll * a[i + j] * inv2) % mod;a[i + j + k] = (1ll * a[i + j + k] * inv2) % mod;}}}}
}void solve()
{memset(b, 0, sizeof b);memset(a, 0, sizeof a);for(int i = 2; i <= m; ++ i)if(vis[i] == 0) b[i] = 1;limit = 1;while(limit <= m) limit <<= 1;XOR(b, limit);
//    for(int i=0;i<=limit;i++){
//    	cout<<"a[i]="<<a[i]<<endl;
//	}
//	for(int i=0;i<=limit;i++){
//    	cout<<"b[i]="<<b[i]<<endl;
//	}for(int i = 0; i <= limit; ++ i) {b[i] = 1ll * qpow(b[i], n) % mod;}XOR(b, limit, -1);printf("%d\n", b[0]);return ;
}int main()
{
//	freopen("data.in", "r", stdin);init(N - 7);while(scanf("%d%d", &n, &m) != EOF) {solve();}
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/315862.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

P2825 [HEOI2016/TJOI2016]游戏

题目描述 详见 P2825 [HEOI2016/TJOI2016]游戏。 solution 套路题。 一般思路是行列建点跑二分图最大匹配。 此题中的#会分隔行列&#xff0c;因此我们把每行的极大联通块设为点&#xff0c;列同理建点。 再对于每个*对应的行的极大联通块与列的最大联通块编号连边&…

P3033 [USACO11NOV]牛的障碍Cow Steeplechase

P3033 [USACO11NOV]牛的障碍Cow Steeplechase 题目描述 详见&#xff1a;P3033 [USACO11NOV]牛的障碍Cow Steeplechase solution 裸题。 对于每一对相交的线段&#xff0c;连一条边&#xff0c;表示这一对线段不能都取。 显然这样会建立一个二分图&#xff08;横的线段不…

Tree Cutting HDU - 5909

Tree Cutting HDU - 5909 题意&#xff1a; 一个无根树&#xff0c;n个点&#xff0c;n-1条边&#xff0c;每个节点有一个权值&#xff0c;一棵树的权值就是其节点(包含本身及其子节点)的权值的异或和&#xff1b;求价值为[0,m)的树有多少颗&#xff1f;(所谓的树其实就是原连…

Redis作者antirez:开源维护者的挣扎

这两天&#xff0c;一篇名为《开源维护者的挣扎》的文章被迅速顶至 Hacker News 首页&#xff0c;这是 Redis 作者 antirez 发布的最新博客。几个月前&#xff0c;一名开源项目的维护者向 antirez 发邮件&#xff0c;倾诉自己苦心维护项目多年&#xff0c;这或多或少带来了一些…

CF662C Binary Table

CF662C Binary Table 题意&#xff1a; 有一个 n 行 m 列的表格&#xff0c;每个元素都是 0/1 &#xff0c;每次操作可以选择一行或一列&#xff0c;把 0/1 翻转&#xff0c;即把 0 换为 1 &#xff0c;把 1 换为 0 。请问经过若干次操作后&#xff0c;表格中最少有多少个 1 …

P4111 [HEOI2015]小Z的房间

P4111 [HEOI2015]小Z的房间 题目描述 详见&#xff1a;P4111 [HEOI2015]小Z的房间 solution Matrix-Tree裸题。 题意为求一个有障碍的网格图的生成树个数。 因此对于每一个点向左边和上面的节点连边&#xff0c;Matrix-Tree定理求解生成树个数即可。 #include<bits/s…

黑科技抢先尝(续) - Windows terminal中WSL Linux 终端的极简美化指南

之前&#xff0c;本人写了两篇文章 黑科技抢先尝 | Windows全新终端初体验(附代码Build全过程) 和 程会玩 | 无需自行编译也能玩转 Windows Terminal&#xff0c;介绍了玩转Windows terminal的两种方式。今天这篇文章&#xff0c;主要介绍如何美化 Windows terminal 中 WSL 的 …

P3203 [HNOI2010]弹飞绵羊

P3203 [HNOI2010]弹飞绵羊 题目描述 详见&#xff1a;P3203 [HNOI2010]弹飞绵羊 solution 这是一道LCT的裸题。 但是我并不想用LCT解决此题&#xff08;In fact 是不会LCT ~QAQ&#xff09; 于是我们开始大力分块。 考虑把弹跳装置分块&#xff0c;我们每次需要知道在一…

L - Two Ants Gym - 102823L

L - Two Ants Gym - 102823L 题意&#xff1a; 有两个线段A&#xff0c;B&#xff0c;两个线段不会超过一个公共点&#xff0c; 你站在线段B上&#xff0c;整个平面你看不到的区域的面积(如图中S所在区域) 题解&#xff1a; 计算几何&#xff0c;恶心题。调了一个小时还是…

C# 8.0 中开启默认接口实现

当你升级到 C# 8.0 和 .NET Core 3.0 之后&#xff0c;你就可以开始使用默认接口实现的功能了。从现在开始&#xff0c;你可以在接口里面添加一些默认实现的成员&#xff0c;避免在接口中添加成员导致大量对此接口的实现崩溃。要写出并且正常使用接口的默认实现&#xff0c;你需…

P4396 [AHOI2013]作业

P4396 [AHOI2013]作业 题目描述 详见&#xff1a;P4396 [AHOI2013]作业 solution 莫队树状数组的裸题&#xff08;莫队分块&#xff0c;CDQ分治都可以&#xff0c;莫队线段树大概需要卡常&#xff09;。 时间复杂度 Code #include<bits/stdc.h> using namespac…

Codeforces Round #732 (Div. 1Div. 2)

Codeforces Round #732 (Div. 2) 题号题目知识点AAquaMoon and Two Arrays贪心BAquaMoon and Stolen String异或CAquaMoon and Strange SortDAquaMoon and ChessEAquaMoon and PermutationsFAquaMoon and Wrong CoordinateG1AquaMoon and Time Stop (easy version)G2AquaMoon …

在编写异步方法时,使用 ConfigureAwait(false) 避免使用者死锁

我在 使用 Task.Wait()&#xff1f;立刻死锁&#xff08;deadlock&#xff09; 一文中站在类库使用者的角度看 async/await 代码的死锁问题&#xff1b;而本文将站在类库设计者的角度来看死锁问题。阅读本文&#xff0c;我们将知道如何编写类库代码&#xff0c;来尽可能避免类库…

P1903 [国家集训队]数颜色 / 维护队列

P1903 [国家集训队]数颜色 / 维护队列 题目描述 详见&#xff1a;P1903 [国家集训队]数颜色 / 维护队列 solution 裸的带修莫队。 多记录一维询问时已经进行多少次修改。 用莫队思想暴力即可。 &#xff08;不开O2时需要卡常&#xff0c;洛谷评测两发一样的代码&#xf…

妄想集合(牛客练习赛90)

妄想集合(牛客练习赛90) 题意&#xff1a; 开始有 n 个可重集合&#xff0c;开始时每一个集合中都有一个数&#xff0c;有 m 个操作。 Quant l r x\text{Quant l r x}Quant l r x&#xff1a;往编号在 l∼rl\sim rl∼r 的每个集合中加入一个数 x。 Ask l r\text{Ask l r}Ask …

韩国政府计划从Windows 7迁移到Linux

据《韩国先驱报》的报道&#xff0c;韩国内政部于上周四表示&#xff0c;韩国政府计划将其用于办公的计算机操作系统从 Windows 迁移至 Linux。当然不是马上全面迁移 —— 内政部将会“第一个吃螃蟹”。他们准备在其 PC 上试运行 Linux&#xff0c;如果没有出现安全问题&#x…

P3332 [ZJOI2013]K大数查询

P3332 [ZJOI2013]K大数查询 题目描述 详见&#xff1a;P3332 [ZJOI2013]K大数查询 solution 整体二分树状数组 code #include<bits/stdc.h> using namespace std; const int MAXN1e550; typedef long long ll; ll n,m,Ans[MAXN],id[MAXN],quel[MAXN],quer[MAXN],f[…

黑科技抢先尝(续2) - Windows terminal中Powershell Tab的美化全攻略

温馨提示: 原文中含有一些外部链接&#xff0c;点击全文左下角的"阅读原文"体验会更佳喔~ 接着之前的文章 黑科技抢先尝(续) - Windows terminal中WSL Linux 终端的极简美化指南&#xff0c;依然假定你安装好了windows terminal预编译版本。这次我的目标是将PowerShe…

CF455B A Lot of Games

CF455B A Lot of Games 题意&#xff1a; Andrew,Fedor和Alex是三个善于创造的人。现在&#xff0c;他们发明了一种字符串双人游戏。 给出n个非空字符串。在游戏中&#xff0c;两位玩家一起建造一个单词。开始时&#xff0c;单词是一个空字符串。两位玩家轮流操作。一位玩家…

P1712 [NOI2016]区间

P1712 [NOI2016]区间 题目描述 P1712 [NOI2016]区间 Solution 尺取法线段树 一个显然的想法是按区间长度排序。 每一次多选取一个区间相当于区间覆盖次数加1&#xff0c;每一次少选取一个区间就有区间覆盖次数减1。 可以用线段树维护区间覆盖次数的最大值。 于是转化成…