OpenJudge百炼-2793-孙子问题-C语言-大衍求一术解一般同余方程组

最新推荐文章于 2022-03-08 19:51:37 发布
最新推荐文章于 2022-03-08 19:51:37 发布
阅读量2.1k 收藏 3
点赞数
分类专栏： C语言 算法 OpenJudge 文章标签： c语言 大衍求一术 数论 孙子问题 一般同余方程组
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/OrdinaryCrazy/article/details/77922767
版权
该博客探讨了《孙子算经》中的孙子问题及其解法，介绍了利用大衍求一术解决一般同余方程组的数论问题。博主分享了一段代码，但表示未能完全理解其巧妙之处，欢迎读者交流解题思路。
摘要由CSDN通过智能技术生成

描述：
我国古代《孙子算经》中，记有如下算题：“今有物不知其数，三三数之剩二，五五数之剩三，七七数之剩二，问物几何？”并给出得数：“答曰：23。”为解决这个问题民间流传了如下歌诀：“三人同行七十稀，五树梅花廿一枝，七子团员正半月，除百零五便得知。”把上面的问题说得明白一点就是：有一个正整数N，除以3的余数是2，除以5的余数是3，除以7的余数是2，要求这个数。民间给出的解法是：把N除以3的余数乘以70，把N除以5的余数乘以21，把N除以7的余数乘以15，把这三个结果加起来，最后把得到的结果除以105得到的就是答案。其实在民间的解法中不除以105得到的也是一个符合题意的答案，而且民间的解法对于已知“除以3的余数，除以5的余数和除以7的余数” 的问题都能得到一个符合要求的答案。比如对于上面的问题，得到的结果是2 * 70 + 3 * 21 + 2 * 15 = 233，这个结果也能满足除以3的余数是2，除以5的余数是3，除以7的余数是2。如果已知的问题是“除以3的余数是1，除以5的余数是4，除以7的余数是4”，民间解法得到的结果1 * 70 + 4 * 21 + 4 * 15 = 214，这个结果也满足除以3的余数是1，除以5的余数是4，除以7的余数是4。把这个问题推广到更普遍的情况：对于给定的正整数a1, a2, ... an，是否存在正整数b1, b2, ... bn，使得对于任意的一个正整数N，

