LGOJ3747 六省联考2017 分手是祝愿
这两天遇到不少这种“人类智慧题”了,感觉都是很巧妙的
Description
现在有 \(n\) 盏灯,设每一次操作控制第 \(i\) 占灯,而改变状态的灯就是 \(i\) 的所有约数
现在给定初始的灯的状态序列,求剩余k次操作,就把灯全部关闭的步数期望\(+k\)和\(n!\) 的乘积
答案对 \(10003\) 取模
\(n \leq 10^5\)
Solution
思路分析
上来我们看到了“期望”,直接想到这题要 \(dp\)
然后定义状态是个难题(下面没有扯淡了)
\(f[i]\) 表示离全关掉还有 \(i\) 步走到离全关掉还有 \(i-1\) 步的期望操作次数。(这里是重点)
转移的时候考虑两种情况:
\(1^0\) 一次性摁对了,这种情况有\(\frac{i}{n}\)的概率
\(2^0\) 一次摁不对,需要转到\(i+1\)的状态
所以转移方程直接给出
\]
整理得:
\]
算法流程
最后给出本题流程:
1.\(O(n \sqrt n)\) 预处理因数的个数
2.从后往前扫一下,看一共需要几次完成游戏(特判如果\(cnt \leq k\),就直接乘上阶乘走人就好)
3.跑一下上面的 \(dp\),\(O(n)\)的,也不用优化
逆元啥的不会先去学板子吧
4.最后记得成阶乘
Code
#include <bits/stdc++.h>
using namespace std;
#define int long long
namespace yspm {
inline int read() {
int res = 0, f = 1;
char k;
while (!isdigit(k = getchar()))
if (k == '-')
f = -1;
while (isdigit(k)) res = res * 10 + k - '0', k = getchar();
return res * f;
}
const int N = 1e5 + 10;
vector<int> vec[N];
int n, k, f[N], mod = 100003, now[N], cnt, fac = 1, ans;
inline void prework() {
for (int i = 1; i <= n; ++i) {
for (int j = i; j <= n; j += i) vec[j].push_back(i);
}
for (int i = n; i >= 1; --i)
if (now[i]) {
++cnt;
int sz = vec[i].size();
for (int j = 0; j < sz; ++j) now[vec[i][j]] = !now[vec[i][j]];
}
return;
}
inline int ksm(int x, int y) {
int res = 1;
for (; y; y >>= 1) {
if (y & 1)
(res *= x) %= mod;
(x *= x) %= mod;
}
return res;
}
inline int inv(int x) { return ksm(x, mod - 2); }
signed main() {
n = read();
k = read(); f[n]=1;
for (int i = 1; i <= n; ++i) now[i] = read(), (fac *= i) %= mod;
prework();
if (cnt <= k)
return cout << cnt * fac % mod << endl, 0;
for (int i = n - 1; i > k; --i) f[i] = (n + (n - i) * f[i + 1] % mod) % mod * inv(i) % mod;
for (int i = cnt; i > k; --i) (ans += f[i]) %= mod;
cout << (ans+k) * fac % mod << endl;
return 0;
}
} // namespace yspm
signed main() { return yspm::main(); }
LGOJ3747 六省联考2017 分手是祝愿的更多相关文章
- BZOJ 4872 luogu P3750 [六省联考2017]分手是祝愿
4872: [Shoi2017]分手是祝愿 Time Limit: 20 Sec Memory Limit: 512 MB[Submit][Status][Discuss] Description ...
- bzoj千题计划266:bzoj4872: [六省联考2017]分手是祝愿
http://www.lydsy.com/JudgeOnline/problem.php?id=4872 一种最优解是 从大到小灯有亮的就灭掉 最优解是唯一的,且关灯的顺序没有影响 最优解 对每个开关 ...
- [BZOJ4872][六省联考2017]分手是祝愿(期望DP)
4872: [Shoi2017]分手是祝愿 Time Limit: 20 Sec Memory Limit: 512 MBSubmit: 516 Solved: 342[Submit][Statu ...
- [BZOJ4872][六省联考2017]分手是祝愿
BZOJ Luogu sol 首先发现肯定有解,又因为每个位置至多操作一次,所以最优解一定是在\([0,n]\)之间 有一种可以在\(O(\sum_{i=1}^{n}\lfloor\frac{n}{i ...
- luoguP3750 [六省联考2017]分手是祝愿 概率期望DP + 贪心
...........真的神状态了,没办法去想的状态................... 考试的时候选择$50$分贪心+$15$分状压吧,别的点就放弃算了........ 令$f[i]$表示从最小步 ...
- [六省联考2017]分手是祝愿 期望DP
表示每次看见期望的题就很懵逼... 但是这题感觉还是值得一做,有可借鉴之处 要是下面这段文字格式不一样的话(虽然好像的确不一样,我也不知道为什么,是直接从代码里面复制出来的,因为我一般都是习惯在代码里 ...
- P3750 [六省联考2017]分手是祝愿 期望DP
\(\color{#0066ff}{ 题目描述 }\) Zeit und Raum trennen dich und mich. 时空将你我分开. B 君在玩一个游戏,这个游戏由 \(n\) 个灯和 ...
- BZOJ4872 [六省联考2017]分手是祝愿 【期望dp】
题目 Zeit und Raum trennen dich und mich. 时空将你我分开.B 君在玩一个游戏,这个游戏由 n 个灯和 n 个开关组成,给定这 n 个灯的初始状态,下标为 从 1 ...
- 洛谷P3750 [六省联考2017]分手是祝愿(期望dp)
传送门 嗯……概率期望这东西太神了…… 先考虑一下最佳方案,肯定是从大到小亮的就灭(这个仔细想一想应该就能发现) 那么直接一遍枚举就能$O(nlogn)$把这个东西给搞出来 然后考虑期望dp,设$f[ ...
随机推荐
- jar类库加载顺序
当我们启动一个tomcat的服务的时候,jar包和claess文件加载顺序: 1. $java_home/lib 目录下的java核心api 2. $java_home/lib/ext 目录下的jav ...
- vmware搭建linux环境
参考博文:https://www.cnblogs.com/wuzm/p/10877990.html
- EUI库 - 快速入口之项目配置
egretProperties.json exmlRoot 指定exml文件存放根目录,该路径必须为相对路径,目录内只能有exml文件 themes 主题文件数组,配置所有主题文件路径,该 ...
- 15. react UI组件和容器组件的拆分 及 无状态组件
1.组件的拆分 组件拆分的前提 当所有的逻辑都出现在一个组件内时 组件会变得非常复杂 不便与代码的维护 所以对组件进行拆分 IU组件 进行页面渲染 容器组件 进行逻辑操作 UI组件的拆分 新建一个 ...
- Canvas基本定义
Android中使用图形处理引擎,2D部分是android SDK内部自己提供,3D部分是用Open GL ES 1.0.今天我们主要要了解的是2D相关的 大部分2D使用的api都在android.g ...
- spring boot 接口返回值封装
Spring Boot 集成教程 Spring Boot 介绍 Spring Boot 开发环境搭建(Eclipse) Spring Boot Hello World (restful接口)例子 sp ...
- day25(025-多线程(下)&GUI)
线程状态图 ###25.01_多线程(单例设计模式)(掌握) 单例设计模式:保证类在内存中只有一个对象. 如何保证类在内存中只有一个对象呢? (1)控制类的创建,不让其他类来创建本类的对象.priva ...
- 2020/2/1 PHP代码审计之任意文件读取及删除漏洞
在开始学习之前先简单记录一下自己现在的思路吧..现在接触的基本都是无防护的漏洞也就是最简单的一些漏洞.我的想法就是以代审思路为主,之前一直在打CTF,白盒的思维我觉得和CTF这种黑盒有很大区别.自己的 ...
- Spring创建Bean的顺序
一直对Spring创建bean的顺序很好奇,现在总算有时间写个代码测试一下.不想看过程的小伙伴可以直接看结论 目录结构: 其中:bean4.bean5包下的class没有注解@Component,测试 ...
- SpringCloud学习之Bus消息总线实现配置自动刷新(九)
前面两篇文章我们聊了Spring Cloud Config配置中心,当我们在更新github上面的配置以后,如果想要获取到最新的配置,需要手动刷新或者利用webhook的机制每次提交代码发送请求来刷新 ...