- 在线时间
- 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 群组: 哈尔滨工业大学建模团 |
百鸡问题
0 q4 Z0 k4 j& B( @, T! l
3 }$ R( s0 g m& U7 d0 ` 今有鸡翁一,值钱伍;鸡母一,值钱三;鸡鶵三,值钱一。凡百钱买鸡百只,问鸡翁、母、鶵各几何?答曰:鸡翁四,值钱二十;鸡母十八,值钱五十四;鸡鶵七十八,值钱二十六。又答:鸡翁八,值钱四十;鸡 母十一,值钱三十三,鸡鶵八十一,值钱二十七。又答:鸡翁十二,值钱六十;鸡母四、值钱十二;鸡鶵八十 四,值钱二十八。”5 A4 Q+ V# h' k& Y4 W) g
编辑本段, K- b/ w5 b, p0 |' r+ H! x: K
原书说明8 h' A, C) I7 T- `
! v, Q* z0 P; s2 H1 i& l# l 原书没有给出解法,只说如果少买7只母鸡,就可多买4只公鸡和3只小鸡。所以只要得出一组答案,就可以推出其余两组答案。中国古算书的著名校勘者甄鸾和李淳风注释该书时都没给出解法,只有约6世纪的算学家谢察微记述过一种不甚正确的解法。到了清代,研究百鸡术的人渐多,1815年骆腾风使用大衍求一术解决了百鸡问题。1874年丁取忠创用一个简易的算术解法。在此前后时曰醇(约1870)推广了百鸡问作《百鸡术衍》,从此百鸡问题和百鸡术才广为人知。百鸡问题还有多种表达形式,如百僧吃百馒,百钱买百禽等。宋代杨辉算书内有类似问题,中古时近东各国也有相仿问题流传。例如印度算书和阿拉伯学者艾布·卡米勒的著作内都有百钱买百禽的问题,且与《张邱建算经》的题目几乎全同。6 o; U& Q( {- m4 G
编辑本段% S+ Z K2 C1 e! q W) e6 U# n
解法
$ e o) L* Y! A- c
- D+ W: N( _8 X8 ^4 m# a: ~6 k 中国古代算书《张丘建算经》中有一道著名的百鸡问题:公鸡每只值5 文钱,母鸡每只值3 文钱,而3 只小鸡值1 文钱。现在用100 文钱买100 只鸡,问:这100 只鸡中,公鸡、母鸡和小鸡各有多少只?
$ e4 Y" b# B' q+ M 这个问题流传很广,解法很多,但从现代数学观点来看,实际上是一个求不定方程整数解的问题。解法如下:
/ o* ~6 p+ c: f$ i$ T! V 设公鸡、母鸡、小鸡分别为x、y、z 只,由题意得:4 b- _3 Z/ ?( g% y$ Z# L
①……x+y+z =100 p7 { d0 |! E0 k
②……5x+3y+(1/3)z =100
0 [- d, u( ~; {7 ]' s 有两个方程,三个未知量,称为不定方程组,有多种解。5 d" h0 Y8 Q+ a
令②×3-①得:7x+4y=100;0 _! N2 } S+ m
所以y=(100-7x)/4=25-2x+x/4, F2 d7 M! n0 c( O1 F+ l {( J3 F
令x/4=t, (t为整数)所以x=4t
& J6 F' c2 l+ ~0 H 把x=4t代入7x+4y=100得到:y=25-7t
. U, ^* @/ f) A$ ]9 ~/ _+ A 易得z=75+3t
, E& P2 D% @" e3 l 所以:x=4t0 g& P3 S- [5 `7 P- s9 m4 ?9 |- N
y=25-7t o0 T6 N) {& J# \
z=75+3t; h* M' i% G2 O1 a/ M% A) w% B
因为x,y,z大于等于00 k- L/ y" g- r9 x i. {" O% P
所以4t大于等于03 r- s. I. i9 b: D& N( @$ N
25-7t大于等于0
( v- ~- I0 y2 ?$ ^6 P8 U# k( u 75+3t大于等于0
5 |9 ` k8 Z# m: b: |3 L% T6 q 解得t大于等于0小于等于25/7 又因为t为整数
- L4 k- U7 }4 Y# n5 D! e0 _ 所以t=0,1,2,3(这里不要忘记t有等于0得可能)
- T* ^. a) l9 I* a 当t=0时* S0 d8 M4 d3 w* A5 |
x=0,y=25,z=75
! i& L+ a$ j4 Z! ^/ b- J 当t=1时; g: H9 M& w( ]. T9 ?8 J1 q: r, |7 q
x =4;y =18;z =78. v) ~& J0 q7 n7 H# `7 b) b- e( X
当t=2时7 l/ X7 \: k5 c0 l2 o" D; T8 |" y
x =8;y =11;z =81
. G: I" Q, Z- n, a2 f0 e) k 当t=3时
+ y l; s" T. T6 |0 @* }9 K x =12;y =4;z =84; ]: Z, A9 Z$ G
编辑本段+ m. e3 f: k; I) N% J
C语言解法
' V0 Y9 N/ H9 R. Y; ?4 n/ E' T+ A% d8 D! g! E& e6 w
% J' T& Q: L4 u* o# i/ K/ N#include <stdio.h>7 N- T$ h' m! ~. m* p5 E
void main()
) S# T# q( E. a% v- \* [6 h; W{
" I4 w5 j7 s; M1 K$ {# b" x' mint cocks=0,hens,chicks;/ C$ W, z- P5 b, r" |, b! f; W
while(cocks<=20)
, J- N6 g$ c9 `4 z% ^{* w& f- R3 L z! W
hens=0;
/ c+ \& ]0 M8 o1 nwhile(hens<=33)1 n6 J9 s9 Z6 E5 w# t
{$ B! {: C: k8 M6 |
chicks=100-cocks-hens;
: b( C# p! J+ tif(5.0*cocks+3.0*hens+chicks/3.0==100.0)1 p# W' J+ |2 i9 H) {
printf("公鸡%d只,母鸡%d只,小鸡%d只\n\n",cocks,hens,chicks);
# E r. [4 S% P! H# L, t# `hens++;
) |: C7 z& j8 k8 s4 b. F V1 q, `}
, f8 ^9 |! E( i% f6 ?% D" `( rcocks++;
5 [5 F9 |/ l: G! J}2 X$ h/ u8 G8 ]# U
}
( |8 S5 Q9 {: t6 l4 j3 x% j _9 u输出结果为:
; h8 G: p1 V; | 公鸡0只,母鸡25只,小鸡75只6 i# ~2 u# ?$ z
公鸡4只,母鸡18只,小鸡78只/ c. j) a' t7 ~# k- [$ N: S: C
公鸡8只,母鸡11只,小鸡81只* `5 F# S4 I+ }# c
公鸡12只,母鸡4只,小鸡84只
4 R( X# B, h; {+ [' `1 M3 a6 @" p编辑本段
; q, ], F+ [0 I5 f9 d: tjava语言解法
6 }& L+ b3 [) K @
2 h' T$ ~" j9 ^, m" d5 c- k6 V public class BaiJiwenti n7 o- V1 J5 ]. e9 e% P
{ K* `4 v g P |# D
public static void main (String [] args); K" I, j# v# X% _9 c
{
$ D+ [- t8 U. y8 B% h for (int x = 0; x <= 19; x++)
% C* d' s! a: X8 Q& q {
5 `9 l6 h4 ?( U. G/ F$ ~ for (int y = 0; y <= 33; y++)
. L0 J K; B' n {& p" l7 |) }) L! Y# q5 J( `
int z = 100 - x - y;
' r0 y$ s1 n- H if((x * 5 + y * 3 + z / 3 == 100 ) && z % 3 == 0)
( p! j4 A& e4 D! s- P {: x4 n, ~; r0 k, N+ ]7 i4 N
System.out.println("可买鸡翁只数:" + x);) L1 u Q8 j* L3 j# f
System.out.println("可买鸡母只数:" + y);0 b" j0 `) k; {' i
System.out.println("可买鸡雏只数:" + z);5 O3 v1 f; Z* l. U! e
}8 b* n- \" Q4 w' m H/ F
}2 [" F& ?3 L' B" K% p5 O
}
/ b2 H7 Q, p% S# q2 j, k0 d$ M1 X }
( V- g% b8 A* y }- W6 z7 X# r% R, f* K
|
zan
|