离散数学真值表c语言实验报告,NJUPT【离散数学】实验报告

实验一:利用真值表法求主析取范式、主合取范式

Using truthTable to find the pdnf and pcnf.cpp

#include

#define MAX 100000

using namespace std;

short trueValue[MAX]; // 真值

short trueTable[MAX][10]; // 真值表

int n, i;

// cnt1 是真值为 T 命题公式的数目

// cnt2 是真值为 F 命题公式的数目

int cnt1 = 0, cnt2 = 0;

// A)主析取范式

void Output_pdnf() {

cout << "主析取范式为 : " << endl;

string ans = "";

int tmpcnt = 0;

// 遍历真值表每一行

for (int i = 0; i < (1 << n); i ++) {

// 若此行命题公式真值为 T

if (trueTable[i][n] == 1) {

tmpcnt ++; // 更新 tmpcnt

ans += '('; // 加左括号 "("

// 遍历此行的各个命题

for (int j = 0; j < n; j ++) {

// 如果命题真值为 F,加关联词 "┐"

if (trueTable[i][j] == 0) {

ans += "┐ ";

}

// 加命题变元 j+'P' 对应的 ASCII码

ans += (j + 'P');

// 若命题变元没遍历完,加关联词 "∧"

if (j != n - 1) {

ans += " ∧ ";

}

}

ans += ')'; // 加右括号 ")"

// 若真值为 T 的命题公式没遍历完,加关联词 "V"

if (tmpcnt < cnt1) {

ans += " V ";

}

}

}

// 输出主析取范式的字符串

cout << ans << endl;

}

// B)主合取范式

void Output_pcnf() {

cout << "主合取范式为 : " << endl;

string ans = "";

int tmpcnt = 0;

// 遍历真值表每一行

for (int i = 0; i < (1 << n); i ++) {

// 若对应命题真值为 F

if (trueTable[i][n] == 0) {

tmpcnt ++; // 更新命题个数

ans += '('; // 加左括号 "("

// 遍历此行的各个命题

for (int j = 0; j < n; j ++) {

// 如果命题真值为 T,加关联词 "┐"

if (trueTable[i][j] == 1) {

ans += "┐ ";

}

// 加命题变元 "P"

ans += (j + 'P');

// 若命题变元没遍历完,加关联词 " V "

if (j != n - 1) {

ans += " V ";

}

}

ans += ')'; // 加右括号 ")"

// 若真值为 F 的命题公式没遍历完,加关联词 "∧"

if (tmpcnt < cnt2) {

ans += " ∧ ";

}

}

}

// 输出主析取范式的字符串

cout << ans << endl;

}

int main() {

cout << "请输入命题变元的个数 :" << endl;

cin >> n;

cout << "请输入命题公式的真值 :" << endl;

for (i = 0; i < (1 << n); i ++) {

// 输入各命题真值

cin >> trueValue[i];

// 若命题真值为 T,更新 cnt1,pcnf

if (trueValue[i]) {

cnt1 ++;

}

// 若命题真值为 F,更新 cnt2,pdnf

else {

cnt2 ++;

}

}

// 给真值表各个命题真值赋值

for (i = 0; i < (1 << n); i ++) {

int tmp = i;

// 保证全为 0 或 1,即 T 或 F

for (int j = n - 1; j >= 0; j--) {

trueTable[i][j] = tmp % 2;

tmp /= 2;

}

}

// 给真值表每行命题公式赋值

for (i = 0; i < (1 << n); i++) {

trueTable[i][n] = trueValue[i];

}

// 输出主析取范式

Output_pdnf();

// 输出主合取范式

Output_pcnf();

return 0;

}

实验二:集合上二元关系性质的判断与实现

Judge properties of binary relation on set.cpp

# include

# include

# include // 控制台输入输出库

// A)读取元素集合

char *get_element(char *p) {

printf("输入集合的(不能有空格):");

// 读取字符串 p

gets(p);

// 清空输入缓冲区,以便下次读取

fflush(stdin);

return p;

}

// B)获取元素在集合中的位置

int get_position(char ch, char *point) {

for(int i = 0; *(point + i); i++) {

// 返回 ch 在 point 中的下标

if(*(point + i) == ch)

return i;

}

return 0;

}

