离散数学（荣誉） 作业一

一、在一阶逻辑中构造下面推理的证明：

1. 学会的成员都有高级职称并且是专家。有些成员是青年人。所以有的成员是青年专家

解：
命题符号化：
$F(x):x$是学会成员 $G(x):x$是专家 $H(x):x$有高级职称 $R(x):x$是青年人
前提：$\forall x(F(x)\to G(x)\wedge H(x)),\exists x(F(x)\wedge R(x))$
结论：$\exists x(F(x)\wedge G(x)\wedge R(x))$
证明：
①$\exists x(F(x)\wedge R(x))$　　　　　　　　前提引入
②$F(y)\wedge R(y)$　　　　　　　　　　　EI规则
③$\forall x(F(x)\to G(x)\wedge H(x))$　　　　前提引入
④$F(x)\to G(x)\wedge H(x)$　　　　　　UI规则
⑤$F(y)$　　　　　　　　　　　　　　　②化简
⑥$G(y)\wedge H(y)$　　　　　　　　　　　④⑤假言推理
⑦$R(y)$　　　　　　　　　　　　　　　②化简
⑧$G(y)$　　　　　　　　　　　　　　　⑥化简
⑨$F(y)\wedge G(y)\wedge R(y)$　　　　　　　合取
⑩$\exists x(F(x)\wedge G(x)\wedge R(x))$　　　　　EG规则

2. 每个喜欢步行的人都不喜欢坐汽车. 每个人或者喜欢坐汽车或者喜欢骑自行车.有的人不喜欢骑自行车. 因而有的人不喜欢步行.（个体域人类集合）

解：
命题符号化：
$F(x):x$喜欢步行 $G(x):x$喜欢坐汽车 $H(x):x$喜欢骑自行车
前提：$\forall x(F(x)\to \neg G(x)),\forall x(G(x)\vee H(x)),\exists x\neg H(x)$
结论：$\exists x\neg F(x)$
证明：
①$\exists x\neg H(x)$ 前提引入
②$\neg H(y)$ EI规则
③$\forall x(G(x)\vee H(x))$ 前提引入
④$G(y)\vee H(y)$ UI规则
⑤$G(y)$ ②④析取三段论
⑥$\forall x(F(x)\to \neg G(x))$ 前提引入
⑦$F(y)\to \neg G(y)$ UI规则
⑧$\neg F(y)$ ⑤⑦拒取式
⑨$\exists x\neg F(x)$ EG规则

3. 东北人都不怕冷。王国端怕冷。所以王国端不是东北人。

解：
命题符号化：
$F(x):x$是东北人 $G(x):x$怕冷 $a:$王国端
前提：$\forall x(F(x)\to \neg G(x)),G(a)$
结论：$\neg F(a)$
证明：
①$\forall x(F(x)\to \neg G(x))$ 前提引入
②$F(a)\to \neg G(a)$ UI规则
③$G(a)$ 前提引入
④$\neg F(a)$ ②③拒取式

4. 每个大学毕业生不是继续读硕士就是找工作就业。每个大学毕业生只有成绩优秀才继续读硕士。有些毕业生成绩优秀，但并非每个毕业生成绩均优秀。因此，有些毕业生找工作就业。

解：
命题符号化：
集合$A：$大学生 $P(x):x$继续读硕士 $Q(x):x$找工作就业 $R(x):x$成绩优秀
前提：$\forall (x\in A)(P(x)\bigvee Q(x))\forall (x\in A)(P(x)\to R(x))\exists (x\in A)(\neg R(x))$
结论：$\exists (x\in A)(Q(x))$
证明：
①$\forall (x\in A)(P(x)\bigvee Q(x))$
②$\forall (x\in A)(P(x)\to R(x))$
③$\exists (x\in A)(\neg R(x))$
④$\neg R(y)$
⑤$\forall (x\in A)(\neg P(x)\bigvee R(x))$
⑥$\neg P(y)\bigvee R(y)$
⑦$\neg P(y)$
⑧$P(y)\bigvee Q(y)$
⑨$Q(y)$
⑩$\exists (x\in A)(Q(x))$

二、证明

1. 已知$a,b\in R,a\ne b$. 证明a和b的算数平均小于其几何平均即 $\frac{a+b}{2}>\sqrt{ab}$

