CF643F-Bears and Juice【组合数学】

正题

题目链接:https://www.luogu.com.cn/problem/CF643F


题目大意

题目有点奇怪就直接放翻译了

nnn 只熊和若干桶果汁和恰好一桶酒,每一天每只熊会选择一些桶(可能不选)并各喝一 杯,喝到酒的熊会去睡觉并不再回来,通过这个信息,熊们想知道哪个桶里是酒。

只有 ppp 个睡 觉的位置,当睡觉的熊超过了 ppp 只或者所有熊都在睡觉时熊们就失败了。

RiR_iRi 表示在 iii 天内桶的数量最多少,使得熊可以成功知道酒的位置。令 Xi=(i×Ri)mod232X_i = (i\times R_i) \bmod 2^{32}Xi=(i×Ri)mod232,你需要求出 X1⊕X2⊕…⊕XqX_1 \oplus X_2 \oplus\ldots \oplus X_qX1X2Xq

1≤n≤1091\leq n\leq 10^91n1091≤p≤1301\leq p\leq 1301p1301≤q≤2×1061\leq q \leq 2\times 10^61q2×106


解题思路

之前在XJ杂题选讲时候的神奇题目

题目比较乱但是我们发现题目问的是最多的数量,而不是最劣情况下的最多数量,所以这个东西是在最优情况下能分辨的数量。

这是我们之前很少接触的一种形式,这里需要用到信息的概念,因为我们是最优的,相当于我们所有的情况都可以去尝试,也就是每种信息都可以为我们选出一个答案,那么显然我们让选出的这些答案两两不同肯定就是最优的,所以这里的RiR_iRi就表示iii天以内我们能够获取的信息的数量

那么我们现在能够得到的信息数就是有多少头熊睡着了,和分别在哪一天睡着的,那么有
Ri=∑j=1min{p−1,n}(nj)ijR_i=\sum_{j=1}^{min\{p-1,n\}}\binom{n}{j}i^jRi=j=1min{p1,n}(jn)ij
也就是组合睡觉的熊,然后每个睡觉的都可以在任意天的时候睡觉

这个东西主要是(nj)\binom{n}{j}(jn)因为没有逆元比较麻烦,但是因为jjj比较小所以我们可以直接暴力枚举上下的因子然后消掉他们的gcdgcdgcd就好了

时间复杂度O(p3log⁡p+q×p)O(p^3\log p+q\times p)O(p3logp+q×p)


code

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
using namespace std;
int n,p,q;
unsigned ans,f[300];
vector<int> a,b;
int main()
{scanf("%d%d%d",&n,&p,&q);p=min(n-1,p);for(int i=0;i<=p;i++){a.clear();b.clear();f[i]=1;for(int j=0;j<i;j++)a.push_back(n-j);for(int j=1;j<=i;j++)b.push_back(j);for(int x=0;x<a.size();x++)for(int y=0;y<b.size();y++){int d=__gcd(a[x],b[y]);a[x]/=d;b[y]/=d;}for(int x=0;x<a.size();x++)f[i]=1u*a[x]*f[i];}for(int i=1,t=1;i<=q;i++){unsigned tmp=0,k=1;for(int j=0;j<=p;j++,k=1u*i*k)tmp+=f[j]*k;ans^=1u*i*tmp;}printf("%u",ans);
}

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

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

相关文章

【模板】Manacher算法

ACM模板 目录初始化以及构造初始化以及构造 数组下标从1开始即cin>>s1 由于奇回文串和偶回文串某些性质不同&#xff0c;我们首先通过init()操作使得新串中所有回文串的长度都变成奇数&#xff0c;返回值是新串的长度&#xff08;原串中的下标i对应新串中的2i&#xff09…

POJ2155 - Matrix(二维树状数组)

