- 在线时间
- 69 小时
- 最后登录
- 2014-11-17
- 注册时间
- 2012-3-15
- 听众数
- 4
- 收听数
- 0
- 能力
- 0 分
- 体力
- 1378 点
- 威望
- 0 点
- 阅读权限
- 40
- 积分
- 505
- 相册
- 0
- 日志
- 5
- 记录
- 6
- 帖子
- 174
- 主题
- 26
- 精华
- 0
- 分享
- 2
- 好友
- 28
升级   68.33% TA的每日心情 | 郁闷 2014-11-17 00:20 |
|---|
签到天数: 106 天 [LV.6]常住居民II
 群组: Matlab讨论组 群组: 学术交流A 群组: 学术交流B 群组: 数学建模培训课堂1 群组: 哈尔滨工业大学建模团 |
百鸡问题* D" f" F! J8 ^- g. H
5 X+ p+ L7 x+ G5 G 今有鸡翁一,值钱伍;鸡母一,值钱三;鸡鶵三,值钱一。凡百钱买鸡百只,问鸡翁、母、鶵各几何?答曰:鸡翁四,值钱二十;鸡母十八,值钱五十四;鸡鶵七十八,值钱二十六。又答:鸡翁八,值钱四十;鸡 母十一,值钱三十三,鸡鶵八十一,值钱二十七。又答:鸡翁十二,值钱六十;鸡母四、值钱十二;鸡鶵八十 四,值钱二十八。” F( W1 U0 ^- g1 ~: j8 R
编辑本段
% Z, X; R0 k& b原书说明
4 o! p; P1 W4 O M& a7 g# {$ t" R
0 }9 v6 t- d: k6 O8 O E/ {. ^ `% o 原书没有给出解法,只说如果少买7只母鸡,就可多买4只公鸡和3只小鸡。所以只要得出一组答案,就可以推出其余两组答案。中国古算书的著名校勘者甄鸾和李淳风注释该书时都没给出解法,只有约6世纪的算学家谢察微记述过一种不甚正确的解法。到了清代,研究百鸡术的人渐多,1815年骆腾风使用大衍求一术解决了百鸡问题。1874年丁取忠创用一个简易的算术解法。在此前后时曰醇(约1870)推广了百鸡问作《百鸡术衍》,从此百鸡问题和百鸡术才广为人知。百鸡问题还有多种表达形式,如百僧吃百馒,百钱买百禽等。宋代杨辉算书内有类似问题,中古时近东各国也有相仿问题流传。例如印度算书和阿拉伯学者艾布·卡米勒的著作内都有百钱买百禽的问题,且与《张邱建算经》的题目几乎全同。# J$ I& R) O6 e' L j$ |* B) X5 g
编辑本段' B9 ^4 o' @* B( `, q7 c
解法7 n4 N$ y& y4 Q5 _ k8 k
8 J5 h% O: C( a; E0 h 中国古代算书《张丘建算经》中有一道著名的百鸡问题:公鸡每只值5 文钱,母鸡每只值3 文钱,而3 只小鸡值1 文钱。现在用100 文钱买100 只鸡,问:这100 只鸡中,公鸡、母鸡和小鸡各有多少只?& L0 ^4 L+ N2 [; A
这个问题流传很广,解法很多,但从现代数学观点来看,实际上是一个求不定方程整数解的问题。解法如下:
, |6 g2 M$ @! ^# P; g$ A9 j- u1 V 设公鸡、母鸡、小鸡分别为x、y、z 只,由题意得:
2 o2 l" G' w% q ①……x+y+z =100
8 V' n9 ]0 f w ②……5x+3y+(1/3)z =1002 U! e6 S4 A0 C4 U' Z
有两个方程,三个未知量,称为不定方程组,有多种解。
6 ]5 ]5 d3 T# M2 g 令②×3-①得:7x+4y=100;
, D' S! n5 R5 w) a 所以y=(100-7x)/4=25-2x+x/4' u& h& H/ d5 [3 h
令x/4=t, (t为整数)所以x=4t
( _; D" `3 w. O, X& w4 W 把x=4t代入7x+4y=100得到:y=25-7t$ n* t$ S6 ~1 q1 @' | {/ K5 W
易得z=75+3t8 i/ {! D0 M0 G) n
所以:x=4t* \( @3 [0 `$ D; T" H+ ]
y=25-7t' r/ l- N! ~8 @. ?. A, [6 d- o
z=75+3t% w6 i* O1 W. j4 F8 P7 Y
因为x,y,z大于等于0; F% p! @ }( S! Y& G( E
所以4t大于等于0
; Z9 l/ K) R2 ^3 m1 S2 f 25-7t大于等于0
7 q. z! J4 K }1 o3 X! a( z% D a5 `# r 75+3t大于等于0
& _: [2 O b$ O+ E8 T 解得t大于等于0小于等于25/7 又因为t为整数0 `5 B% K& J [/ w
所以t=0,1,2,3(这里不要忘记t有等于0得可能)
8 x$ r$ E) l" Y" o5 h 当t=0时6 n& `4 M" h1 ?1 t Y
x=0,y=25,z=75
' Z4 R% {# ~' b1 l2 C8 x 当t=1时
/ g/ H) u' f; \ x =4;y =18;z =78
% h( P8 i7 _1 e: f 当t=2时; R6 H/ n$ q1 q S% o) U7 K
x =8;y =11;z =81
/ j2 e' Y* B8 L) f8 y 当t=3时6 I( X* F/ F' C- z
x =12;y =4;z =84; b6 T, O/ B3 N6 a. Q
编辑本段
( q/ X) d9 J2 R# @C语言解法% J q7 L; G# j, g
0 R4 ]' c, \ X; a% v; m
8 d1 u3 S0 O3 T* z/ Y#include <stdio.h>
1 }4 j4 g+ w! O; G) e; wvoid main(); g! X, x8 r( Z( Q$ P/ q0 B) [
{
/ b) [% k5 U1 z( }& j3 ?int cocks=0,hens,chicks;, {( }/ ]1 d" Z3 q; t
while(cocks<=20)
. b; {" [, l6 N3 y! E/ I{ @9 |- d8 J8 d$ L7 y; i
hens=0;. I2 B! V$ ~4 B) q% D5 {8 I1 c4 q
while(hens<=33)7 y3 F& X4 a1 A6 o# \9 o1 b
{
3 _9 K) T5 V' h9 x9 Q7 r! O4 g. Ichicks=100-cocks-hens;' p" h4 d5 r# O9 _' F( @
if(5.0*cocks+3.0*hens+chicks/3.0==100.0): R% h2 A; @* d( M1 g
printf("公鸡%d只,母鸡%d只,小鸡%d只\n\n",cocks,hens,chicks);: f) J1 _ ~# L
hens++;
4 C/ n( C- ~" w' f$ X. J. {. B8 \}7 A$ _5 q. T6 Z8 C
cocks++;
# U, W) L' M" j0 u( C2 Y+ d `}
9 d6 \# S5 x7 r# C. m; }- G R- B}
# _6 v* ]' p! u4 r% y2 a2 F输出结果为:% G- Z% d T" v& L2 p
公鸡0只,母鸡25只,小鸡75只
" L) `. |* C" P) y. d" U3 w 公鸡4只,母鸡18只,小鸡78只
/ q0 j; J6 b2 ^& Q8 k 公鸡8只,母鸡11只,小鸡81只$ s N6 f. a2 M! R6 @( d
公鸡12只,母鸡4只,小鸡84只
6 d% y5 n6 ]- @/ e编辑本段
5 _& G! y% |$ q. P# ~+ o* F! Jjava语言解法0 ]8 [8 d/ Z5 {, `) H O
$ S+ y8 E, L k# K; h- q6 a public class BaiJiwenti
: x8 T2 t: |0 \( E {
& u% |/ t' Z, l" N; L% S' @& D public static void main (String [] args)
% _0 j. _$ E" x- I! ? ~% ?+ v J {
4 j7 Z* \5 O+ B9 ~1 ~ for (int x = 0; x <= 19; x++)
8 x$ @, O. h' ^: ^4 p {# R) L0 k9 W# a ?" a0 g+ ~
for (int y = 0; y <= 33; y++)
' n: r! P6 f5 x6 v! d5 e; y2 _* l {; O* b" M/ \5 ^: r
int z = 100 - x - y;
5 L$ V6 f( h9 M" m" O% I7 L Y if((x * 5 + y * 3 + z / 3 == 100 ) && z % 3 == 0)3 n" r( R3 y3 c) e" ^% e! `, i
{
( e/ f' Q9 a$ E. N$ C. g" p$ N System.out.println("可买鸡翁只数:" + x);
0 d( X2 \, o4 N5 a7 n, S2 [ System.out.println("可买鸡母只数:" + y);
( K g6 l g6 s System.out.println("可买鸡雏只数:" + z);% S% J1 c7 Y; j: A
}, N1 R) ~9 r1 n0 }% F
}- m/ _# X% w+ Q# a. ]: ~2 U; q4 g
}
, J) H2 |$ ^$ @ L. Z J) K }
7 D6 u$ Q" s" |; e6 h }
, O; p7 R- o' q4 j4 D, e |
zan
|