<题目链接> A为一个方阵,则Tr A表示A的迹(就是主对角线上各项的和),现要求Tr(A^k)%9973。 Input 数据的第一行是一个T,表示有T组数据。 每组数据的第一行有n(2 <= n <= 10)和k(2 <= k < 10^9)两个数据。接下来有n行,每行有n个数据,每个数据的范围是 ...
分类:
其他好文 时间:
2018-08-08 22:31:00
阅读次数:
200
一、快速幂 快速幂是用于解决类似$a^b$ $mod$ $p$值类型的问题的。使用普通的方法是从$1$循环至$b$,再逐次累乘,逐次取模。但这种方法对于$b$很大的时候却可能会超时。那么,这时候我们就需要使用快速幂了。 快速幂是基于以下式子: 若$b$ $mod$ $2=1$,则$a^b=a^\fr ...
分类:
其他好文 时间:
2018-08-04 22:33:08
阅读次数:
226
1 /* 2 题意:求a(n)=a(n-1)+a(n-2)+a(n-3)+......a(n-k)+... k 8 using namespace std; 9 10 typedef long long LL; 11 const int MAXN = 100005; 12 const LL MOD7... ...
分类:
其他好文 时间:
2018-07-31 23:25:18
阅读次数:
136
const int maxn = 105; struct Matrix { int m[maxn][maxn]; }ans,res; /**矩阵乘法**/ Matrix mul(Matrix a,Matrix b,int n){ Matrix tmp; for(int i=1;i<=n;i++){ ...
分类:
其他好文 时间:
2018-07-31 17:17:19
阅读次数:
150
P1962 斐波那契数列 题目背景 大家都知道,斐波那契数列是满足如下性质的一个数列: ? f(1) = 1 ? f(2) = 1 ? f(n) = f(n-1) + f(n-2) (n ≥ 2 且 n 为整数) 题目描述 请你求出 f(n) mod 1000000007 的值。 输入输出格式 输入 ...
分类:
其他好文 时间:
2018-07-29 20:04:21
阅读次数:
96
题目 分析 对着Trie图搞快速幂。 为什么这样是对的呢? 详见:http://www.matrix67.com/blog/archives/276 有些地方还不是很理解。。。。。。为什么节点还要往下扩展? 有空再来搞一搞 代码 1 //#include <bits/stdc++.h> 2 #inc ...
分类:
其他好文 时间:
2018-07-29 17:58:38
阅读次数:
150
题目链接 矩阵快速幂+二分图匹配, 对于对角线上的每个点看看能不能换到就行, 但是一开始$dicnic$写挂了 只好写的匈牙利 ...
分类:
其他好文 时间:
2018-07-27 01:23:53
阅读次数:
199
矩阵A乘矩阵B是A的第i行向量乘以B的第j列向量的值放在结果矩阵的i行j列。因为矩阵乘法满足结合律,所以它可以与一般的快速幂算法同理使用。注意矩阵在乘的时候取模。 ...
分类:
其他好文 时间:
2018-07-24 22:21:30
阅读次数:
193
题意 题目给出m(m<=10)个仅仅由A,T,C,G组成的单词(单词长度不超过10),然后给出一个整数n(n<=2000000000),问你用这四个字母组成一个长度为n的长文本,有多少种组成方法可以使得它不含任何一个给出的单词。 分析 当时一看以为是跟训练指南上(UVA11468)一样的题,感觉只有 ...
分类:
其他好文 时间:
2018-07-21 11:46:04
阅读次数:
169
//矩阵快速幂,本题二阶矩阵#include<cstdio>#include<cstring>#include<cstdlib>#include<iostream>#include<algorithm>using namespace std; const int mod = 1e9+7;struct ...
分类:
其他好文 时间:
2018-07-18 14:03:41
阅读次数:
204