CUDA编程（十） 使用Kahan’s Summation Formula提高精度 上一次我们准备去并行一个矩阵乘法。然后我们在GPU上完毕了这个程序，当然是非常单纯的把任务分配给各个线程。也没有经过优化。终于我们看到，执行效率相当的低下，可是更重要的是出现了一个我们之前做整数立方和没遇到的问题，那 ...
分类：
其他好文 时间：
2017-07-19 20:29:24
阅读次数：
186
算(7+4*sqrt(3))^n的整数部分（mod 1e9+7）。 容易想到矩乘快速幂，但是怎么算整数部分呢？ (7+4*sqrt(3))^n一定可以写成a+b*sqrt(3)，同理(7-4*sqrt(3))^n一定可以写成a-b*sqrt(3)，于是， (7+4*sqrt(3))^n = (7+4 ...
分类：
其他好文 时间：
2017-07-18 22:14:14
阅读次数：
436
数论第一章：矩阵乘法。 包含知识点： 1. 矩阵乘法 2. 矩阵快速幂 3. 根据题目所给数据范围确认使用矩阵快速幂 4. 如何将一些问题转换成矩阵快速幂（难点） 5. 如何根据题目推导出矩阵的表达式（难点） 矩阵乘法的最基本思想：将n-1的矩阵推导到n的矩阵（做题的基本原则） 数论第二章：容斥原理 ...
分类：
其他好文 时间：
2017-07-16 20:14:34
阅读次数：
112
题目描述 这是一道模板题。 分别给定 n×p n \times pn×p 和 p×m p \times mp×m 的两个矩阵 A AA 和 B BB，求 A×B A \times BA×B。 这是一道模板题。 分别给定 n×p n \times pn×p 和 p×m p \times mp×m 的两 ...
分类：
其他好文 时间：
2017-07-11 19:12:11
阅读次数：
252
【BZOJ4861】[Beijing2017]魔法咒语 题意：别看BZ的题面了，去看LOJ的题面吧~ 题解：显然，数据范围明显的分成了两部分：一个是L很小，每个基本词汇长度未知；一个是L很大，每个基本词汇的长度是1或2。看来只能写两份代码了。 对于L很小的，我们先将禁忌串建成一个AC自动机，然后预处 ...
分类：
其他好文 时间：
2017-07-11 15:47:29
阅读次数：
225
P3390 【模板】矩阵快速幂 题目背景 矩阵快速幂 题目描述 给定n*n的矩阵A，求A^k 输入输出格式 输入格式： 第一行，n,k 第2至n+1行，每行n个数，第i+1行第j个数表示矩阵第i行第j列的元素 输出格式： 输出A^k 共n行，每行n个数，第i行第j个数表示矩阵第i行第j列的元素，每个 ...
分类：
其他好文 时间：
2017-07-11 01:04:50
阅读次数：
184
Fibonacci Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 15215 Accepted: 10687 Description In the Fibonacci integer sequence, F0 = 0, F1 = ...
分类：
其他好文 时间：
2017-07-11 01:03:44
阅读次数：
226
Numpy 数组的创建 数组的访问方法 蜜汁用法 布尔 数组计算 矩阵乘法&转置 求和 各种函数 http://link.zhihu.com/?target=http%3A//docs.scipy.org/doc/numpy/reference/routines.array-manipulation ...
分类：
编程语言 时间：
2017-07-10 00:32:24
阅读次数：
245
【BZOJ2510】弱题 Description 有M个球，一开始每个球均有一个初始标号，标号范围为1～N且为整数，标号为i的球有ai个，并保证Σai = M。 每次操作等概率取出一个球（即取出每个球的概率均为1/M），若这个球标号为k（k < N），则将它重新标号为k + 1；若这个球标号为N，则 ...
分类：
其他好文 时间：
2017-07-07 10:05:16
阅读次数：
153
题目描述 一个由自然数组成的数列按下式定义： 对于i <= k：ai = bi 对于i > k: ai = c1ai-1 + c2ai-2 + ... + ckai-k 其中bj和 cj （1<=j<=k）是给定的自然数。写一个程序，给定自然数m <= n, 计算am + am+1 + am+2 + ...
分类：
其他好文 时间：
2017-07-03 22:32:37
阅读次数：
206