解：
若证明：$\frac{a+b}{2}>\sqrt{ab}$，只需证$a+b>2\sqrt{ab}$
若证明：$a+b>2\sqrt{ab}$，只需证$(a+b{)}^2>(2\sqrt{ab}{)}^2$
化简得$a^2-2ab+b^2>0$，即$(a-b{)}^2>0$
显然成立，故得证。

2. 已知n∈N, 且n不能被2或3整除，请证明$n^2-1$能被24整除。

∵对于连续三个整数$n-1、n、n+1$。$n$不能被3整除，则$n-1$或$n+1$有且仅有一个数能够被3整除。
又∵n不能被2整除，则$n-1$和$n+1$都必能被2整除。
又∵连续两个2的倍数，必有一个数能被4整除。
∴$n-1$或$n+1$有且仅有一个数能够被4整除。
∴$(n-1)(n+1)$必能被24整除（$2\ast 3\ast 4=24$）
又∵$n^2-1=(n-1)(n+1)$
∴$n^2-1$能被24整除
故得证

3. 证明$\sqrt{3}$是无理数。

假设$\sqrt{3}$为有理数，则$\sqrt{3}$为有理数。根据有理数定义有$\sqrt{3}=m/n$（其中m与n互质）。
对等式两侧平方有：$3n^2=m^2$ 。即m为3的倍数。不妨令$m=3k$（k为整数）
化简得$n^2=3k^2$ 即n也为3的倍数。这与前提m与n互质矛盾，故假设不成立。即$\sqrt{3}$为无理数。故得证。

4. 证明存在某个整数，可以用两种不同的方式将其表示成整数的立方和。

即$a^3+b^3=c^3+d^3$（a,b,c,d互不相等）
通过计算机编程得当a,b,c,d都小于30时有4组符合条件的解。
分别如下：
$1^3+12^3=1729$ 　　 $9^3+10^3=1729$
$2^3+16^3=4104$ 　　 $9^3+15^3=4104$
$2^3+24^3=13832$ 　　 $18^3+20^3=13832$
$2^3+34^3=39312$ 　　 $15^3+33^3=39312$
故得证

5. 证明存在无理数x，y, 使得$x^y$为有理数。

当$x=e$且$y=\ln 2$时，x与y均为无理数 $x^y$为有理数。
故得证。

6. 证明可用4分和5分邮票组成n分邮资, n≥12.

∵12=3×4 13=2×4+5 14=2×5+4 15=3×5
任何大于12的整数都可以转化为12或13或14或15与若干个4的和。故得证。

三、选作题

1. 了解什么是几何平均数，平方平均数，调和平均数，算术平均数，证明它们之间的大小关系，举例说明这几种平均数在做数据统计时不同用处。

①几何平均数（$G_n$）：$G_n=\sqrt[n]{x_1{x}_2{x}_3\cdots x_n}$
几何平均数是对各变量值的连乘积开项数次方根。如果总水平、总成果等于所有阶段、所有环节水平、成果的连乘积总和时，求各阶段、各环节的一般水平、一般成果，此时将使用几何平均数。
用途：
（1）对比率、指数等进行平均；
（2）计算平均发展速度
（3）复利下的平均年利率；
（4）连续作业的车间求产品的平均合格率。
特点：
（1）几何平均数受极端值的影响较算术平均数小；
（2）如果变量值有负值，计算出的几何平均数可能会成为负数或虚数；
（3）它仅适用于具有等比或近似等比关系的数据；
（4）几何平均数的对数是各变量值对数的算术平均数。
例：假定某地储蓄年利率（按复利计算）：5%持续1.5年，3%持续2.5年，2.2%持续1年。求此5年内该地平均储蓄年利率。
解：平均储蓄年利率：$G=\sqrt[1.5+2.5+1]{1.05^{1.5}\times 1.03^{2.5}\times 1.022^1}-1\approx 3.43\%$