最低0.47元/天 解锁文章
OrdinaryCrazy
关注
专栏目录
6取余11c语言,中国剩余定理“大衍求一术”手算方法及四个习题
weixin_32233571的博客
05-22 1656
中国剩余定理“大衍求一术”手算方法及四个习题2016-01-24佛山下雪今天广东佛山下雪。有一些雪花还真有点像鹅毛。这是我退休后移居广东20年来，第一次看到的景像。一中国剩余定理解题过程“中国剩余定理”是公元5-6世纪、我国南北朝时期的一部著名算术著作《孙子算经》中的一个“物不知数”的解法问题：今有物，不知其数，三三数之剩二，五五数之剩三，七七数之剩二。问物几何？答曰：二十三。术曰：「三三数之剩二...
【算法】大衍数算法
最新发布
Xiaoxin的博客
07-05 1025
大衍数算法，特别是“大衍求一术”和“大衍总数术”，是中国古代数学中的杰出成就，尤其体现在南宋数学家秦九韶的《数书九章》中。
一次同余方程与大衍求一术
dianshu1593的博客
02-27 2151
一次同余方程 前面已经提到，剩余类可以看作一个特殊的“数”，剩余类环可以看作定义了剩余类加法和乘法的“数集”.类似于实数集情形，我们也可以在剩余类环中解方程或方程组。 例如，在模6的剩余类环中解方程[5][x]=3，这里[x]是模6的剩余类环中的未知剩余类，注意到 $$[5][x] = [3]\Leftrightarrow [5x]=[3]\Leftrightarrow 6|5x-...
c语言孙子算法,孙子算法
weixin_39532466的博客
05-24 2110
《孙子算经》中的题目：有物不知其数，三个一数余二，五个一数余三，七个一数又余二，问该物总数几何？《孙子算经》中的解法：三三数之，取数七十，与余数二相乘；五五数之，取数二十一，与余数三相乘；七七数之，取数十五，与余数二相乘。将诸乘积相加，然后减去一百零五的倍数。解析：现在有三个数abc，满足a%3=2，b%5=3，c%7=2。如果1.b和c都是3的倍数，那么(a+b+c)%3依然对于2。同理2.a和...
三个线性同余方程组的计算机解决方案（C程序）
qq_37444385的博客
02-15 1309
线性同余方程组又称一次同余方程组，其形如ax≡b(mod n)。此方程有解当且仅当b能够被a与n的最大公约数整除，记作 gcd(a,n) | b。这时，如果 x0 是方程的一个解，那么所有的解可以表示为: {x0+kn/d|(k∈z)} 其中 d 是a 与 n 的最大公约数。在模 n 的完全剩余系 {0,1,…,n-1} 中，恰有 d 个解。 然而传统的数理解决方案对于数学基础较差的人来说不易理解。因此本文将提供一种基于C程序的计算工具，用于解决三个线性同余方程组。只需要输入解的范围，除数与余数即可。
POJ2115 C Looooops 解一元线性同余方程
AC_Gibson的专栏
06-05 1052
题目链接：http://poj.org/problem?id=2115 题目大意：有代码段for(variable=A;variable!=B;variable+=C) statement;已知所用到数都是mod 2^k的数，问你statement语句一共执行了多少次，如果循环为死循环，输出FOREVER. 分析：同余方程求解 相当于求A+Cx≡B（mod 2^
OpenJudge百炼-2747-数字方格-C语言-枚举
OrdinaryCrazy的博客
09-17 2922
描述： 如上图，有3个方格，每个方格里面都有一个整数a1，a2，a3。已知0 输入： 一行，包含一个整数n (0 输出： 一个整数，即a1 + a2 + a3的最大值。 样例输入： 3 样例输出： 5 /********************************************************** **文件名：百炼-2747 **Copyrig
OpenJudge百炼-2981-大整数加法-C语言-高精度计算
OrdinaryCrazy的博客
08-12 2284
描述： 求两个不超过200位的非负整数的和。 输入： 有两行，每行是一个不超过200位的非负整数，可能有多余的前导0。 输出： 一行，即相加后的结果。结果里不能有多余的前导0，即如果结果是342，那么就不能输出为0342。 样例输入： 22222222222222222222 33333333333333333333 样例输出： 55555555555555555555
OpenJudge百炼-2712-细菌繁殖-C语言-日期处理
OrdinaryCrazy的博客
08-10 2923
描述： 一种细菌的繁殖速度是每天成倍增长。例如：第一天有10个，第二天就变成20个，第三天变成40个，第四天变成80个，……。现在给出第一天的日期和细菌数目，要你写程序求出到某一天的时候，细菌的数目。 输入： 第一行有一个整数n，表示测试数据的数目。其后n行每行有5个整数，整数之间用一个空格隔开。第一个数表示第一天的月份，第二个数表示第一天的日期，第三个数表示第一天细菌的数目，第四个数表示要
OpenJudge百炼-2809-计算2的N次方-C语言-高精度计算
OrdinaryCrazy的博客
09-10 6584
描述： 任意给定一个正整数N(N 输入： 输入一个正整数N。 输出： 输出2的N次方的值。 样例输入： 5 样例输出： 32 /************************************************ **文件名：百炼-2809 **Copyright 2015-2025 (c) OrdinaryCrazy **创建人：OrdinaryCrazy *
中国剩余定理cpp代码
12-30
#include #include using namespace std; typedef int LL; typedef pair PLL; LL inv(LL t, LL p) {//求t关于p的逆元 if (t >= p) t = t%p; return t == 1 ? 1 : (p - p / t) * inv(p % t, p) % p; } LL gcd(LL a, LL b){ return b == 0 ? a : gcd(b, a%b); } PLL linear(LL A[], LL B[], LL M[], int n) {//求解A[i]x = B[i] (mod M[i]),总共n个线性方程组 LL x = 0, m = 1; for (int i = 0; i < n; i++) { LL a = A[i] * m, b = B[i] - A[i] * x, d =gcd(M[i], a); if (b % d != 0) return PLL(0, -1);//答案不存在，返回-1 LL t = b / d * inv(a / d, M[i] / d) % (M[i] / d); x = x + m*t; m *= M[i] / d; } x = (x % m + m) % m; return PLL(x, m);//返回的x就是答案，m是最后的lcm值 } int main() { int n; scanf_s("%d", &n); LL a[2017], b[2017], m[2017]; for (int i = 0; i<n; i++) { scanf_s("%d%d%d", &a[i], &b[i], &m[i]); } PLL pa = linear(a, b, m, n); printf("%lld\n", pa.first); } 设计思路： 有这样一道算术题：“今有物不知其数，三三数之剩二，五五数之剩三，七七数之剩二，问物几何？” 解这题，先构造一个答案 5*7*inv(5*7, 3) % 3 = 1 3*7*inv(3*7, 5) % 5 = 1 3*5*inv(3*5, 7) % 7 = 1 然后两边同乘你需要的数 2 * 5*7*inv(5*7, 3) % 3 = 2 3 * 3*7*inv(3*7, 5) % 5 = 3 2 * 3*5*inv(3*5, 7) % 7 = 2 令 a = 2 * 5*7*inv(5*7, 3) b = 3 * 3*7*inv(3*7, 5) c = 2 * 3*5*inv(3*5, 7) 那么 a % 3 = 2 b % 5 = 3 c % 7 = 2 其实答案就是a+b+c 因为 a%5 = a%7 = 0 因为a是5的倍数，也是7的倍数 b%3 = b%7 = 0 因为b是3的倍数，也是7的倍数 c%3 = c%5 = 0 因为c是3的倍数，也是5的倍数 所以 (a + b + c) % 3 = (a % 3) + (b % 3) + (c % 3) = 2 + 0 + 0 = 2 (a + b + c) % 5 = (a % 5) + (b % 5) + (c % 5) = 0 + 3 + 0 = 3 (a + b + c) % 7 = (a % 7) + (b % 7) + (c % 7) = 0 + 0 + 2 = 2 答案a+b+c完全满足题意 但是答案，不只一个，有无穷个，每相隔105就是一个答案（105 = 3 * 5 * 7） a=2*5*7*2=140 b=3*3*7*1=63 c=2*3*5*1=30 140+63+30=233 2335 = 23 如果题目问你最小的那个答案，那就是23了。 当 1*x=2(%3) 1*x=3(%5) 1*x=2(%7) 输入： 3 1 2 3 1 3 5 1 2 7 输出： 23
孙子问题
yinshi2009的专栏
11-10 1862
描述 我国古代《孙子算经》中，记有如下算题：“今有物不知其数，三三数之剩二，五五数之剩三，七七数之剩二，问物几何？”并给出得数：“答曰：23。” 为解决这个问题民间流传了如下歌诀：“三人同行七十稀，五树梅花廿一枝，七子团员正半月，除百零五便得知。” 把上面的问题说得明白一点就是：有一个正整数N，除以3的余数是2，除以5的余数是3，除以7的余数是2，要求这个数。 民间给出的解法是：把N除以
模素数幂的二次同余方程的程序C语言,二次同余方程的解
weixin_36281932的博客
05-24 533
今天要讨论的问题是解方程，其中是奇质数。 引理： 证明：由费马小定理， 引理：方程有解当且仅当 定理：设满足不是模的二次剩余，即无解，那么是二次 剩余方程的解。 证明：由，前面的等号用二项式定理和，后面的等 号用了费马小定理和是模的二次非剩余。然后 在算法实现的时候，对的选择可以随机，因为大约有一半数是模的二次非剩余，然后快速幂即可。 题意：求二次同余方程的解。代码：#include #incl...
POJ 2115 C Looooops(exgcd—解一元线性同余方程)
Chunzhen的博客
08-19 1171
C Looooops Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 23820   Accepted: 6591 Description A Compiler Mystery: We are given a C-language style for loop of
C语言编程求大衍数列,四川省三台中学实验学校2020学年高二数学上学期入学考试试题...
weixin_35715235的博客
05-19 412
四川省三台中学实验学校2020学年高二数学上学期入学考试试题Ⅰ卷(选择题)一、选择题(每小题5分,共60分)1．已知直线的倾斜角为，则为A． B． C． D． 不存在2. 设，，R，且＞，则A． B． C． D．3. 设，向量且 ，则A． B. C. D.4. 在△ABC中内角A，B，C所对各边分别为，，，...
大衍数列-基础算法（python / c）
qq_45024553的博客
03-08 3705
大衍数列一级目录c 一级目录 """ 中国古代文献中，曾记载过“大衍数列”, 主要用于解释中国传统文化中的太极衍生原理。 它的前几项是：0、2、4、8、12、18、24、32、40、50 … 其规律是：对偶数项，是序号平方再除2，奇数项，是序号平方减1再除2。 这是前一百位 """ for i in range(1, 101): if i % 2 == 0: print(int(i**2/2), " ", end="|") elif i % 2 != 0:
大衍求一术的计算方法
banqu3475的博客
08-04 6156
引言 大衍求一术是秦九韶发明的一种求特殊一次同余式的方法，以下摘自百度百科 秦九韶，字道古，生活于南宋时期，自幼喜好数学，经过长期积累和苦心钻研，于公元1247年写成《数书九章》。这部中世纪的数学杰作，在许多方面都有创造，其中求解一次同余组的“大衍求一术”和求高次方程数值解的“正负开方术”，更是具有世界意义的成就。 进入正题 对于形如ax ≡ 1 (mod n)，(a,...
蓝桥杯/第五届/大衍数列
weixin_30634661的博客
03-25 624
【问题描述】 中国古代文献中，曾记载过“大衍数列”, 主要用于解释中国传统文化中的太极衍生原理。 它的前几项是：0、2、4、8、12、18、24、32、40、50 ... 其规律是：对偶数项，是序号平方再除2，奇数项，是序号平方减1再除2。 以下的代码打印出了大衍数列的前 100 项。 for(int i=1; i<100; ...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值