- 在线时间
- 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 群组: 哈尔滨工业大学建模团 |
百鸡问题
( ~- ^3 P7 b2 j, r7 f$ D/ ]" A7 @% O% J* M
今有鸡翁一,值钱伍;鸡母一,值钱三;鸡鶵三,值钱一。凡百钱买鸡百只,问鸡翁、母、鶵各几何?答曰:鸡翁四,值钱二十;鸡母十八,值钱五十四;鸡鶵七十八,值钱二十六。又答:鸡翁八,值钱四十;鸡 母十一,值钱三十三,鸡鶵八十一,值钱二十七。又答:鸡翁十二,值钱六十;鸡母四、值钱十二;鸡鶵八十 四,值钱二十八。”
" o5 z% j) J7 L" G [; o- q编辑本段9 }, I1 T& g& J. @; g
原书说明) H; M* }/ K. g$ S7 _% T
) g& A& r7 k: ]7 a8 M8 ?
原书没有给出解法,只说如果少买7只母鸡,就可多买4只公鸡和3只小鸡。所以只要得出一组答案,就可以推出其余两组答案。中国古算书的著名校勘者甄鸾和李淳风注释该书时都没给出解法,只有约6世纪的算学家谢察微记述过一种不甚正确的解法。到了清代,研究百鸡术的人渐多,1815年骆腾风使用大衍求一术解决了百鸡问题。1874年丁取忠创用一个简易的算术解法。在此前后时曰醇(约1870)推广了百鸡问作《百鸡术衍》,从此百鸡问题和百鸡术才广为人知。百鸡问题还有多种表达形式,如百僧吃百馒,百钱买百禽等。宋代杨辉算书内有类似问题,中古时近东各国也有相仿问题流传。例如印度算书和阿拉伯学者艾布·卡米勒的著作内都有百钱买百禽的问题,且与《张邱建算经》的题目几乎全同。
( N9 N; }% N8 p. j) ~编辑本段
9 h8 N1 c% i% n' |解法) s, {& F5 @0 E3 a V# p: `
6 `+ A( S# [$ a( O0 ` 中国古代算书《张丘建算经》中有一道著名的百鸡问题:公鸡每只值5 文钱,母鸡每只值3 文钱,而3 只小鸡值1 文钱。现在用100 文钱买100 只鸡,问:这100 只鸡中,公鸡、母鸡和小鸡各有多少只?" U% z! L4 ?5 |: Q
这个问题流传很广,解法很多,但从现代数学观点来看,实际上是一个求不定方程整数解的问题。解法如下:
3 T M$ Q5 `7 g# u4 ?3 O 设公鸡、母鸡、小鸡分别为x、y、z 只,由题意得:
0 I8 t4 B( f* k9 c- p" `4 f# { ①……x+y+z =1009 U; O/ n. z* ^) a7 f
②……5x+3y+(1/3)z =1001 L: s) R8 E0 n/ d3 `, Z) G5 I
有两个方程,三个未知量,称为不定方程组,有多种解。! ^ b+ q6 m5 y. \
令②×3-①得:7x+4y=100;
; T$ |" \" z% d2 @& k4 n9 H$ A: \ 所以y=(100-7x)/4=25-2x+x/4
$ D8 j* t; K4 d9 P3 m. C. J 令x/4=t, (t为整数)所以x=4t! u: p, o7 Q9 I: h6 K8 n- q
把x=4t代入7x+4y=100得到:y=25-7t) a- a& @1 P% f' \+ [2 y. ]
易得z=75+3t
3 \: _. j0 M% M# K 所以:x=4t: L3 n3 ^ K" w# X- t
y=25-7t
( k! { e5 ~$ g6 ^8 }4 W& h z=75+3t" F& c0 i5 t) D. t, g6 |. ^( m
因为x,y,z大于等于0! a: \( [; u! K4 p1 O, W- w( q& U
所以4t大于等于0$ V6 {( I* e( B0 F5 E b Y
25-7t大于等于08 \7 f% V: R8 o0 ]# i( c6 [7 E% N
75+3t大于等于0$ P* k+ r0 i" |3 a$ O/ L
解得t大于等于0小于等于25/7 又因为t为整数
- g/ @4 J8 C, K! ~ 所以t=0,1,2,3(这里不要忘记t有等于0得可能)
4 C8 L3 E+ i/ y" F" z6 P7 ]: t 当t=0时
% r) J9 e" l. P2 t6 v x=0,y=25,z=75
3 {3 ?' h) B0 T$ f 当t=1时0 p3 J E# ]. M5 \& L
x =4;y =18;z =78. C/ h" n, D5 ~, }7 a
当t=2时( D; S) P9 X6 ]. ~
x =8;y =11;z =81
$ L P9 w$ n! E! `" r: ? 当t=3时
6 N3 T" m( H9 z) H+ y, k x =12;y =4;z =84 C& _& @' u9 w: s4 K
编辑本段
" @; G3 g0 N0 H8 K- YC语言解法
3 }" l' d P- n _" v m7 p4 v' c; {; ^
8 k+ T9 y( C" A#include <stdio.h># q, n% q# K) U
void main()# h& n9 N4 h$ F& Y5 k) U
{' \7 y% |/ H8 O, [9 U B/ N
int cocks=0,hens,chicks;7 P2 T* t8 }5 L }2 c7 n$ t
while(cocks<=20)
1 i3 |! W: E& K% @{
/ P! ?! J" G% @ c! P- fhens=0;- Q" A3 s* b: _3 x$ G. e1 `/ z7 B
while(hens<=33)
0 ~& W; b) k; @$ V: L1 N{
+ b# j6 u* v* Nchicks=100-cocks-hens;# W# Y3 i# _ B' d f: }( E
if(5.0*cocks+3.0*hens+chicks/3.0==100.0)
- b, I E2 D0 g$ s' N1 ]. Jprintf("公鸡%d只,母鸡%d只,小鸡%d只\n\n",cocks,hens,chicks);
' b4 M- C8 f( w% s: Y4 E- Fhens++;2 @9 w/ U* H q$ O4 C$ j; a
}
" y( L+ l" U3 a" [/ F% M% O0 Tcocks++;2 n# f1 [9 C& E- N3 C! w: M. V
}0 L- Z a, |; _1 _
}
8 K# ?7 h2 [, F* v# g输出结果为:: ?3 l5 x1 ^- H$ e3 I# O# p h
公鸡0只,母鸡25只,小鸡75只
& X6 i) l! n# M' \+ s$ D: ~ 公鸡4只,母鸡18只,小鸡78只9 i0 {. k6 D0 u! P2 g# [
公鸡8只,母鸡11只,小鸡81只5 c- R% q$ X, H3 Z
公鸡12只,母鸡4只,小鸡84只
% w8 D9 d1 c! P* ]6 n编辑本段
9 U6 _* {, s% X& {4 h6 J8 t3 r- u. Wjava语言解法1 A# {2 i. p! `4 b5 P: A0 |
. D5 C5 G1 \' t
public class BaiJiwenti
. d* c5 S! E) x& P, y) w6 ]2 A {
. P4 M4 d y; ?& F8 H' k: k public static void main (String [] args)% `. ^/ g# f, y0 }9 u' `4 G
{! J8 [* l) ~( D. H, H- i& V {% C! d
for (int x = 0; x <= 19; x++)9 e* j7 a! R$ ]
{
% H- a V! E( P5 [/ [6 ? for (int y = 0; y <= 33; y++)
# J {( V& }- L. x7 n {: O+ R4 X3 v1 M; P& f" L
int z = 100 - x - y;
1 W) _. e: `3 E8 t if((x * 5 + y * 3 + z / 3 == 100 ) && z % 3 == 0)
( Q" \/ T$ k6 Z- M {
, O3 l& h( u$ N* R System.out.println("可买鸡翁只数:" + x);9 |+ n0 k5 X6 ]6 U2 x9 Z
System.out.println("可买鸡母只数:" + y);
H3 Y9 @. z5 v& Y+ J System.out.println("可买鸡雏只数:" + z);( ~8 Q3 F' D3 |( m1 K5 o& U
}
. u) W; c, T. V5 e# ]3 ]- W! H }
6 I5 C7 K! `9 e- _: b" ` }3 G2 s7 `( }1 X& X
}
1 U, Y! z( f" |& c }; S7 r# S% o' m; ~' }
|
zan
|