- 在线时间
- 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 群组: 哈尔滨工业大学建模团 |
百鸡问题
! ]/ v5 \7 R! Z/ U( b y
3 L8 F2 a6 p' J 今有鸡翁一,值钱伍;鸡母一,值钱三;鸡鶵三,值钱一。凡百钱买鸡百只,问鸡翁、母、鶵各几何?答曰:鸡翁四,值钱二十;鸡母十八,值钱五十四;鸡鶵七十八,值钱二十六。又答:鸡翁八,值钱四十;鸡 母十一,值钱三十三,鸡鶵八十一,值钱二十七。又答:鸡翁十二,值钱六十;鸡母四、值钱十二;鸡鶵八十 四,值钱二十八。”; B2 z+ u$ n0 S/ T- v( @
编辑本段! Q. t. L' \" x
原书说明
g0 S s! F2 w5 ^# p
% q+ x: K; h- z7 O# ]' r 原书没有给出解法,只说如果少买7只母鸡,就可多买4只公鸡和3只小鸡。所以只要得出一组答案,就可以推出其余两组答案。中国古算书的著名校勘者甄鸾和李淳风注释该书时都没给出解法,只有约6世纪的算学家谢察微记述过一种不甚正确的解法。到了清代,研究百鸡术的人渐多,1815年骆腾风使用大衍求一术解决了百鸡问题。1874年丁取忠创用一个简易的算术解法。在此前后时曰醇(约1870)推广了百鸡问作《百鸡术衍》,从此百鸡问题和百鸡术才广为人知。百鸡问题还有多种表达形式,如百僧吃百馒,百钱买百禽等。宋代杨辉算书内有类似问题,中古时近东各国也有相仿问题流传。例如印度算书和阿拉伯学者艾布·卡米勒的著作内都有百钱买百禽的问题,且与《张邱建算经》的题目几乎全同。/ ?& F- m$ g- |* Q; v- L
编辑本段, d- W9 |( y0 C, Q+ @, d
解法+ p' x1 K8 @- A6 E) y) ~, T
0 b4 T1 v9 d5 s/ B9 {0 i4 w
中国古代算书《张丘建算经》中有一道著名的百鸡问题:公鸡每只值5 文钱,母鸡每只值3 文钱,而3 只小鸡值1 文钱。现在用100 文钱买100 只鸡,问:这100 只鸡中,公鸡、母鸡和小鸡各有多少只?
2 g: o0 L5 _- t2 R3 T 这个问题流传很广,解法很多,但从现代数学观点来看,实际上是一个求不定方程整数解的问题。解法如下:
, y5 c: {% c! ?6 ?& b1 { 设公鸡、母鸡、小鸡分别为x、y、z 只,由题意得:
! t% r8 [6 L* W- S ①……x+y+z =100( F8 o% y" I. k6 T+ N x' B3 r
②……5x+3y+(1/3)z =100+ ` r% g2 q+ o: I" A) r, J+ I2 G
有两个方程,三个未知量,称为不定方程组,有多种解。5 `- ^6 a$ G4 p: [4 }
令②×3-①得:7x+4y=100;
2 N) E2 N X3 \ 所以y=(100-7x)/4=25-2x+x/4
6 w/ n/ k6 K$ N 令x/4=t, (t为整数)所以x=4t6 g% S$ z5 v% k
把x=4t代入7x+4y=100得到:y=25-7t
3 u6 f/ ^) q+ x2 P8 M6 B0 p2 l 易得z=75+3t
) w; {/ H) ]! @; F( o 所以:x=4t5 e) m# j- c5 @. }; C
y=25-7t
8 C+ T8 K& E _! e9 J' U9 o3 S3 G% L z=75+3t; ~' J* [( h8 t8 z& U! H
因为x,y,z大于等于00 ^0 @0 U4 G+ s
所以4t大于等于01 ?5 S4 z `" @1 \! N0 r8 f; l
25-7t大于等于0/ G* [6 T1 I% z, q! W
75+3t大于等于0- ~0 K* V4 k" O2 F# s# P/ Q
解得t大于等于0小于等于25/7 又因为t为整数$ O5 S$ a. M% O# z+ ^0 V; ^) ?
所以t=0,1,2,3(这里不要忘记t有等于0得可能)
2 T3 I- {- o+ T" c1 b 当t=0时2 Y8 u) ], G4 x) D6 m
x=0,y=25,z=75
3 n; @+ D0 }6 s# m 当t=1时
! S j: e5 ]* M! r6 D; z1 R x =4;y =18;z =78
$ D2 G2 p% Y* [ 当t=2时
# H' i0 @, l* u4 `$ V% k- V7 H' B x =8;y =11;z =81
; i* K$ ?% C6 H8 U$ z0 c. ` 当t=3时7 z8 Z# H# [3 C; f. r
x =12;y =4;z =84
- R: D3 m3 B5 R8 j$ K9 V/ L% s编辑本段1 W2 O& k1 V8 {
C语言解法) B) S9 ?5 _0 G' W
1 a: `% S4 T" H& G
# ^( ^6 p" J0 I5 x1 r) i#include <stdio.h>8 X3 e6 X9 P0 R/ u# T4 f0 Z Z) R: S
void main()2 {3 b$ | |- z- z
{
8 N: p( f! g' F' U- G/ yint cocks=0,hens,chicks;2 ]+ z1 Y6 X3 f, A
while(cocks<=20)1 h5 @; c+ c( p& P/ I1 b$ {
{
/ {6 i' P( P8 O* e+ j' w2 fhens=0;
" ~7 s c* y/ A/ b2 z7 e. `: z: Pwhile(hens<=33)
' B J7 Z& V! U. h{1 [1 V5 b! Z, [, \6 o5 I
chicks=100-cocks-hens;7 ?0 d/ |+ g: c( Q) ?" w3 m4 t
if(5.0*cocks+3.0*hens+chicks/3.0==100.0)4 n6 s& A3 J" U8 N, _
printf("公鸡%d只,母鸡%d只,小鸡%d只\n\n",cocks,hens,chicks);
2 D$ z! V, S. H; `& zhens++;/ z5 e! P: s. \
}
( j" Q1 C% n8 z1 F* b! @cocks++;
; s5 Y, _' R3 z9 D# K2 n" S}* G- ]2 T# J% F3 b9 j
}
5 \; f7 f3 o* C) B# m! Q9 b E输出结果为:
$ _* V; n$ a4 Z6 e3 y! E2 A: n3 S 公鸡0只,母鸡25只,小鸡75只
" r! U/ _" A6 A2 m 公鸡4只,母鸡18只,小鸡78只
) P+ v) X3 c8 Q+ ^" B, P 公鸡8只,母鸡11只,小鸡81只
( D* s, b" A, L: U3 {. H. E$ Z' P 公鸡12只,母鸡4只,小鸡84只
. z3 ^7 ?' J, X) C- ?5 x/ m编辑本段
+ ~) ]) Y3 h7 W& gjava语言解法
3 L6 @2 I- t" g! u; w+ I. m, Z0 Z( I/ |* f
public class BaiJiwenti/ k2 J: \ W* m
{
6 o) E& c2 o N; n0 q public static void main (String [] args)
" h$ v3 Q7 y0 `1 D+ d; w: \ {- {6 p4 `3 z4 N7 l6 [. {$ o
for (int x = 0; x <= 19; x++)
5 G; o4 [" ]6 q) L {) b8 p* [5 h# L7 f4 c( e! d
for (int y = 0; y <= 33; y++)4 L3 f; ]% C' a" T6 G }1 b
{& [0 p( m2 e7 Q0 h) u0 ?8 F7 `
int z = 100 - x - y;- a3 V5 S' U, ~, l0 e
if((x * 5 + y * 3 + z / 3 == 100 ) && z % 3 == 0)$ x9 h* x# K# b/ X9 c/ n
{# `: r* t( \- V' [/ P/ I
System.out.println("可买鸡翁只数:" + x);6 o; |) N- W6 p q1 V- G: T8 P
System.out.println("可买鸡母只数:" + y);4 O6 P* e2 S k+ t5 }+ |, I9 v
System.out.println("可买鸡雏只数:" + z);- W+ y5 d+ x- }$ b2 P2 g" G6 j& m
}
1 X) r) L6 S$ s% C* y0 m }+ Y; u* I6 h1 P8 R* d7 a" v. l
}( j7 ~! g. a9 W9 r( L
}/ L% D3 b, z3 n! L- \. V+ g; ?9 H
}0 K* x0 y6 e. |1 J- I
|
zan
|