②调和平均数（$H_n$）：$H_n=\frac{n}{\frac{1}{x_1}+\frac{1}{x_2}+\frac{1}{x_3}+\cdots +\frac{1}{x_n}}$
调和平均数又称倒数平均数，是总体各统计变量倒数的算术平均数的倒数。调和平均数是平均数的一种。但统计调和平均数，与数学调和平均数不同，它是变量倒数的算术平均数的倒数。由于它是根据变量的倒数计算的，所以又称倒数平均数。
特点：
（1）调和平均数易受极端值的影响，且受极小值的影响比受极大值的影响更大。
（2）只要有一个标志值为0，就不能计算调和平均数。
（3）当组距数列有开口组时，其组中值即使按相邻组距计算，假定性也很大，这时的调和平均数的代表性很不可靠。
（4）调和平均数应用的范围较小。在实际中，往往由于缺乏总体单位数的资料而不能直接计算算术平均数，这时需用调和平均法来求得平均数。

例：三个分别为10Ω，20Ω，30Ω的电阻并联，求等效电阻
解：$H=\frac{1}{\frac{1}{10}+\frac{1}{20}+\frac{1}{30}}=\frac{60}{11}\Omega$

③算数平均数（$A_n$）：$A_n=\frac{x_1+x_2+x_3+\cdots +x_n}{n}$
算术平均数又称均值，是统计学中最基本、最常用的一种平均指标，分为简单算术平均数、加权算术平均数。它主要适用于数值型数据，不适用于品质数据。根据表现形式的不同，算术平均数有不同的计算形式和计算公式。
特点：
（1）算术平均数是一个良好的集中量数，具有反应灵敏、确定严密、简明易解、计算简单、适合进一步演算和较小受抽样变化的影响等优点。
（2）算术平均数易受极端数据的影响，这是因为平均数反应灵敏，每个数据的或大或小的变化都会影响到最终结果。

例：某销售小组有5名销售员，元旦一天的销售额分别为520元、600元、480元、750元和500元，求该日平均销售额。
解：$A=\frac{520+600+480+750+500}{5}=570$

④平方平均数（$Q_n$）：$Q_n=\sqrt{\frac{x_1^2+x_2^2+x_3^2+\cdots +x_n^2}{n}}$
平方平均数，又名均方根，是指一组数据的平方的平均数的算术平方根。
特点：
（1）方均根值并非所有模型均适用， 只有在数值分布呈现正态分布时才适用。
（2）如果分布呈现方波、三角波，那就要用其他的公式， 否则失真会很大。

例：已知某次期中考试学生成绩为64,80,95,60,42,100,95,92,97,79求标准差
解：标准差$Q=\sqrt{\frac{(64-M{)}^2+(80-M{)}^2+(95-M{)}^2+\cdots +(79-M{)}^2}{10}}\approx 18.4$

⑤大小关系：
调和平均数≤几何平均数≤算术平均数≤平方平均数（$H_n\le G_n\le A_n\le Q_n$）
$$\frac{n}{\frac{1}{x_1}+\frac{1}{x_2}+\frac{1}{x_3}+\cdots +\frac{1}{x_n}}\le \sqrt[n]{x_1{x}_2{x}_3\cdots x_n}\le \frac{x_1+x_2+x_3+\cdots +x_n}{n}\le \sqrt{\frac{x_1^2+x_2^2+x_3^2+\cdots +x_n^2}{n}}$$

证明如下：
对于$\sqrt[n]{x_1{x}_2{x}_3\cdots x_n}\le \frac{x_1+x_2+x_3+\cdots +x_n}{n}$有：
当$n=1,2$时，不等式显然成立；
当$n=2^k(k\in N^{\ast })$ 时，不等式为$\sqrt{ab}\le \frac{a+b}{2}$的直接推论
当$n\ne 2^k(k\in N^{\ast })$时，不妨设$l\in N^{\ast }$，使得$2^{l-1}<n<2^l$。令
$$\sqrt[n]{x_1{x}_2{x}_3\cdots x_n}=\bar{x}$$
在$x_1,x_2,x_3\dots x_n$后补上$(2^l-n)$个$\bar{x}$，将其扩充为$2^l$个数。此时有：
$$\frac{1}{2^l}[x_1+x_2+x_3+\cdots +x_n+(2^l-n)\bar{x}]\ge (x_1{x}_2{x}_3\cdots x_n{\bar{x}}^{2^l-n}{)}^{\frac{1}{2^l}}=\bar{x}$$
整理得：
$$\sqrt[n]{x_1{x}_2{x}_3\cdots x_n}\le \frac{x_1+x_2+x_3+\cdots +x_n}{n}$$
对于上式，不妨用$\frac{1}{x_i}$替换$x_i$此时有：
$$\sqrt[n]{\frac{1}{x_1}\frac{1}{x_2}\frac{1}{x_3}\cdots \frac{1}{x_n}}\le \frac{\frac{1}{x_1}+\frac{1}{x_2}+\frac{1}{x_3}+\cdots +\frac{1}{x_n}}{n}$$
因为$x_1,x_2,x_3\cdots x_n$都大于0，故取倒数并整理后有：
$$\frac{n}{\frac{1}{x_1}+\frac{1}{x_2}+\frac{1}{x_3}+\cdots +\frac{1}{x_n}}\le \sqrt[n]{x_1{x}_2{x}_3\cdots x_n}$$
依琴生不等式，当$f\left(x\right)=\sqrt{x}$得：

