- 博客(25)
- 资源 (6)
- 收藏
- 关注
RSS订阅原创 CodeForces 14C Four Segments
题意：给出四条线段的起点，终点坐标，四条线段平行或垂直坐标轴，问这四条线段能否构成矩形 矩形的特征就两点，对边等长平行，有四条边，抓住这几点就可以了#include #include #include using namespace std;struct node{ __int64 x1,y1,x2,y2,len;}a[4];struct node1{ __in
2015-05-31 21:41:30
813
原创 ZOJ 3547 The Boss on Mars（容斥）
题意：求在1~n中与n互质的数的四次方的和 容斥不会的建议先看一下入门题 解题思路：先素数打表，求出n的质因子，求出1~n的n个数的四次方的和，容斥原理，减去奇数个质因子相乘的数以及它的倍数的四次方，加上偶数个相乘的数以及它的倍数的四次方1~n的四次方求和公式：sum = n*(n+1)*(6*n*n*n+9*n*n+n-1)/30注意：2^4+4^4+6^4+8^
2015-05-26 23:53:53
537
原创 HDU 1588 Gauss Fibonacci（矩阵 + 斐波那契）
题意：已知g( i ) = k*i+b,f( i )是斐波那契数列输入k，b，n，mod，当0 解题思路：直接一项一项求，那肯定得超时啊，凭直觉，肯定得用矩阵，可关键就是如何去构造显然，g( i )是一个等差数列，那f( g(i) )也肯定是有某种规律的，我们可以大胆假设它是等差或等比，如果想跟矩阵结合，那么等差的可能性就不大，那就可能是等比拿第一组数据凑一下，k=2
2015-05-26 23:05:09
686
2
原创 hdu 1568 Fibonacci
题意：输出第n个 Fibonacci数的前四位 解题思路： Fibonacci数列的公式 用log做，举个例子，log10（1007） = 3+log10（1.007）,那么令a = log10（1.007）= 0.0030294710^a = 1.007那么可先求b = 然后求出b的小数部分，取出10^b的前四位即可#include #include
2015-05-25 21:01:52
569
原创 标准Fibonacci序列（简单介绍）
1、定义：F（0）=0，F（1）=1，F（n）=F(n-1)+F(n-2)（n≥2，n∈N*）（0, 1, 1, 2, 3, 5, 8, 13, 21, 34, 55, 89, 144, 233，377，610，987，1597，2584，4181，6765，10946，17711，28657，46368……）特别指出：第0项是0，第1项是第一个1。 2、通项公式： 3、特
2015-05-20 22:02:51
1038
原创 poj 3070 Fibonacci（标准Fibonacci序列 入门 水）
结合矩阵来写#include #include #include #include using namespace std;#define mod 10000struct mat{ int a[2][2];};mat mul(mat a,mat b){ mat ret; memset(ret.a,0,sizeof(ret.a)); for(int i=0;i<2
2015-05-20 21:58:50
502
原创 fzu 1683 纪念SlingShot （矩阵快速幂）
解题思路：求和需要小小的处理一下构造矩阵，一开始定义sum=f[0]+f[1]sum 1 1 0 0 sumf[3] 0 3 2 7 f[2]f[2] = 0 1 0 0 * f[1]f[1]
2015-05-19 23:36:04
717
原创 hdu 2604 Queuing（矩阵快速幂）
题意：一个字符串长为L，由f和m组成，字符串包含fff和fmf的叫做O队列，其余叫E队列，求有多少个E队列（输出对m取余后的值） 解题思路：字符串编号1~L，第L位有两种可能，m和f，用F[ i ]表示到第i位时的种类假设第L位是m，那么此时字符串的种类有F[L-1]假设第L位是f，那么最后三位有四种可能：mmf, mff, fmf, fff，其中fff和fmf不符合要求，
2015-05-19 20:48:01
473
原创 hdu 2256 Problem of Precision(矩阵快速幂 推公式)
题意：给出n，求（（sqrt（2）+sqrt（3））^（2*n））%1024 解题思路：（sqrt（2）+sqrt（3））^2 = （5+2*sqrt（6））(5+2*sqrt(6))^n = Xn+Yn*sqrt( 6 ) = ( Xn-1+Yn-1*sqrt( 6 ) ) * ( 5+2*sqrt( 6 ) )
2015-05-18 20:17:29
636
1
原创 FZU 1896 神奇的魔法数（ 组队题 ）
解题思路：用dp记录，一位以为算过去，一个数最高只有12位，不用担心超时问题#include #include #include #include #include using namespace std;#define LL __int64LL dp[15][15];int ai[15];void init(int l,int m){ for(int i=0;
2015-05-17 14:19:51
569
原创 FZU 1895 整除45问题
先考虑末位，要能被45整除，那么肯定能被5和9整除，被5整除，末位定是0和5，被9整除，所有数字之和是9的倍数抓住这两点，DFS一边即可，注意剪枝 #include #include #include using namespace std;const int N = 1005;const int M = 15;const int dir[5]= {0,5};int
2015-05-16 23:15:52
758
原创 HDU 5072 Hackthon （容斥 同色三角形变形）
题意：给出几个数字，求选出的数全部互质或全不互质的种数 解题思路：同色三角形原理，一个平面上有n个点，两个点之间可能是红线，也可能是连黑线，问一共能连出多少个同色三角形？从一个点出发，有n-1条线，假设红线有的d[ i ]条，黑线有n-1-d[ i ]条，那么从这一点出发不同色的三角形个数为d[ i ]*(n-1-d[ i ])个，所以总的不同色三角形个数为sum =d[
2015-05-16 22:24:16
568
原创 hdu 1757 A Simple Math Problem（矩阵快速幂 水~）
题意：If x If x >= 10 f(x) = a0 * f(x-1) + a1 * f(x-2) + a2 * f(x-3) + …… + a9 * f(x-10);输入k，m以及a0,a1,a2.....a9，计算f（k）%m 解题思路：构造矩阵矩阵a
2015-05-15 22:07:45
434
原创 HDU 1575 Tr A（矩阵快速幂 水 入门）
题意：求m个n*n的矩阵相乘后的矩阵的对角线之和 简单的矩阵相乘，可当入门题，套套模板就可以了 #include #include #include #include using namespace std;#define mod 9973int n;struct matrix{ int a[15][15];};matrix mul(matrix a,ma
2015-05-14 22:16:41
411
原创 POJ 3695 Rectangles（容斥）
题意：n个矩阵m次访问，求每次访问的总面积 参考http://http://blog.csdn.net/sdj222555/article/details/7983902#include #include #include #include #include using namespace std;#define MAXN 111111#define inf 1000
2015-05-14 21:56:20
389
原创 HDU 1796 How many integers can you find（容斥）
题意：输入n，m以及m个数，输出在1~n中能被m个数的其中一个整除的数的数量 解题思路：用一个数组记录这m个数，以及其中的2个，3个。。。m个数的最小公倍数能被奇数个数的最小公倍数整除的用+，偶数个用减 #include using namespace std;#define ll __int64ll a[2000],b[15];int gcd(int a,in
2015-05-13 19:44:36
432
原创 HDU 1018 Big Number（斯特林公式）
题意：给出一个数，输出这个数的阶乘的位数 Stirling公式: n! = ((2*pi*n)^(1/2))*((n/e)^n); 前提是n > 3；n！的位数就是以10为底的对数+1；那么这个式子就能转化为lg(n!) = (lg(2*pi)+lg(n))/2 + n*(lg(n)-lg(e)); #include #include using namesp
2015-05-12 21:57:47
823
原创 CodeForces 97B Superset （分治）
题意：给出n对坐标，可以在其中添加一些点，使得平面上任意两点满足以下三个条件的其中一个1、两个点在同一水平面上2、两个点在同一竖线上3、这两个点所围成的矩形里面有点存在最后输出满足条件的点的总数，并输出这些点 先排序，再解题 #include #include #include #include #include #include using nam
2015-05-11 20:10:36
952
原创 ZOJ 2573 DNA Problem（矩阵快速幂）
题意：基因题意：N个基因 经过M次突变 A代表S1,B代表S2，比如第一组数据，它的突变顺序就是S1,S2,S2,S1,S2,S2....我们要输出经过M次变异后的基因接下来不是有两个N*N的东西吗，第一个就是代表S1，第二个代表S2,举个例子0 1 00 1 10 1 0（横着是i，竖着是j）当i=0时，如何判断第i个基因是alive还是dead，就
2015-05-11 00:37:45
490
原创 hdu 4965 Fast Matrix Calculation（矩阵快速幂）
题意：矩阵A是n*k，矩阵B是k*n，矩阵C = A*B，先求C^(n*n)，再对矩阵C的每个元素对6取余，输出矩阵C的元素和 解题思路：直接按照题目顺序来求会超时 （ A*B ）^( n*n ) = A*((B*A)^( n*n -1) )*B先算P = (B*A)^( n*n -1)再算A*P*B这样转化的好处是在用矩阵快速幂计算时可以转化为6*
2015-05-09 23:54:07
370
原创 hdu 2276 Kiki & Little Kiki 2 (矩阵快速幂)
题意：一个时间t，一串由0和1组成的字符串s，长为n，编号为1~n，1代表灯是开的，0代表灯是关的，每次交换时，如果它的左边是1，那么它本身的状态改变，否则不变，1号灯的左边是n号灯，求t秒以后的状态（每一秒整个字符串的状态改变一次） 解题思路：首先就是要先发现规律第i个状态 = （ 前一个状态+它本身状态 ）%2发现这个一切好办，如果长为5，那么矩阵相乘 前一秒的
2015-05-09 16:31:40
502
原创 矩阵快速幂的写法（模板）
首先一般会定义一个结构体struct matrix{ int a[n][m];}; 2*2的矩阵直接可以用两个formatrix mul(matrix a, matrix b){ matrix ret; for(int i=0;i<2;i++) { for(int j=0;j<2;j++) { ret.a[i][j] = (a.a[i][0]*
2015-05-08 20:49:57
678
原创 CodeForces 300C Beautiful Numbers（排列组合 乘法逆元）
题意：给出a，b，n三个数一串由a，b组成的长度为n的数，各个数位相加和为sum，若sum也是由a，b组成的，那么ans++求最终的ans值 思路：0~n枚举过去，暴力若此时的sum满足条件，那么那就算出这个条件下的组合数，加到ans上去c( n,m ) = f[ n ]/(f[ m ]*f[ n-m ])乘法逆元：当mod = prime时，则（a/b）%
2015-05-04 22:12:50
1151
原创 HDU 4135 Co-prime (容斥)
题意：t组数据，每组输入a，b，n在a - b中找出有几个与n互质的 思路：找出n的素数因子，假设它的素数因子a1，a2，a3 用ans1记录1 - （ a-1 ）中与n不互质的个数（ans1 = （a-1）-a1的倍数-a2的倍数-a3的倍数+a1*a2的倍数+a1*a3的倍数+a2*a3的倍数-a1*a2*a3的倍数） ==> 不难发现奇数个用减，偶数个用加
2015-05-04 18:53:13
488
原创 CodeForces 296B Yaroslav and Two Strings (容斥)
题意：两串长为n的字符串，若满足s[ i ] >w[ i ] && s[ j ] 给出的字符串是由0-9以及？组成的，？是未知的意思，可代替0-9任意一个数，问一共有多少个non-comparable 解题思路：数字的处理很简单，两者比一下就可以了，主要是？难处理，提及容斥原理，那么不难想到用ans记录所有的情况，那么每遇到一个？就乘以10（不考虑大小，0-9十个数均可
2015-05-03 23:37:07
682
sourcetree中文版 Windows
2018-02-28
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人