POJ2155 - Matrix 文章目录题目题解&#xff1a;代码题目 给你一个二维矩阵&#xff0c;初始化为0&#xff0c;然后可以进行两次操作&#xff1a; C:x,y,x1,y2 对该小矩阵内的数进行取反 Q:查询某个点是0还是1 题解&#xff1a; C是区间修改&#xff0c;Q是单点查询&#xf…

牛客挑战赛48C-铬合金之声【Prufer序列】

正题 题目链接:https://ac.nowcoder.com/acm/contest/11161/C 题目大意 nnn个点加mmm条边使得不存在环&#xff0c;每种方案的权值是所有联通块的大小乘积。 求所有方案的权值和。 1≤n≤109,1≤m≤1071\leq n\leq 10^9,1\leq m\leq 10^71≤n≤109,1≤m≤107 解题思路 就是分…

F - Sugoroku2(期望dp)

F - Sugoroku2 期望dp&#xff1a;fuf_ufu​表示到达终点的期望于是有 fu{0,u≥nf0,u∈Afu1⋯fumm1f_u\begin{cases} 0,u\ge n \\f_0,u\in A\\ \frac{f_{u1}\dotsf_{um}}{m}1 \end{cases}fu​⎩⎪⎨⎪⎧​0,u≥nf0​,u∈Amfu1​⋯fum​​1​ 由于存在返回f0f_0f0​的操作&…

.NET Core2使用Azure云上的Iot-Hub服务

基于工业4.0大背景下的工业物联网是近几年内热门的话题&#xff0c;依靠信息化技术企业可以实现数字化转型&#xff0c;生产可以实现智能化制造&#xff0c;设备可以实现自动化运作。然而&#xff0c;海量的数据采集是整个建设过程的基础环节&#xff0c;如何处理与利用这海量的…

Matrix Subtraction(小米icpc邀请赛第一场)

Matrix Subtraction 题意&#xff1a; 一个给定的矩阵&#xff0c;然后给定一个子矩阵的大小&#xff0c;子矩阵可以 将覆盖矩阵的区域的值减1&#xff0c;问能否将矩阵全部减为0 题解&#xff1a; 思路和下面这个链接讲的题十分相似 传送 本质就是二维树状数组差分求解 用…

E - Rotate and Flip(线性代数矩阵坐标变换)

E - Rotate and Flip 顺时针 [010−100001][xy1][y−x1]\begin{bmatrix} 0 & 1 & 0\\ -1&0&0\\0 &0&1 \end{bmatrix}\begin{bmatrix} x\\y\\1 \end{bmatrix}\begin{bmatrix} y\\-x\\1 \end{bmatrix}⎣⎡​0−10​100​001​⎦⎤​⎣⎡​xy1​⎦⎤​⎣…

CF755G-PolandBall and Many Other Balls【倍增FFT】

正题 题目链接:https://www.luogu.com.cn/problem/CF755G 题目大意 nnn个东西排成一排&#xff0c;每个组可以选择一个单独的物品或者两个连续的物品&#xff0c;一个物品不同同时在两个组里&#xff0c;但是可以不在组里。对于i∈[1,k]i\in[1,k]i∈[1,k]求分成iii组的方案数。…

Intelligent Warehouse(小米邀请赛)

题目 题意&#xff1a; n个数字&#xff0c;问存在的最长的一组数&#xff0c;使得其中任意两个数的都是倍数关系&#xff0c;问最长的长度是多少 题解&#xff1a; 暴力。。。 没想到暴力就能做&#xff0c;当时就该交上去试试的 用dp[i]表示当期选的所有数都是i的约数且符…

SmartCode.ETL 这不是先有鸡还是蛋的问题!

继国庆节 SmartCode 正式版&#xff08;SmartCode.Generator&#xff09;发布之后&#xff0c;SmartCode 迎来了新的能力 SmartCode.ETL !SmartCode 正式版从开始发布就从未说过自己仅仅是个代码生成器&#xff0c;这点上从我第一次宣布SmartCode正式开源的文章就可以说明&…

codeforces1469 E. A Bit Similar