$$\frac{\sqrt{x_1}+\sqrt{x_2}+\sqrt{x_3}+\dots +\sqrt{x_n}}{n}\le \sqrt{\frac{x_1+x_2+x_3+\cdots +x_n}{n}}$$

即：
$$\frac{x_1+x_2+x_3+\cdots +x_n}{n}\le \sqrt{\frac{x_1^2+x_2^2+x_3^2+\cdots +x_n^2}{n}}$$
因此：

$$\frac{n}{\frac{1}{x_1}+\frac{1}{x_2}+\frac{1}{x_3}+\cdots +\frac{1}{x_n}}\le \sqrt[n]{x_1{x}_2{x}_3\cdots x_n}\le \frac{x_1+x_2+x_3+\cdots +x_n}{n}\le \sqrt{\frac{x_1^2+x_2^2+x_3^2+\cdots +x_n^2}{n}}$$

即调和平均数≤几何平均数≤算术平均数≤平方平均数（$H_n\le G_n\le A_n\le Q_n$）

2、自学prolog语言，编程实现，找出谁是凶手，并提交完整代码的实验报告。

已知条件：

1. 六个嫌疑犯是三男（George、John、Robert）三女（Barbara、Christine、Yolanda）。
2. 六个嫌疑犯分别待在六个房间：浴室（Bathroom）、饭厅（Dining Room）、厨房（Kitchen）、起居室（Living Room）、 储藏室（Pantry）、书房（Study）。每间房间都有一件可疑的物品，可以当作凶器：包（Bag）、火枪（Firearm）、煤气（Gas）、刀（Knife）、毒药（Poison）、绳索（Rope）。
3. 厨房里面是一个男人，那里的凶器不是绳索、刀子、包和火枪。
4. Barbara 和 Yolanda 在浴室和书房。
5. 带包的那个人不是 Barbara 和 George，也不在浴室和饭厅。
6. 书房里面是一个带绳子的女人。
7. 起居室里面那件凶器，与 John 或 George 在一起。
8. 刀子不在饭厅。
9. 书房和食品储藏室里面的凶器，没跟 Yolanda 在一起。
10. George 所在的那间屋子有火枪。
11. Boddy 先生死在食品储藏室里，那里的凶器是煤气。

man(george). man(john). man(robert).
woman(barbara). woman(christine). woman(yolanda).
person(X):- man(X).
person(X):- woman(X).

location(bathroom). location(dining). location(kitchen).
location(livingroom). location(pantry). location(study).
weapon(bag). weapon(firearm). weapon(gas). 
weapon(knife). weapon(poison). weapon(rope).

uniq_ppl(A,B,C,D,E,F):- person(A), person(B), person(C), person(D), person(E), person(F),  \+A=B, \+A=C, \+A=D, \+A=E, \+A=F, \+B=C, \+B=D, \+B=E, \+B=F, \+C=D, \+C=E, \+C=F, \+D=E, \+D=F, \+E=F.

writeanswers(Bathroom, Dining, Kitchen, Livingroom, Pantry, Study, Bag, Firearm, Gas, Knife, Poison, Rope):- 
    write("Bathroom: "), write(Bathroom), nl,
    write("Dining: "), write(Dining), nl,
    write("Livingroom: "), write(Livingroom), nl, 
    write("Pantry: "), write(Pantry), nl,
    write("Study: "), write(Study), nl,
    write("Kitchen: "), write(Kitchen), nl,

    write("Knife: "), write(Knife), nl,
    write("Gas: "), write(Gas), nl,
    write("Rope: "), write(Rope), nl, 
    write("Bag: "), write(Bag), nl,
    write("Poison: "), write(Poison), nl,
    write("Firearm: "), write(Firearm), nl.