// C)读取序偶,构建关系矩阵

void get_relation(int (*a)[100], char *p) {

int k1, k2;

char ch1, ch2;

printf("输入关系的各个序偶(以 结束):\n");

while(1) {

printf("

// 输入后立即读取 ch1,不以回车结束

ch1 = getche();

printf(",");

// 输入后立即读取 ch2,不以回车结束

ch2 = getche();

printf(">\n");

if(ch1 == '*' && ch2 == '*') // 输入序偶 读取结束

break;

k1 = get_position(ch1, p); // 取得 ch1 在 p 中的位置

k2 = get_position(ch2, p); // 取得 ch2 在 p 中的位置

// 将关系矩阵相应元素置 1

a[k1][k2] = 1;

}

}

// D)输出关系矩阵

void output_relat_array(int (*a)[100], int arry_w) {

for(int i = 0; i

for(int j = 0; j

printf("%4d", a[i][j]);

printf("\n");

}

}

// E)输出序偶集合

void output_relate(int (*a)[100], int arry_w,char *p) {

int count=0;

printf("{ ");

// 遍历关系矩阵

for(int i = 0; i < arry_w; i++)

for(int j = 0; j

// 若相应元素为 1,则返回对应序偶

if(a[i][j] == 1)

printf(",",*(p+i), *(p+j));

count++;

// \b表示删除一个字符

printf("\b }");

printf("\n");

}

/*

1)自反性

*/

int ZF(int (*a)[100], int n) {

// 初始化判断

int flag1 = 1;

for(int i = 0; i < n; i++) {

// 若存在对角元素为 0

if(!a[i][i]) {

// 则不具有自反性

flag1 = 0;

break;

}

}

return flag1;

}

/*

2)反自反性

*/

int FZF(int (*a)[100], int n) {

// 初始化判断

int flag2 = 1;

for(int i = 0; i < n; i++) {

// 若存在对角元素为 1

if(a[i][i]) {

// 则不具有反自反性

flag2 = 0;

break;

}

}

return flag2;

}

/*

3)对称性

*/

int DC(int (*a)[100], int n) {

// 初始化判断

int flag3 = 1;

for(int i = 0; i < n; i++)

for(int j = 0; j < n; j++) {

// 若存在对称元素不相等

if(a[i][j] != a[j][i] && i != j)

// 则不具有对称性

flag3 = 0;

break;

}

return flag3;

}

/*

4)反对称性

*/

int FDC(int (*a)[100],int n) {

int flag4 = 1;

for(int i = 0; i < n; i++)

for(int j = 0; j < n; j++)

// 若存在对称元素都为 1

if(a[i][j] && a[j][i] && i != j) {

// 则不具有反对称性

flag4 = 0;

break;

}

return flag4;

}

/*

5)传递性

*/

int CD(int (*a)[100], int n) {

// 初始化判断

int flag5 = 1;

// 遍历序偶所有可能的传递组合

for(int i = 0; i < n; i++)

for(int j = 0; j < n; j++)

for(int k = 0; k < n; k++)

// 若不满足传递定义

if(a[i][j] && a[j][k] && !a[i][k]) {

// 则不具有传递性

flag5 = 0;

break;

}

return flag5;

}

int main() {

int a[100][100] = {0}; // 初始化关系矩阵所有元素为 0

char point[100]; // 元素集合 point

int stlen; // 元素集合个数 stlen

char *p; // 元素集合指针 p

p = get_element(point); // 读取 point,得到 p

stlen = strlen(point);

get_relation(a, p); // 读取序偶,构建关系矩阵

printf("序偶集合为:");

output_relate(a, stlen, p); // 输出序偶集合

printf("关系矩阵为:\n");

output_relat_array(a,stlen); // 输出关系矩阵

printf("该关系具有的性质:");

if(ZF(a, stlen)) {

printf("自反性 ");

}

if(FZF(a,stlen)) {

printf("反自反性 ");

}

if(DC(a,stlen)) {

printf("对称性 ");

}

if(FDC(a,stlen)) {

printf("反对称性 ");

}

if(CD(a,stlen)) {

printf("传递性 ");

}

return 0;

}

实验三:偏序关系中求盖住关系,格伦中判定有补格

Covering relation and Judging complementation.cpp

#include

#include

#include

#include

using namespace std;

int count = 0; // 从 0 开始计数

int n; // 正整数

int factors[100]; // 存因子的数组

int matrixs[100][100] = {0}; // 关系矩阵初始化为 0

void factor(){

cout << "The factors of " << n << " are: ";

for(int i = 1; i <= n/2; i++){

if(n % i == 0){

factors[count++] = i;

cout << i << ", ";

}

}

factors[count] = n;

cout << n << endl << endl;

}

void cover(){

for(int i = 0; i <= count; i++){

for(int j = 0; j <= count; j++){

if(factors[j] % factors[i] == 0){ // 如果满足整除,设为 1

matrixs[i][j] = 1;

}

}

}

// 输出偏序关系

cout << "The partial relation is: { ";

for(int i = 0; i <= count; i++){

for(int j = 0; j <= count; j++){

if(matrixs[i][j]){

cout << " ";

}

}

}

cout << " }" << endl << endl;

for(int i = 0; i <= count; i++){

for(int j = 0; j <= count; j++){

for(int k = 0; k <= count; k++){

matrixs[k][k] = 0; // 去掉自反性,设为 0

if(matrixs[i][j] && matrixs[j][k]){

matrixs[i][k] = 0; // 去掉传递性,设为 0

}

}

}

}

// 输出盖住关系

cout << "The covering relation is: { ";

for(int i = 0; i <= count; ++i){

for(int j = 0; j <= count; ++j){

if(matrixs[i][j]){

cout << " ";

}

}

}

cout << " }" << endl << endl;

}

// 求最大公因子(辗转相除法)

int gcd(int x, int y){

int m;

while(m != 0){

m = x % y;

x = y;

y = m;

}

return x;

}

void complemented_lattice(){

bool flag;

int Gcd, Lcm;

for(int i = 1; i < count; i++){

flag = false; // 初始化 flag = false

for(int j = 1; j < count; j++){

if(i == j)

continue;

Gcd = gcd(factors[i], factors[j]); // 最大公约数,即最大下界

Lcm = factors[i]*factors[j] / Gcd; // 最小公倍数,即最小上界

if(Gcd == 1 && Lcm == n){ // 若是补元,flag = true

flag = true;

break;

}

if(!flag){ // 若存在元素没有补元,则此格不是有补格

cout << "This is not complemented lattice!" << endl;

return;

}

}

}

// 若所有元素均有补元,则此格是有补格

cout << "This is complemented lattice!" << endl;

return;

}

int main(){

cout << "Please input a positive integer: ";

cin >> n;

cout << endl;

factor(); // 求 n的因子,并保存

cover(); // 求盖住关系,并输出

complemented_lattice(); // 判断有补格

return 0;

}

实验四:图的随机生成,欧拉(回)路的确定

Random generation of graphs and Euler loop.cpp

#include

#include

#include

#include

#define Size 1000

using namespace std;

int n, m; // 图的结点,图的边

int G[Size][Size]; // 图的邻接矩阵

void Generate(){

printf("\n>>> 正在生成...\n",n,m);

int cnt=0;

srand(time(NULL));

while(cnt

int x=rand()%n;

int y=rand()%n;

if(x!=y && G[x][y]==0){

G[x][y]=1;

G[y][x]=1;

cnt++;

}

}

printf(">>> 生成完成!!!\n\n");

printf(">>> 随机生成图的邻接矩阵 G=\n");

for(int i=0;i

for(int j=0;j

printf("%d ",G[i][j]);

}

printf("\n");

}

printf("\n");

}

int Judge() {

int flag=0;

for(int i=0;i

int cnt=0;

for(int j=0;j

if(G[i][j]==1){

cnt++;

}

}

if(cnt%2==1){ // 若第 i行有奇数个 1

flag++;

}

}

if(flag==0){ // 若没有行有奇数个 1

return 0; // 欧拉回路

}

else if(flag==2){ // 若有两行有奇数个 1

return 1; // 欧拉路

}

else{

return -1;

}

}

int P[Size][Size]; // 可达性矩阵

int T[Size][Size]; // 存放 G

int TT[Size][Size]; // 存放 G的 K次方

int JudgeLianTong(){

for(int i=0;i

for(int j=0;j

P[i][j]=G[i][j];

T[i][j]=G[i][j];

}

}

for(int k=2;k<=n;k++){ //n的4次方复杂度,计算可达性矩阵

for(int i=0;i

for(int j=0;j

int t=0;

for(int a=0;a

t+=T[i][a]*G[a][j];

}

if(t==0){

TT[i][j]=0;

}

else{

TT[i][j]=1;

}

}

}

for(int i=0;i

for(int j=0;j

T[i][j]=TT[i][j];

}

}

for(int i=0;i

for(int j=0;j

if(T[i][j]>0 || P[i][j]>0 ){

P[i][j]=1;

}

}

}

}

printf(">>> 随机生成图的可达性矩阵 P=\n");

for(int i=0;i

for(int j=0;j

printf("%d ",P[i][j]);

}

printf("\n");

}

printf("\n");

for(int i=0;i

for(int j=0;j

if(i!=j && P[i][j] ==0){

return 0;

}

}

}

return 1;

}

int choice,has,cnt; // cnt为欧拉回路个数

int vis[Size][Size];

int record[Size];

void FindLu(int cur){

if(choice==1 && has==1) return;

if(cnt==m+1){

for(int i=0;i

if(i==0) printf("%d",record[i]);

else printf("->%d",record[i]);

}

printf("\n");

has=1;

}

else{

for(int i=0;i

if(G[cur][i]==1 && vis[cur][i]==0 ){

vis[i][cur]=vis[cur][i]=1;

record[cnt++]=i;

FindLu(i);

cnt--;

vis[i][cur]=vis[cur][i]=0;

}

}

}

}

int main(){

do{

printf(">>> 请输入无向图的结点个数:"); scanf("%d",&n);

printf(">>> 请输入无向图的边的个数:"); scanf("%d",&m);

if(m>n*(n-1)/2){

printf("\n>>> %d个结点的无向图最多有%d条边!\n\n",n,n*(n-1)/2);

}

}while(m>n*(n-1)/2);

Generate(); //随机生成图

if(JudgeLianTong()==0){ //判断连通性

printf(">>> 这个图不是一个连通图,所以也不是欧拉图和半欧拉图\n");

}

else{

printf("\n>>> 这个图是连通图\n");

int tmp=Judge(); //判断欧拉图

if(tmp==0){

printf(">>> 这个图是一个欧拉图\n");

printf("\n-----------------1.打印一个欧拉回路------------------\n");

printf("-----------------2.打印所有欧拉回路------------------\n\n");

printf(">>> 请输入你的选择:"); scanf("%d", &choice);

if(choice==1){

printf(">>> 其中一欧拉回路为:\n");

record[cnt++]=0;

FindLu(0); // 输出

cnt--;

}

else if(choice==2){

printf(">>> 所有的欧拉回路为:\n");

for(int i=0;i

record[cnt++]=i;

FindLu(i); // 输出

cnt--;

}

}

}

else if(tmp==1){

printf(">>> 这个图是一个半欧拉图\n");

printf("\n-----------------1.打印一个欧拉路------------------\n");

printf("-----------------2.打印所有欧拉路------------------\n\n");

printf(">>> 请输入你的选择:"); scanf("%d",&choice);

if(choice==1){

for(int i=0;i

int t=0;

for(int j=0;j

if(G[i][j]==1){

t++;

}

}

if(t%2==1){ // 若此行有奇数个 1

record[cnt++]=i;

FindLu(i); // 输出,循环中断

cnt--;

break;

}

}

}

else if(choice==2){

printf(">>> 所有的欧拉路为:\n");

for(int i=0;i

int t=0;

for(int j=0;j

if(G[i][j]==1){

t++;

}

}

if(t%2==1){ // 若此行有奇数个 1

record[cnt++]=i;

FindLu(i); // 输出,循环继续

cnt--;

}

}

}

}

else{

printf(">>> 这个图既不是欧拉图,也不是半欧拉图\n");

}

}

return 0;

}

你可能感兴趣的:(离散数学真值表c语言实验报告)