解题思路:分析需要不少时间,比较懒,直接把别人的分析贴在这里, 然后贴上自己写的代码:K相当之大。所以逐一递推的算法无法胜任。这时我们就不得不运用矩阵加速。首先来讲一下矩阵乘法:若一矩阵的列数与另一矩阵的行数相等,则可定义这两个矩阵的乘积。如 A 是 m×n 矩阵和 B 是 n×p矩阵,它们是乘积....
分类:
其他好文 时间:
2015-11-07 21:56:19
阅读次数:
239
矩阵乘法是可以分块的,而且幂的和也是具有线性的。不难得到Si = Si-1+A*Ai-1,Ai = A*Ai-1。然后矩阵快速幂就可以了。/********************************************************** ----------...
分类:
其他好文 时间:
2015-11-07 17:24:04
阅读次数:
252
vijosP1067Warcraft III 守望者的烦恼链接:https://vijos.org/p/1067【思路】 矩阵乘法。 可以得出递推式: f[i]=sum{ f[n-1],f[n-2]…f[n-k] } 矩阵乘法加速转移如下:1、 原始矩阵F 1 x k: | 1,0,0,0,0...
分类:
其他好文 时间:
2015-10-30 10:26:29
阅读次数:
172
1、线性代数矩阵乘法 dot 函数x= np.array([[1,2,3],[4,5,6]])y=np.array([[6,23],[-1,7],[8,9]])xOut[16]: array([[1, 2, 3], [4, 5, 6]])yOut[17]: array([[ 6, 23]...
分类:
编程语言 时间:
2015-10-27 22:07:05
阅读次数:
387
最近做作业,发现了一个一直不理解的问题终于明白了,高兴!block只是用来划分task。block和device memory 访问没有对应关系。 block可以访问任一device memory的内容。 之前,以为block只能访问划分给它的那块数据!!比如以前提过的两矩阵乘法, 矩阵小块可以除了...
分类:
其他好文 时间:
2015-10-23 06:43:56
阅读次数:
255
洛谷1349 广义斐波那契数列题目描述广义的斐波那契数列是指形如an=p*an-1+q*an-2的数列。今给定数列的两系数p和q,以及数列的最前两项a1和a2,另给出两个整数n和m,试求数列的第n项an除以m的余数。输入输出格式输入格式:输入包含一行6个整数。依次是p,q,a1,a2,n,m,其中在...
分类:
其他好文 时间:
2015-10-15 14:15:31
阅读次数:
215
头文件 结构的定义stdafx.h// stdafx.h : 标准系统包含文件的包含文件,// 或是经常使用但不常更改的// 特定于项目的包含文件//#pragma once#include "targetver.h"#include #include // TODO: 在此处引用程序需要的其他头....
分类:
编程语言 时间:
2015-10-13 20:43:36
阅读次数:
237
题意:给你一个N*N的矩阵,不用算矩阵乘法,但是每次询问一个子矩形的第K小数思路:整体二分+二维树状数组二分询问的答案mid,将数值小等mid的全部插入二维树状数组然后查询每个矩阵内的元素个数,若数量>K-1则放左边,否则放右边继续向下分治,左边二分l-mid,右边mid-r代码:#include#...
分类:
其他好文 时间:
2015-10-11 23:12:53
阅读次数:
415
在做矩阵乘法并行化测试的时候,在利用<time.h>的clock()计时时出现了一点问题。 首先看串行的程序: //?matrix_cpu.c
#include?<stdio.h>
#include?<stdlib.h>
#include?<time.h>
#define?NUM?2048
void?...
分类:
编程语言 时间:
2015-10-09 23:07:38
阅读次数:
537
1. 矩阵乘法(2015.09.28)答:向量点乘函数dot_product(A,B),矩阵乘法函数matmul(A,B),矩阵乘法包括向量点乘。例如,program matpro implicit none real :: a(3),b(3),c(3,1) a(1)=1.0 a(...
分类:
其他好文 时间:
2015-09-28 22:11:04
阅读次数:
237