"传送门" 好久没写数论题了写一次调了1h 首先发现递推式是一个乘方的形式，线性递推和矩阵快速幂似乎都做不了，那么是否能够把乘方运算变成加法运算和乘法运算呢？ 使用原根！学过$NTT$的都知道$998244353$的原根$G=3$。 使用原根之后，可以得到一个等价的新递推式：$G^{g_i} = \ ...
分类：
其他好文 时间：
2019-02-01 11:29:08
阅读次数：
178
题目描述 加里敦星球的人们特别喜欢喝可乐。因而，他们的敌对星球研发出了一个可乐机器人，并且放在了加里敦星球的1号城市上。这个可乐机器人有三种行为： 停在原地，去下一个相邻的城市，自爆。它每一秒都会随机触发一种行为。现 在给加里敦星球城市图，在第0秒时可乐机器人在1号城市，问经过了t秒，可乐机器人的行 ...
分类：
其他好文 时间：
2019-01-29 23:16:41
阅读次数：
303
题目链接：http://poj.org/problem?id=3734 题意：给出n个排成一列的方块，用红、蓝、绿、黄四种颜色给它们染色，求染成红、绿的方块个数同时为偶数的方案数模10007的值。 题解：这是WC2019第二课堂生成函数的题，实际上可以不用生成函数，我们考虑如下状态：a[i]表示前i ...
分类：
其他好文 时间：
2019-01-27 16:32:35
阅读次数：
148
题面（权限题） 题解 一道概率$dp$，可以设$f[i][j]$表示第$i$次操作后，标号为$j$的小球的期望个数，那么有： $$ \begin{aligned} &f "i][j]=(1 \frac 1m)f[i 1][j]+\frac1mf[i 1][j 1" \\ &f[i][0]=(1 \f ...
分类：
其他好文 时间：
2019-01-26 17:11:10
阅读次数：
178
矩阵加速，专门用来解决一些递推的关系，其原理和矩阵运算的法则有关 由于矩阵的乘法有交换律和结合律，所以我们可以通过矩阵快速幂来快速求解递推关系，一般时间复杂度是O(nlogn)。 矩阵快速幂很简单，写一下模板就会了，但是推导单位矩阵是个难题。 一般地，我们推导单位矩阵时，有这几个步骤。 1. 确定递 ...
分类：
其他好文 时间：
2019-01-23 12:28:54
阅读次数：
182
No more tricks, Mr Nanguo Time Limit: 3000/1000 MS (Java/Others) Memory Limit: 65535/32768 K (Java/Others)Total Submission(s): 576 Accepted Submission ...
分类：
其他好文 时间：
2019-01-20 14:00:53
阅读次数：
160
struct Matrix { __int64 v[110][110]; Matrix() { memset(v, 0, sizeof(v)); } Matrix operator *(const Matrix B) // 重载的速度比写独立的函数慢点。 { int i, j, k; Matri..... ...
分类：
其他好文 时间：
2019-01-20 14:00:05
阅读次数：
199
题目描述 给定一个n*n的矩阵A以及一个正整数k，计算S=A1+A2+A3+?+Ak 输入 输入只有一组测试数据。输入的第一行包括三个正整数 n，k，m。接下来的 n 行每行包括 n 个非负整数(在[0,32768)之间)，按照行优先的顺序输入矩阵 A 的元素。 输出 输出 S 中每一个元素 mod ...
分类：
其他好文 时间：
2019-01-20 13:42:52
阅读次数：
181
引入问题：给定一个对角线非零的上三角矩阵$M$，求$M^k$，满足$M$的阶$\le 500$，$k\le 10^9$。 对998244353取模。 一个显而易见的算法是矩阵快速幂，然而是$O(N^3\log k)$的，无法通过本题。 ~~一开始我想，既然是上三角矩阵，那么特征多项式一定不难求，那么 ...
分类：
其他好文 时间：
2019-01-19 21:15:02
阅读次数：
139
题目地址 "题目链接" 题解 get一个新操作，序列的循环节 这据说这一种模数很小的序列的话都会有循环节？对这玩意有点萌币。 如果$Fib$序列中在$mod\space n$的意义下，出现了连续的0和1，就说明出现了循环节，对这个循环节进行取膜就行了 本来还打了矩阵快速幂的结果没用上 include ...
分类：
其他好文 时间：
2019-01-18 10:54:55
阅读次数：
244
