1.5.1 数论数学初步
Bryce1010_贤哉回也
github: https://github.com/Bryce1010
\n
e-mail: bryceyx@gmail.com
展开
-
UVA - 11582 Colossal Fibonacci Numbers!
/*题意:输入两个非负整数a、b和正整数n(0<=a,b<2^64,1<=n<=1000),你的任务是计算f(a^b)除以n的余数,其中f(0)=f(1)=1;且对于所有非负整数i,f(i+2)=f(i+1)+f(i)。思路:假设n=3,序列F[i]的前10项为1,1,2,0,2,2,1,0,1,1猜想这是有规律的,根据递推公式,第十一项会等于第3项,第1...原创 2018-07-13 16:11:46 · 172 阅读 · 0 评论 -
【唯一分解定理】LightOJ - 1236 Pairs Forming LCM
【唯一分解定理】LightOJ - 1236 Pairs Forming LCM题意:求1=<i<=n,1<=j<=n,lcm(i,j)=n的方案数。思路:n = p1 ^ e1 * p2 ^ e2 *…*pn ^ enfor i in range(1,n): ei 从0取到ei的所有组合必能包含所有n的因子。现在取n的两个因子a,ba=p1 ^ a1 * ...原创 2018-10-03 10:32:41 · 540 阅读 · 0 评论 -
【快速幂】LightOJ - 1213 Fantasy of a Summation
【快速幂】LightOJ - 1213 Fantasy of a Summationhttps://vjudge.net/problem/LightOJ-1213题意:计算k*n^(k-1).#include<bits/stdc++.h>using namespace std;#define ll long long ll MOD;ll powmod(ll a,ll n...原创 2018-10-03 10:29:57 · 217 阅读 · 0 评论 -
【数论】质因数分解+唯一分解定理
任何一个大于1的正整数都能唯一分解为有限个质数的乘积。试除法可以扫描2-sqrt(n)的每个数d,若d能整除N,则从N中除掉所有的因子d,同时累计除去d的个数。//质因数分解试除法void divide(int n){ m=0; for(int i=2;i&amp;amp;lt;=sqrt(n);i++) { if(n%i==0) { ...原创 2018-09-04 02:00:14 · 2447 阅读 · 0 评论 -
2018 Multi-University Training Contest 10----HDU - 6434 Problem I. Count
Bryce1010模板https://vjudge.net/problem/HDU-6434题意: 计算 ∑ni=1∑i−1j=1[gcd(i−j,i+j)=1]∑i=1n∑j=1i−1[gcd(i−j,i+j)=1] \sum{_{i=1}^n}\sum{_{j=1}^{i-1}}[gcd(i-j,i+j)=1]公式化简后就是计算1-n奇数项欧拉函数除以2,偶数项欧拉函数不变...原创 2018-08-24 14:05:11 · 387 阅读 · 0 评论 -
UVA - 10791 Minimum Sum LCM
Bryce1010模板 /*题意:求最小公倍数的最小和思路:唯一分解定理。*/#include&lt;bits/stdc++.h&gt;using namespace std;#define maxn 30000000int cnt=0;int v[maxn]={0};int a,b,c,n;int init(){ for (c=1;c&lt;=ma...原创 2018-07-13 16:25:39 · 240 阅读 · 0 评论 -
UVA - 12716 GCD XOR
/*题意:求最小公倍数的最小和思路:唯一分解定理。*/#include&lt;bits/stdc++.h&gt;using namespace std;#define maxn 30000000int cnt=0;int v[maxn]={0};int a,b,c,n;int init(){ for (c=1;c&lt;=maxn/2;c++) ...原创 2018-07-13 16:23:22 · 202 阅读 · 0 评论 -
UVA - 10375 Choose and divide
思路: 唯一分解定理。#include<bits/stdc++.h>using namespace std;const int MAXN=10000+10;int prime[MAXN];int e[MAXN];void getPrime(){ memset(prime,0,sizeof(prime)); for(int i=2;i<...原创 2018-07-13 16:21:09 · 162 阅读 · 0 评论 -
UVA - 12169 Disgruntled Judge
题意: 找3个整数x1,a,b;然后按照递推公式xi=(axi−1+b)mod10001xi=(axi−1+b)mod10001x_i=(a x_{i-1}+b )mod 10001计算出一个长度为2T的序列,其中T为测试数据的组数。然后把奇数项当做输入,偶数项当做输出。/*UVa12169 - Disgruntled Judge--a,b都在10000以内,这样可以枚举a的值,然后根...原创 2018-07-13 16:18:46 · 382 阅读 · 0 评论 -
【欧拉定理+思维】UVA - 11426 GCD - Extreme (II)
【欧拉定理+思维】UVA - 11426 GCD - Extreme (II)https://vjudge.net/problem/UVA-11426题意:计算∑i=1i&lt;N∑j=i+1j&lt;=NGCD(i,j)\sum{_{i=1}^{i&lt;N}} \sum{_{j=i+1}^{j&lt;=N}}GCD(i,j)∑i=1i<N∑j=i...原创 2018-10-05 20:00:57 · 358 阅读 · 0 评论