E. A Bit Similar 代码抄的这份题解 找到一个长度为k的串使得在s长度为k的子串中&#xff0c;对于所有子串答案串都至少有一位相同。也就是如果将s长度为k的子串全部按位取反&#xff0c;那么答案串不能与这些反串匹配。 现在问题转化成&#xff0c;找到一个答案串&#xff…

对拍

参考博客 晚上不想刷题。。无聊的开始研究对拍 以前oi的时候看过对拍&#xff0c;奈何当时水平太差根本用不上对拍&#xff0c;&#xff08;现在水平也不咋地。。&#xff09; 上面那个博文讲的非常好。。 对拍最终的就是 对拍程序.bat 首先第一步&#xff1a;生成一组输入数据…

P7444-「EZEC-7」猜排列【dp】

正题 题目链接:https://www.luogu.com.cn/problem/P7444 题目大意 一个长度为nnn的排列&#xff0c;已知每个cic_ici​表示那个排列中mexmexmex为iii的区间个数。求满足条件的排列个数 1≤n≤5105,ci≥0,∑i1ncin(n1)2−11\leq n\leq 5\times 10^5,c_i\geq 0,\sum_{i1}^nc_i\…

dubbo线程池为什么耗尽

文章概述 大家可能都遇到过DUBBO线程池打满这个问题&#xff0c;报错如下&#xff0c;本文我们就一起分析DUBBO线程池打满这个问题。 cause: org.apache.dubbo.remoting.RemotingException: Server side(10.0.0.100,20881) thread pool is exhausted, detail msg:Thread pool …

手把手教你写DI_0_DI是什么?

DI是什么&#xff1f;Dependency Injection 常常简称为&#xff1a;DI。它是实现控制反转&#xff08;Inversion of Control – IoC&#xff09;的一个模式。fowler 大大大神 “几十年”前的经典文章 https://www.martinfowler.com/articles/injection.html 说的很清楚。“几十…

C++顺序创建txt文件

今天下午在对拍数据的时候感觉好麻烦。。 一次次手调文件名称&#xff0c;突然想起可以直接写段程序集中操作 顺序输出 #include <iostream> #include <fstream> #include<string> using namespace std;int main () {ofstream File;for(int i0;i<3;i){str…

【模板】字符串哈希

ACM模板 目录构建应用构建 字符串哈希就是将字符串映射成一个数&#xff0c;哈希冲突是不可避免的&#xff0c;我们需要选用合适的base尽可能使得哈希冲突可能性降低 unsigned long long溢出后相当于取模&#xff0c;相当于模264−12^{64}-1264−1 get(l,r)函数返回字符串下标…

P5137-polynomial【倍增】

正题 题目链接:https://www.luogu.com.cn/problem/P5137 题目大意 TTT组数据给出n,a,b,pn,a,b,pn,a,b,p求 (∑01naibn−i)%p\left(\sum_{01}^na^ib^{n-i}\right)\%p(01∑n​aibn−i)%p 1≤T≤105,1≤n,a,b,p≤10181\leq T\leq 10^5,1\leq n,a,b,p\leq 10^{18}1≤T≤105,1≤n,a…

使用Visual Studio Code开发.NET Core看这篇就够了

在本文中&#xff0c;我将带着大家一步一步的通过图文的形式来演示如何在Visual Studio Code中进行.NET Core程序的开发&#xff0c;测试以及调试。尽管Visual Studio Code的部分功能还达不到Visual Studio的水平&#xff0c;但它实际上已经足够强大来满足我们的日常开发。而且…

【模板】最大权闭合图

ACM模板 目录概念建图证明模板题概念 闭合图中所有的点的出边必须指向内部的点 建图 原图的边在网络流中的边容量是INF&#xff0c;如果点权是正&#xff0c;那么源点向其连边&#xff0c;容量是点权&#xff1b;否则它向汇点连边&#xff0c;容量是点权绝对值 证明 考虑最…