2.23日刷数论题后总结(题目整理自SCUT
第一道:Rightmost digit
求N^N次最后一个数字
快速幂mod10咯
代码如下:
#include <cstdio>
#define ll long long
using namespace std;
const int mod = ; int qm(ll a,ll b) {
ll res = ;
while (b) {
if (b & ) res = res * a % mod;
a = a * a % mod; b >>= ;
}
return (int)res;
} int main() {
int T; scanf("%d", &T);
while (T--) {
int n;
scanf("%d", &n);
printf("%d\n", qm((ll)n, (ll)n));
}
return ;
}
好像可以打表找规律。
第二道:又见GCD
从2*b开始枚举c就可以啦
代码如下:
#include <cstdio>
using namespace std; int gcd(int a, int b) {
return b == ? a : gcd(b, a % b);
} int main() {
int n; scanf("%d", &n);
while (n--) {
int a, b;
scanf("%d%d", &a, &b);
for (int i = ; ; i++) {
int c = i * b;
if (gcd(a, c) == b) {
printf("%d\n", c);
break;
}
}
}
return ;
}
第三道:A / B
求 (A / B) % 9973
就是求A * INV(B) % 9973咯
由费马小定理 a^(p-1) ≡ 1(mod p)
所以 a * a^(p-2) ≡ 1(mod p)
所以INV(B) = b ^ (MOD-2)
代码如下:
#include <cstdio>
#define ll long long
using namespace std;
const int mod = ; int qm(ll a,ll b) {
ll res = ;
while (b) {
if (b & ) res = res * a % mod;
a = a * a % mod; b >>= ;
}
return (int)res;
} int main() {
int T; scanf("%d", &T);
while (T--) {
int a, b;
scanf("%d%d", &a, &b);
printf("%d\n", a * qm((ll)b, (ll)mod - ) % mod);
}
return ;
}
第四道:Fibonacci
最基本的矩阵快速幂(今天才学(可见刷题数量之少,以及之菜
和普通快速幂一个道理,就是加上矩阵相乘,以及定义ans变量的时候要定义成单位矩阵的形式
#include <cstdio>
using namespace std;
const int mod = ;
struct Matrix {
int m[][];
} ans, base; Matrix mul(Matrix a, Matrix b) {
Matrix temp;
for (int i = ; i < ; i++) {
for (int j = ; j < ; j++) {
temp.m[i][j] = ;
for (int k = ; k < ; k++) {
temp.m[i][j] = (temp.m[i][j] + a.m[i][k] * b.m[k][j]) % mod;
}
}
}
return temp;
} int qm(int n) {
base.m[][] = base.m[][] = base.m[][] = ;
base.m[][] = ;
ans.m[][] = ans.m[][] = ;
ans.m[][] = ans.m[][] = ;
while (n) {
if (n & ) ans = mul(ans, base);
base = mul(base, base);
n >>= ;
}
return ans.m[][];
} int main() {
int n;
while (~scanf("%d", &n) && n != -) {
printf("%d\n", qm(n));
}
return ;
}
扩展欧几里得算法
分为x % k == 0 和 不等于 0 两种情况
如果x % k == 0 直接找一个p + q = k的解打印就好了
如果不整除
令 a = [x / k] 即求 x = a * p + (a + 1) * q 一组整数解
又因为GCD(a, a + 1) = 1 所以就可以先用扩展欧几里得算法求a * p + (a + 1) * q = 1
p q再分别乘以x就好了
#include <cstdio>
#define ll long long
using namespace std; void extgcd(ll a, ll b, ll& d, ll& x, ll& y) {
if (b) {
extgcd(b, a % b, d, y, x);
y -= (a / b) * x;
} else {
d = a;
x = , y = ;
}
} int main() {
int T; scanf("%d", &T);
while (T--) {
ll x, k;
scanf("%lld%lld", &x, &k);
if (x % k == ) {
printf("%d %lld\n", , k);
continue;
}
ll p, q;
ll a = x / k;
ll d;
extgcd(a, a + , d, p, q);
p *= x / d, q *= x / d;
printf("%lld %lld\n", p, q);
}
return ;
}
第六道:zhx's contest
给定一个数n,求先增后减或先减后增或单调增或单调减的排列个数
先看 1,2,...,n
考虑把这个序列变成先增后减的情况 那么就有一个地方得容纳n
从左边n-1个数取0个放到n的右边 C(n-1, 0)
取1个 C(n-1, 1)
取两个 C(n-1, 2) 这两个元素取出来后放到n的右边只有一种排列方式, 就是从大到小排 因为要先增后减 到了n已经停止增 所以只能减
取3个 C(n-1, 3)
......
取n-1个 C(n-1, n-1)
那么可能数就为 C(n-1,0) + C(n-1,1) + ... + C(n-1, n-1) = (1 + 1) ^ (n - 1) = 2 ^ (n - 1)
那么对印的先减后增也是 2 ^ (n - 1)
但是会重复算单调增和单调减的序列
所以答案就是 2^n - 2
n, p <= 1e18
在快速幂的过程中很有可能会溢出long long
所以加上快速乘法(用加法模拟乘法,每一次都取模,这样保证不会溢出)
特判一下1的情况
代码如下:
#include <cstdio>
#define ll long long
using namespace std; ll quick_mul(ll a, ll b, ll mod) {
ll ans = ;
while (b) {
if (b & ) ans = (ans + a) % mod;
a = (a + a) % mod;
b >>= ;
}
return ans;
} ll quick_mod(ll a, ll b, ll mod) {
ll ans = ;
while (b) {
if (b & ) {
ans = quick_mul(ans, a, mod) % mod;
}
a = quick_mul(a, a, mod) % mod;
b >>= ;
}
return ans;
} int main() {
ll n, q;
while (~scanf("%lld%lld", &n, &q)) {
if (n == ) printf("%lld\n", % q);
else
printf("%lld\n", (quick_mod(, n, q) - + q) % q);
}
return ;
}
2.23日刷数论题后总结(题目整理自SCUT的更多相关文章
- 苹果越狱后必备软件,总有你需要的!11月23日追加14个,支持【iOS4】
http://bbs.dospy.com/thread-7398730-1-301-2.html越狱后必备软件,总有你需要的!11月23日追加14个,支持[iOS4] 背景自定义插件
- 北京Uber优步司机奖励政策(11月23日~11月29日)
用户组:人民优步"关羽组"(适用于11月23日-11月29日)奖励政策: 滴快车单单2.5倍,注册地址:http://www.udache.com/ 如何注册Uber司机(全国版最 ...
- 北京Uber优步司机奖励政策(12月23日)
滴快车单单2.5倍,注册地址:http://www.udache.com/ 如何注册Uber司机(全国版最新最详细注册流程)/月入2万/不用抢单:http://www.cnblogs.com/mfry ...
- Linux自用指令——2019年10月23日
1.ls ls命令是列出目录内容(List Directory Contents)的意思.运行它就是列出文件夹里的内容,可能是文件也可能是文件夹. ls -a 列出目录所有文件,包含以.开始的隐藏文件 ...
- Dataworks批量刷数优化方案探讨
Dataworks批量刷数优化方案探讨 在数据仓库的日常使用中,经常会有批量补数据,或者逻辑调整后批量重跑数据的场景. 批量刷数的实现方式,因调度工具差异而各有不同. Dataworks调度批量刷数局 ...
- [分享] 从定制Win7母盘到封装详细教程 By BILL ( 10月23日补充说明 )
[分享] 从定制Win7母盘到封装详细教程 By BILL ( 10月23日补充说明 ) billcheung 发表于 2011-10-23 00:07:49 https://www.itsk.com ...
- Week16(12月23日):复习
Part I:提问 =========================== 1.声明强类型视图时,使用关键字( ) A.ViewBag B.model C.Type D.Tit ...
- FJUT-这还是一道数论题
这还是一道数论题 TimeLimit:4000MS MemoryLimit:128MB 64-bit integer IO format:%lld Special Judge Problem D ...
- 北京Uber优步司机奖励政策(4月23日)
滴快车单单2.5倍,注册地址:http://www.udache.com/ 如何注册Uber司机(全国版最新最详细注册流程)/月入2万/不用抢单:http://www.cnblogs.com/mfry ...
随机推荐
- Mybatis学习总结(九)——查询缓存
一.什么是查询缓存 mybatis提供查询缓存,用于减轻数据压力,提高数据库性能.mybaits提供一级缓存和二级缓存. 1.一级缓存是sqlSession级别的缓存.在操作数据库时需要构造sqlSe ...
- Java调用FFmpeg进行视频处理及Builder设计模式的应用
1.FFmpeg是什么 FFmpeg(https://www.ffmpeg.org)是一套可以用来记录.转换数字音频.视频,并能将其转化为流的开源计算机程序.它用来干吗呢?视频采集.视频格式转化.视频 ...
- 【转】MySQL中的行级锁,表级锁,页级锁
在计算机科学中,锁是在执行多线程时用于强行限制资源访问的同步机制,即用于在并发控制中保证对互斥要求的满足. 在数据库的锁机制中介绍过,在DBMS中,可以按照锁的粒度把数据库锁分为行级锁(INNODB引 ...
- Python _内置函数3_45
reversed: #reversed() l = [1,2,3,4,5] l.reverse() print(l) #改变了原来的列表 l = [1,2,3,4,5] l2 = reversed(l ...
- Python_迭代器-生成器-复习-习题_41
# 迭代器和生成器# 迭代器 # 可迭代协议 —— 含有iter方法的都是可迭代的 # 迭代器协议 —— 含有next和iter的都是迭代器 # 特点 # 节省内存空间 # 方便逐个取值,一个迭代器只 ...
- vue的高阶组件
探索Vue高阶组件 探索Vue高阶组件的使用 Vue高阶组件的使用方法 高阶组件应用-组件重新实例化 深入理解React 高阶组件 探索Vue高阶组件 2018-01-05 探索Vue高阶组件 Vue ...
- 关于iframe页面里的重定向问题
最近公司做的一个功能,使用了iframe,父页面内嵌子页面,里面的坑还挺多的,上次其实就遇到过,只不过今天在此描述一下. 请允许我画个草图: 外层大圈是父级页面,里层是子级页面,我们是在父级引用子级页 ...
- 利用php查看某个服务的进程数
查看进程就是使用ps命令而已,只不顾ps的参数太多了. 使用php查询的话,必须要开启几个函数(可以执行外部程序的函数),参考官网:http://php.net/manual/zh/book.exec ...
- centos6 yum 安装memcached
centos6 yum 安装memcached - 像块石头 - 博客园http://www.cnblogs.com/rockee/archive/2012/08/01/2619160.html yu ...
- centos 7 aufs
Docker storage drivers | Docker Documentationhttps://docs.docker.com/storage/storagedriver/select-st ...