clue1(Bathroom, Dining, Kitchen, Livingroom, Pantry, Study, Bag, Firearm, Gas, Knife, Poison, Rope) :-
   man(Kitchen), 
   \+Kitchen=Rope, \+Kitchen=Knife, \+Kitchen=Bag, \+Kitchen=Firearm.

clue2(Bathroom, Dining, Kitchen, Livingroom, Pantry, Study, Bag, Firearm, Gas, Knife, Poison, Rope) :-
    barbara=Bathroom, yolanda=Study.
clue2(Bathroom, Dining, Kitchen, Livingroom, Pantry, Study, Bag, Firearm, Gas, Knife, Poison, Rope) :-
    barbara=Study, yolanda=Bathroom.

clue3(Bathroom, Dining, Kitchen, Livingroom, Pantry, Study, Bag, Firearm, Gas, Knife, Poison, Rope) :-
    \+barbara=Bag, \+george=Bag, \+Bag=Bathroom, \+Bag=Dining.

clue4(Bathroom, Dining, Kitchen, Livingroom, Pantry, Study, Bag, Firearm, Gas, Knife, Poison, Rope) :-
    woman(Rope), Rope=Study.

clue5(Bathroom, Dining, Kitchen, Livingroom, Pantry, Study, Bag, Firearm, Gas, Knife, Poison, Rope) :-
    Livingroom = john.
clue5(Bathroom, Dining, Kitchen, Livingroom, Pantry, Study, Bag, Firearm, Gas, Knife, Poison, Rope) :-
    Livingroom = george.

clue6(Bathroom, Dining, Kitchen, Livingroom, Pantry, Study, Bag, Firearm, Gas, Knife, Poison, Rope) :-
    \+Knife=Dining.

clue7(Bathroom, Dining, Kitchen, Livingroom, Pantry, Study, Bag, Firearm, Gas, Knife, Poison, Rope) :-
    \+yolanda=Pantry, \+yolanda=Study.

clue8(Bathroom, Dining, Kitchen, Livingroom, Pantry, Study, Bag, Firearm, Gas, Knife, Poison, Rope) :-
    Firearm=george.

murderer(X) :-
   uniq_ppl(Bathroom, Dining, Kitchen, Livingroom, Pantry, Study),
   uniq_ppl(Bag, Firearm, Gas, Knife, Poison, Rope),
   clue1(Bathroom, Dining, Kitchen, Livingroom, Pantry, Study, Bag, Firearm, Gas, Knife, Poison, Rope),
   clue2(Bathroom, Dining, Kitchen, Livingroom, Pantry, Study, Bag, Firearm, Gas, Knife, Poison, Rope),
   clue3(Bathroom, Dining, Kitchen, Livingroom, Pantry, Study, Bag, Firearm, Gas, Knife, Poison, Rope),
   clue4(Bathroom, Dining, Kitchen, Livingroom, Pantry, Study, Bag, Firearm, Gas, Knife, Poison, Rope),
   clue5(Bathroom, Dining, Kitchen, Livingroom, Pantry, Study, Bag, Firearm, Gas, Knife, Poison, Rope),
   clue6(Bathroom, Dining, Kitchen, Livingroom, Pantry, Study, Bag, Firearm, Gas, Knife, Poison, Rope),
   clue7(Bathroom, Dining, Kitchen, Livingroom, Pantry, Study, Bag, Firearm, Gas, Knife, Poison, Rope),
   clue8(Bathroom, Dining, Kitchen, Livingroom, Pantry, Study, Bag, Firearm, Gas, Knife, Poison, Rope),
   X = Gas, X = Pantry,
   
   write("KILLER IS :"), write(X), nl, 
   writeanswers(Bathroom, Dining, Kitchen, Livingroom, Pantry, Study, Bag, Firearm, Gas, Knife, Poison, Rope).

运行结果：

KILLER IS :christine
Bathroom: yolanda
Dining: george
Livingroom: john
Pantry: christine
Study: barbara
Kitchen: robert
Knife: yolanda
Gas: christine
Rope: barbara
Bag: john
Poison: robert
Firearm: george
X = christine

即凶手是christine