基础数学
小天位
这个作者很懒,什么都没留下…
展开
-
poj 2115 C Looooops(扩展欧几里德算法)
题意: 对于for(i=A ; i!=B ;i+=C)循环语句,问在k位存储系统中循环几次才会结束。 比如:当k=4时,存储的数 i 在0-15之间循环。(本题默认为无符号) 若在有限次内结束,则输出循环次数。 否则输出死循环。二,思路:本题利用扩展欧几里德算法求线性同余方程,设循环次数为 x ,则解方程 (A + C*x) % 2^k = B;求原创 2017-11-16 21:11:09 · 231 阅读 · 0 评论 -
hdu4549 M斐波那契数列
题意:F[0] = a ,F[1] = b,F[n] = F[n-1] * F[n-2] ( n > 1 )现在给出a, b, n,你能求出F[n]的值想到点子上就很简单了,可当时做的时候都没有向找递推式的方向去思考1.找F[n]的递推式附:斐波那契数矩阵公式 Fn+1 Fn 1 1 = ...原创 2018-10-29 20:42:38 · 188 阅读 · 0 评论 -
数论定理总结
0.威尔逊定理、判断一个数是不是素数得充分必要条件: p为素数,必有: (p-1)! ≡ -1 mod p 即:(p-1) !+ 1 ≡ 0 mod p 引证: 令 p=7 则有(7-1)!= 6! = 1 * 2 * 3 * 5 * 6 重排乘积中的因子并把乘积是模7的逆的分成一组 即:2*4 ≡ 1 mod 7 、3*5 ≡ 1 mod 7 ...原创 2018-10-29 22:20:29 · 358 阅读 · 0 评论 -
【light 1341Aladdin and the Flying Carpet】
给你矩形的面积(矩形的边长都是正整数),让你求最小的边大于等于b的矩形的个数。肯定是和面积的因子有关系了,边长都是面积值的 因子。那么就可以算出面积值得所有因子,因为是让求得矩形得个数,当然是要<<1了,还有个限制条件就是最小得边长要大于等于b得,那么把小于b得因子减去,得到的值就是符合条件的对数了210 212 2Sample OutputCase 1: 1C...原创 2018-12-01 13:01:27 · 152 阅读 · 0 评论 -
【2015沈阳现场A】
题意:求个一个n 求2^n-1能被7整除得有多少个,(1,2,....n)简单得同余方程#include<bits/stdc++.h>#include <iostream>#include <cmath>#include <cstdio>#include <stdlib.h>#include <ctime>...原创 2018-12-01 13:26:19 · 248 阅读 · 0 评论 -
【HDU1582 HDU1452 HDU1098 HDU3524 HDU1005 HDU2623 HDU2674】
题意:求2008^n的所有因子和对k取余,得到m,然后在求2008^m对给出的k取余。k是任意数(k>0)分析: 1. 用素因子唯一分解定理,对2008分解。 2. 素因子求和公式之后得到m 3. 因为要(a*b/250)%k 对于给定的k和250不确定是否互质,所以不能用求逆元的方法 4. 公式a/b%k=a%(b*k)/b 5. 快速幂求解最终结...原创 2018-12-02 09:08:52 · 473 阅读 · 0 评论 -
【HDU2683 TCE-frep number system 完全数+二项展开式】
g(n)是n的因子和A a b 查询a b区间有多少个n满足上式。Q a 查询a满不满足式子 参考#include <bits/stdc++.h>//#include <iostream>//#include <cstdio>#define X 10005#define inF 0x3f3f3f3f#define PI 3.141592653...原创 2018-12-02 08:28:26 · 299 阅读 · 0 评论 -
【HDU4497 GCD and LCM】
题意:已知l,g其中g=gcd(x,y,z),l=lcm(x,y,z),问有x,y,z多少种组合使得关系成立。 题解:已知x%g=y%g=z%g=0,l%x=l%y=l%z=0,所以l%g=0。这个可以判定是否存在(x,y,z)符合条件。gcd(a,b,c)=p; => gcd(a/p,b/p,c/p)=1 要不然gcd(a,b,c)>p且 lcm(a/p,b/...原创 2018-12-02 09:08:25 · 238 阅读 · 0 评论 -
【几种特殊的矩阵:对角矩阵、上下三角矩阵、正态分布随机矩阵、魔方矩阵、希尔伯特矩阵、托普利兹矩阵】
(这两个最不特殊了,线代中学过不过有点忘)1、对角矩阵: 对角矩阵(diagonal matrix)是一个主对角线之外的元素皆为0的矩阵,常写为diag(a1,a2,...,an) 。对角矩阵可以认为是矩阵中最简单的一种,值得一提的是:对角线上的元素可以为 0 或其他值,对角线上元素相等的对角矩阵称为数量矩阵;对角线上元素全为1的对角矩阵称为单位矩阵。对角矩阵的运算包括和、差运算、数乘...原创 2018-12-17 23:04:45 · 13111 阅读 · 2 评论 -
二进制相关
258. 各位相加给定一个非负整数num,反复将各个位上的数字相加,直到结果为一位数。10 111 212 313 414 515 617 818 9 019 120 221 322 423 524 6class Solution {public: int addDigits(int num) { return (num-1)%9+...原创 2019-03-25 18:15:48 · 827 阅读 · 0 评论 -
hdu1999 不可摸数 好题.
题意:s(n)是正整数n的真因子之和,即小于n且整除n的因子和.例如s(12)=1+2+3+4+6=16.如果任何数m,s(m)都不等于n,则称n为不可摸数. (2<=n<=1000)问一个数n是否是不可触摸数1.算出每个数作为其他数的因子时为因子和做出的贡献2.好题 #include <bits/stdc++.h>using namespace ...原创 2018-10-29 19:46:02 · 482 阅读 · 0 评论 -
lightoj-1028 Trailing Zeroes (I)
题意:给你一个数,问你有多少种进制对n的表示,存在后导零;1.挺简单的,画画一个数转化为二进制是怎么转化过来的,就可以看出来求出非1的因子个数有多少个就可以了用1进制表示没意义#include <bits/stdc++.h>#define inF 0x3f3f3f3f#define X 10005#define PI 3.1415926535897932384626...原创 2018-10-29 17:13:23 · 156 阅读 · 0 评论 -
poj1061
escription两只青蛙在网上相识了,它们聊得很开心,于是觉得很有必要见一面。它们很高兴地发现它们住在同一条纬度线上,于是它们约定各自朝西跳,直到碰面为止。可是它们出发之前忘记了一件很重要的事情,既没有问清楚对方的特征,也没有约定见面的具体位置。不过青蛙们都是很乐观的,它们觉得只要一直朝着某个方向跳下去,总能碰到对方的。但是除非这两只青蛙在同一时间跳到同一点上,不然是永远都不可能碰面的原创 2017-11-18 14:53:18 · 306 阅读 · 0 评论 -
poj 3233 待提交 Matrix Power Series
Matrix Power SeriesTime Limit: 3000MS Memory Limit: 131072KTotal Submissions: 16403 Accepted: 6980DescriptionGiven a n × n matrix A and a positive integ原创 2017-10-24 20:04:05 · 196 阅读 · 0 评论 -
HDU (1575)Tr A ---矩阵快速幂
Tr AProblem DescriptionA为一个方阵,则Tr A表示A的迹(就是主对角线上各项的和),现要求Tr(A^k)%9973。Input数据的第一行是一个T,表示有T组数据。每组数据的第一行有n(2 Output对应每组数据,输出Tr(A^k)%9973。Sample Input22 21 00 13 9原创 2017-10-23 20:10:42 · 227 阅读 · 0 评论 -
A Simple Math Problem HDU - 5974
Given two positive integers a and b,find suitable X and Y to meet the conditions: X+Y=a Least Comm原创 2017-10-14 18:02:01 · 286 阅读 · 0 评论 -
poj 3150 Cellular Automaton(迷糊,但原理是用的快速幂)
http://blog.csdn.net/guard_mine/article/details/44351103http://www.cnblogs.com/xin-hua/archive/2013/07/29/3223984.htmlhttp://www.cnblogs.com/Blue233333/p/7445402.html原创 2017-10-24 16:56:38 · 322 阅读 · 0 评论 -
hdu 1576
n=A%9973 A/B%9973=xA/B=x+9973kA=Bx+9973Bk则 n=Bx%9973 gcd(B,9973) = 1Bx=n+9973yBx+(-9973y)=n(不定方程)Bx≡n(mod 9973)(同余式)#include #define mod 9973using namespace std;int ex_gcd(int a,int b原创 2017-11-30 20:45:59 · 185 阅读 · 0 评论 -
poj 3070 Fibonacci
FibonacciTime Limit: 1000MS Memory Limit: 65536KTotal Submissions: 11123 Accepted: 7913DescriptionIn the Fibonacci integer sequence, F0 = 0, F1 = 1, and Fn = Fn − 1 + Fn − 2 for n ≥ 2. For example, th...原创 2017-10-24 15:42:19 · 217 阅读 · 0 评论 -
POJ-1006 Biorhythms (生物节律)
POJ1006 BiorhythmsTime Limit: 1000MS Memory Limit: 10000K Total Submissions: 139159 Accepted: 44664 DescriptionSome people believe that there are three cycles in a perso...原创 2017-09-19 20:32:02 · 971 阅读 · 2 评论 -
poj 2891
转载:/**********************一般模线性方程组***********************/同样是求这个东西。。X mod m1=r1X mod m2=r2.........X mod mn=rn首先,我们看两个式子的情况X mod m1=r1……………………………………………………………(1)X mod m2=r2…原创 2017-12-02 14:22:28 · 180 阅读 · 0 评论 -
cf552 G Minimum Possible LCMn个数,求最小得一对lcm
You are given an arrayaconsisting ofnintegersa1,a2,…,anYour problem is to find such pair of indicesi,j (1≤i<j≤n ) thatlcm(ai,aj) is minimum possible.lcm(x,y) is the least common multi...原创 2019-04-17 11:57:36 · 373 阅读 · 0 评论