- 在线时间
- 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 群组: 哈尔滨工业大学建模团 |
百鸡问题
& t; z% |9 k0 L3 W" ?1 p7 p! [& y" g
今有鸡翁一,值钱伍;鸡母一,值钱三;鸡鶵三,值钱一。凡百钱买鸡百只,问鸡翁、母、鶵各几何?答曰:鸡翁四,值钱二十;鸡母十八,值钱五十四;鸡鶵七十八,值钱二十六。又答:鸡翁八,值钱四十;鸡 母十一,值钱三十三,鸡鶵八十一,值钱二十七。又答:鸡翁十二,值钱六十;鸡母四、值钱十二;鸡鶵八十 四,值钱二十八。”6 k5 \. g+ J3 V8 q0 k+ q% j
编辑本段5 Z2 T% M* I. w. t7 G4 |
原书说明
7 x# |. D3 \) ?+ Y" a4 g! k! b3 v% S3 H/ V, y; I V
原书没有给出解法,只说如果少买7只母鸡,就可多买4只公鸡和3只小鸡。所以只要得出一组答案,就可以推出其余两组答案。中国古算书的著名校勘者甄鸾和李淳风注释该书时都没给出解法,只有约6世纪的算学家谢察微记述过一种不甚正确的解法。到了清代,研究百鸡术的人渐多,1815年骆腾风使用大衍求一术解决了百鸡问题。1874年丁取忠创用一个简易的算术解法。在此前后时曰醇(约1870)推广了百鸡问作《百鸡术衍》,从此百鸡问题和百鸡术才广为人知。百鸡问题还有多种表达形式,如百僧吃百馒,百钱买百禽等。宋代杨辉算书内有类似问题,中古时近东各国也有相仿问题流传。例如印度算书和阿拉伯学者艾布·卡米勒的著作内都有百钱买百禽的问题,且与《张邱建算经》的题目几乎全同。 \# P/ ?! B+ Q2 B& v; V$ o
编辑本段
: U5 Z9 T/ k+ P6 D解法
6 D# h8 e- G( R* O6 z" W8 v, T, w0 e- o1 }4 j* ^
中国古代算书《张丘建算经》中有一道著名的百鸡问题:公鸡每只值5 文钱,母鸡每只值3 文钱,而3 只小鸡值1 文钱。现在用100 文钱买100 只鸡,问:这100 只鸡中,公鸡、母鸡和小鸡各有多少只?
( f; c% x9 W' n) @5 W 这个问题流传很广,解法很多,但从现代数学观点来看,实际上是一个求不定方程整数解的问题。解法如下:
7 @! J& n2 {# f/ y, F, B 设公鸡、母鸡、小鸡分别为x、y、z 只,由题意得:
9 G- q% B) Q |/ [ ①……x+y+z =100* m8 I; r( K/ Z; K% q
②……5x+3y+(1/3)z =1001 X: C7 h2 y$ v2 s, B
有两个方程,三个未知量,称为不定方程组,有多种解。
1 g I% t" M2 j( f$ }6 I+ k9 x 令②×3-①得:7x+4y=100;/ `" e' y* L1 Y9 a7 E; p
所以y=(100-7x)/4=25-2x+x/4: x! h" d. r. E) _0 s8 D
令x/4=t, (t为整数)所以x=4t, X2 ]) ~/ ^' ~) j- v- Z
把x=4t代入7x+4y=100得到:y=25-7t
9 m; }# p; Y) i- l' @ 易得z=75+3t+ N: b1 E. D& ^( o
所以:x=4t
. z6 Q% a1 X, Y* ^* G y=25-7t* ^5 o; t$ i) ~4 [6 Z7 e
z=75+3t" \: `: ~" p! g9 X3 N7 t
因为x,y,z大于等于0
, T' b* D* \; H3 }% f 所以4t大于等于0$ Q# w+ R- }7 ?9 a
25-7t大于等于08 q3 \8 e/ `( W+ j, G
75+3t大于等于0# ?: {0 b6 D5 K5 _1 `! U" i c
解得t大于等于0小于等于25/7 又因为t为整数
, a |4 _, v' @# s- W- P% h" \ 所以t=0,1,2,3(这里不要忘记t有等于0得可能)
4 ?- x, r& a$ h& P- m& Z. t 当t=0时
0 T N( p" O# i, w9 C x=0,y=25,z=756 z7 `. g" Y" Y; C7 N& ^
当t=1时8 q% z* G; I+ V) W, b3 y
x =4;y =18;z =78- X3 S9 G- B0 T( m
当t=2时! b. h7 f" i# v+ f/ O
x =8;y =11;z =81
2 k& Z( b( j' _& A% T- j9 y 当t=3时9 Y. F$ `3 c- w" y1 S. t- S0 T3 \$ {
x =12;y =4;z =84
+ Q: Z6 v7 y, f6 X* f3 q编辑本段
/ i# s; e9 ?7 w9 f8 _4 @C语言解法
' v3 s& i0 P& U6 e! U; ^6 x+ w0 x
% }- s4 q# s1 g) A) u }5 C. M% Y$ R5 \! Q% s5 ?
#include <stdio.h>
- x% ~$ p- }: e* uvoid main()4 _& q1 ]5 a8 y& P. L4 T
{
; L& @* R( M: K7 G. o2 ~' k) Yint cocks=0,hens,chicks;
( {0 t7 V. H: z1 w6 j% ^while(cocks<=20)" V- |3 ~ G+ y. c, h5 _8 v. D. F
{
$ a# Y7 R7 f6 X& qhens=0;' T! L y& x/ w: s6 }6 N# N, r
while(hens<=33)
; ^& s0 W* G( I( P0 A0 Q{
8 t- j/ \ u5 n8 Y) b" Tchicks=100-cocks-hens;% N# w; i3 M4 B- Z; }3 R) B, j! @
if(5.0*cocks+3.0*hens+chicks/3.0==100.0); U8 b4 D; }7 V; V% X) o
printf("公鸡%d只,母鸡%d只,小鸡%d只\n\n",cocks,hens,chicks);
0 V6 q* a1 e0 m0 n3 P5 }9 x, [hens++;4 V. j8 k$ \8 c5 T4 O
}) f+ @2 p, q$ I: G
cocks++;! Z" u' A( N! [
}) v% {1 w! \2 T# R
}
7 d% e& `5 [. u+ c$ v, V输出结果为:
( K. l6 z) I) T$ J) [7 ` 公鸡0只,母鸡25只,小鸡75只) c4 ~9 f; A* H5 G
公鸡4只,母鸡18只,小鸡78只
6 U) ?% E$ j& i, m 公鸡8只,母鸡11只,小鸡81只6 K! E, k" R3 J- ?2 \; d. A
公鸡12只,母鸡4只,小鸡84只4 q9 ^: K! u! E
编辑本段
( E4 o1 {( S4 X/ E1 l3 J5 ejava语言解法
1 D, U' G" B% F, ^
( c( |! F% j" s: T2 H public class BaiJiwenti/ z- p. B! M X/ O7 I4 z6 h- M! L
{0 q9 u$ b ?, a: ^( e( M5 z
public static void main (String [] args): H! V0 t2 h; l4 x" J) Q; ]
{
* P2 c4 Q7 d: ?& K. b4 y9 n for (int x = 0; x <= 19; x++)
) Y8 \5 K1 M, X' _ {
" }. K9 g2 r# f" o4 I for (int y = 0; y <= 33; y++)
2 a8 O Q, _& P4 a5 j5 Z {
* Y0 F9 i# J @ int z = 100 - x - y;
7 q9 v. r, q) R% ~ if((x * 5 + y * 3 + z / 3 == 100 ) && z % 3 == 0)
9 E$ P6 Z' D: w6 T7 y5 Z9 B {) s, l2 m4 g% |7 p3 K
System.out.println("可买鸡翁只数:" + x);* e7 Z" u! Q0 d
System.out.println("可买鸡母只数:" + y);
+ Y! }" K/ R2 O, s3 A' u System.out.println("可买鸡雏只数:" + z);
9 Y% m/ n3 z) N3 ^( d# n& a }( i3 t3 p: H- w4 Y/ |
}" n, z! o% v; Q: g7 L" C0 L
}
& s: d0 Q8 b9 F5 e0 W2 @ }: c2 Y3 r) r S1 \9 b5 x
}
; @. u9 w1 ~2 s& o' W8 S2 w0 \$ q. Z |
zan
|