QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3132|回复: 1
打印 上一主题 下一主题

百鸡问题

[复制链接]
字体大小: 正常 放大
marchboy 实名认证       

26

主题

4

听众

505

积分

升级  68.33%

  • TA的每日心情
    郁闷
    2014-11-17 00:20
  • 签到天数: 106 天

    [LV.6]常住居民II

    新人进步奖

    群组Matlab讨论组

    群组学术交流A

    群组学术交流B

    群组数学建模培训课堂1

    群组哈尔滨工业大学建模团

    跳转到指定楼层
    1#
    发表于 2012-6-9 00:10 |只看该作者 |倒序浏览
    |招呼Ta 关注Ta
    百鸡问题2 ^9 _: R/ U; b" p( g/ S# p- ~
    ; |5 h# d4 I" b# I8 U
      今有鸡翁一,值钱伍;鸡母一,值钱三;鸡鶵三,值钱一。凡百钱买鸡百只,问鸡翁、母、鶵各几何?答曰:鸡翁四,值钱二十;鸡母十八,值钱五十四;鸡鶵七十八,值钱二十六。又答:鸡翁八,值钱四十;鸡 母十一,值钱三十三,鸡鶵八十一,值钱二十七。又答:鸡翁十二,值钱六十;鸡母四、值钱十二;鸡鶵八十 四,值钱二十八。”  r* f$ ~6 ]3 D) O# l
    编辑本段) Z  l- R2 J7 |' q
    原书说明
    9 i. w( Z) S* z' ?3 v$ E! y! g8 {( C+ O7 j& e. T
      原书没有给出解法,只说如果少买7只母鸡,就可多买4只公鸡和3只小鸡。所以只要得出一组答案,就可以推出其余两组答案。中国古算书的著名校勘者甄鸾和李淳风注释该书时都没给出解法,只有约6世纪的算学家谢察微记述过一种不甚正确的解法。到了清代,研究百鸡术的人渐多,1815年骆腾风使用大衍求一术解决了百鸡问题。1874年丁取忠创用一个简易的算术解法。在此前后时曰醇(约1870)推广了百鸡问作《百鸡术衍》,从此百鸡问题和百鸡术才广为人知。百鸡问题还有多种表达形式,如百僧吃百馒,百钱买百禽等。宋代杨辉算书内有类似问题,中古时近东各国也有相仿问题流传。例如印度算书和阿拉伯学者艾布·卡米勒的著作内都有百钱买百禽的问题,且与《张邱建算经》的题目几乎全同。
    ! T# S( _7 k: ?) r4 r编辑本段$ ]) e. L* y. b) k$ j0 v; X( C
    解法
    % M  W# ?" D9 C9 \4 C( M$ _! U$ P* o4 q) N* Q
      中国古代算书《张丘建算经》中有一道著名的百鸡问题:公鸡每只值5 文钱,母鸡每只值3 文钱,而3 只小鸡值1 文钱。现在用100 文钱买100 只鸡,问:这100 只鸡中,公鸡、母鸡和小鸡各有多少只?1 y. ~' O/ J$ h8 g; E# L8 [
      这个问题流传很广,解法很多,但从现代数学观点来看,实际上是一个求不定方程整数解的问题。解法如下:
    # \! t: W3 G* o' g5 D$ u5 \: M  设公鸡、母鸡、小鸡分别为x、y、z 只,由题意得:2 B+ d% z6 M! k" P7 {7 \% a
      ①……x+y+z =100
    7 Q$ @: O3 N5 A; Q  ②……5x+3y+(1/3)z =100
    ; W7 j% n. h* B0 c  有两个方程,三个未知量,称为不定方程组,有多种解。5 P# k" P* i9 B2 S" w5 a( P( n
      令②×3-①得:7x+4y=100;" ^* x# {& l4 t/ Y" {1 N0 u* d( V
      所以y=(100-7x)/4=25-2x+x/4
    , {; p& m' X. n$ y  令x/4=t, (t为整数)所以x=4t
    ; d7 z8 y+ c3 ^( y/ Q8 f  把x=4t代入7x+4y=100得到:y=25-7t4 W: i8 r  v6 ]$ \
      易得z=75+3t
    8 L" N  d1 b. T0 l  所以:x=4t- z9 U1 X0 M5 l/ e
      y=25-7t
    - Y9 r8 ?0 a# y6 l- g- l9 u- b$ R9 T  z=75+3t+ F- K( i) G2 m, w; R6 ^  _  G" {$ x
      因为x,y,z大于等于0
    7 o( v% @3 ^: A  a& X  所以4t大于等于0: V6 i- S+ L, B1 j; |6 k
      25-7t大于等于0
    ; E: V# z8 R' m' ?" X* I) x  75+3t大于等于0
    6 m7 J- C# z$ B  K  解得t大于等于0小于等于25/7 又因为t为整数
    * Z3 W5 q5 o7 ~* R" v! z( r; s  所以t=0,1,2,3(这里不要忘记t有等于0得可能)
    / @3 K, T3 l# N# f* L3 p& D% k; R  当t=0时
    % `/ ~' y, W( h! x1 ]( E  x=0,y=25,z=75
    3 R2 v7 i" P  F4 S3 \3 ?' {8 x  当t=1时1 C& n  t0 ^. M- Q
      x =4;y =18;z =78
    " L3 e$ f' w0 B4 I$ ]  当t=2时9 r  S& t1 b- M: e
      x =8;y =11;z =81
    * H1 {7 a; a$ t  当t=3时
    ; H0 f9 {3 K) `/ N2 n# S$ B  x =12;y =4;z =84/ G  }1 p( f0 Y3 d5 c  i+ ~' ?# C9 w
    编辑本段
    ' _+ F5 g6 J( V2 ?% k& U9 yC语言解法
    7 X2 b3 v2 o# t" C1 `" [
    7 d/ C2 w* E$ V$ C  K) c- w, e" L$ ^  
    ' d5 q' ~; k' o#include <stdio.h>
    3 a; b) {1 K* w! \void main()5 Q9 c* W. [4 M* E5 m& X- I" q
    {' Y" M3 Y5 L& \! ~
    int cocks=0,hens,chicks;
    + f4 q  {" c: ~6 R6 p% k0 Ywhile(cocks<=20)
      h: W0 T  i- ^; Y4 T{
    , p: [8 q# O% E; t5 z: Thens=0;
    % b, j' Z: H0 w2 T8 r* {2 w; Rwhile(hens<=33)
    3 b9 |' M' @: g, ~. T, W8 t{5 l- K7 h1 O2 u1 b( D
    chicks=100-cocks-hens;
    - U% i8 ?1 l3 N8 ?9 }- Zif(5.0*cocks+3.0*hens+chicks/3.0==100.0)
    6 v. u/ d% R) c% a! Yprintf("公鸡%d只,母鸡%d只,小鸡%d只\n\n",cocks,hens,chicks);
    ) f; i  ]8 H3 I5 r% L: n2 ehens++;
    " D; W: _' \: Z- K}: `; l5 T$ f, V& D
    cocks++;
    / u# S- d# ^/ g1 `/ v}0 k' b7 e) v4 A* \, I2 n2 j
    }9 u( v8 n" F  }' _! Z+ I
    输出结果为:
    0 c6 j& y2 A" V" d9 ^  公鸡0只,母鸡25只,小鸡75只
    ) `' A6 V. N! m" [) j$ y& p  公鸡4只,母鸡18只,小鸡78只
    : X5 z* `+ k( a, g/ a' ]  W9 J" C  公鸡8只,母鸡11只,小鸡81只
    4 W  I6 ?! P* I  公鸡12只,母鸡4只,小鸡84只
    ! ~$ c" v# G" y8 v: w编辑本段
    " [5 o+ R$ E9 P( A; P/ Gjava语言解法6 Y' t& z: V7 a, w* m/ P% I$ M2 L

    , L  I# ^  N* F" z( v; L  public class BaiJiwenti2 ^/ P* [& R3 _
      {
    + {9 F  W1 R' G1 {7 l, `  public static void main (String [] args)& e  S& J' ?. U6 [, v! o' O3 d
      {4 \! M& h' |; u- u, o- D% L
      for (int x = 0; x <= 19; x++)' p6 [+ m; ]9 T7 U- U/ d
      {8 i% ]2 S5 r+ O) ]$ v# p
      for (int y = 0; y <= 33; y++): ?9 a% o' X0 x7 e
      {
    $ [( ?+ S9 M: z* u3 }  int z = 100 - x - y;% u+ M) M- B/ c3 A3 j  p
      if((x * 5 + y * 3 + z / 3 == 100 ) && z % 3 == 0)( D) F$ o8 N5 G! q" V7 |
      {
    , x- ?& `  I7 ?3 p" v% `: b  System.out.println("可买鸡翁只数:" + x);# r. Z1 G4 |7 V0 C$ |: g. S; }
      System.out.println("可买鸡母只数:" + y);% r# }0 T: X: d( C
      System.out.println("可买鸡雏只数:" + z);# _' x0 e/ g2 X2 d6 I6 F" T; f" j
      }
    4 \: s; N9 T6 U' W  }
    2 V, B) U3 w# m' ~2 f  @2 N  }* Z: a/ N9 y: Y9 X3 a6 v
      }
    4 Z9 B9 r( d/ R* Q( m  }
    . H+ @+ }& I$ s; T
    zan
    转播转播0 分享淘帖0 分享分享1 收藏收藏0 支持支持0 反对反对0 微信微信

    33

    主题

    10

    听众

    1691

    积分

    升级  69.1%

  • TA的每日心情
    开心
    2014-7-8 08:29
  • 签到天数: 201 天

    [LV.7]常住居民III

    发帖功臣 新人进步奖

    群组PLC和单片机

    群组2012第三期美赛培训

    群组MCM优秀论文解析专题

    群组沈阳理工应用技术学院

    群组学术交流B

    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-4-14 01:17 , Processed in 0.434459 second(s), 61 queries .

    回顶部