離散數(shù)學(xué)實(shí)驗(yàn) C ++關(guān)系的運(yùn)算(冪運(yùn)算,閉包運(yùn)算)

1、實(shí)驗(yàn)2 關(guān)系的運(yùn)算（1） 關(guān)系的冪運(yùn)算輸入：集合A，二元關(guān)系集合R，冪次n輸出：R的n次冪要求：盡量使運(yùn)算的計(jì)算量最?。?） 關(guān)系閉包的計(jì)算輸入：集合A，二元關(guān)系集合R輸出：R的傳遞閉包t(R)要求：(a) 采用Warshall 算法（89頁(yè)）(b) 編寫(xiě)代碼判斷輸出t(R)為傳遞閉包程序代碼:#include<iostream>#include<sstream>#include<vector>using namespace std;typedef vector< vector <int> > Mat;class Relationve

2、ctor<int>s;/集合Mat A;/關(guān)系矩陣Mat B;Mat C;Mat E;Mat D100; /用來(lái)存儲(chǔ)矩陣int n;public:void inputs();/將集合存入向量中void inputa();/將讀入的關(guān)系轉(zhuǎn)化為關(guān)系矩陣void print();/輸出關(guān)系矩陣void mi();int Warshall(); ;/定義類(lèi)int n,m;/全局變量，下文中使用void Relation:inputs()cout<<"輸入集合"for(int a;cin>>a;)s.push_back(a);if(getchar(

3、)='n')break;/將集合存入向量中void Relation:inputa()/將讀入的關(guān)系轉(zhuǎn)化為關(guān)系矩陣cout<<"輸入關(guān)系"int i,j,e,r;for(i=0;i<s.size();i+)vector<int> u;for(j=0;j<s.size();j+)int ia=0;u.push_back(ia);A.push_back(u);B.push_back(u);C.push_back(u);E.push_back(u);/創(chuàng)建二維向量，初始化，是每個(gè)元素為0for(int h,z;cin>&g

4、t;h>>z;) if(h=0&&z=0)break;for(i=0;i<s.size();i+)if(si=h) e=i;if(si=z) r=i;Aer=1;Ber=1;Eer=1;/Cer=1;/讀入關(guān)系，將關(guān)系對(duì)應(yīng)的矩陣中的位置元素變?yōu)?if(getchar()='n')break;void Relation:print()for(int i=0;i<s.size();i+)for(int j=0;j<s.size();j+)cout<<Aij<<" "cout<<en

5、dl;/輸出關(guān)系矩陣void Relation:mi()int a,b,i,c;cin>>n; /讀入冪次if(n=0) /0次冪for(int k=0;k<s.size();+k)for(int j=0;j<s.size();+j)if(k=j)cout<<"1 " /對(duì)角線(xiàn)上元素為1elsecout<<"0 "cout<<endl;elsefor(i=1;i<n;+i)for(int h=0;h<s.size();+h)for(int d=0;d<s.size();+d)i

6、nt m=0;for(int x=0;x<s.size();+x)m=m+Bhx*Axd; /第h行第d列的元素對(duì)應(yīng)相乘的和Chd=m;if(i>1)for(a=0;a<s.size();+a)for(b=0;b<s.size();+b)if(Cab!=D0ab) break;if(b!=s.size()break; /檢驗(yàn)是否重復(fù)if(a=s.size()&&b=s.size()break;/重復(fù)則跳出不再冪乘for(int k=0;k<s.size();k+)for(int j=0;j<s.size();j+)Bkj=Ckj;Di-1=B

7、;c=i;if(a=s.size()&&b=s.size()int q;q=(n-i)%c; /找出結(jié)果位置if(q=0) q=c;for(int e=0;e<s.size();e+)for(int f=0;f<s.size();f+)cout<<Dq-1ef<<" " /輸出cout<<endl;return;else/1次冪for(int h=0;h<s.size();h+)for(int n=0;n<s.size();n+)cout<<Bhn<<" "

8、;cout<<endl;int Relation:Warshall()for(int i=0;i<s.size();+i)for(int j=0;j<s.size();+j)if(Aji=1)for(int k=0;k<s.size();+k)Ajk=Ajk+Aik;if(Ajk!=0&&Ajk!=1)Ajk=1;print();int a=1;int b=1;/for(int p=0;p<s.size();+p)for(int l=0;l<s.size();+l)if (Apl=0)for (int x=0;x<s.size()

9、;+x)if(Apx*Axl=1)a=0;if(a=0)cout<<"wrong!"<<endl;elsefor(int p=0;p<s.size();+p)for(int l=0;l<s.size();+l)if(Apl=1&&Epl=0)Apl=0; /再判斷傳遞性for(int p=0;p<s.size();+p)for(int l=0;l<s.size();+l)if (Apl=0)for (int x=0;x<s.size();+x)if(Apx*Axl=1)b=0;if(b=1)cout<<"wrong!"<<endl;return 0;Apl=1;cout<<"right!"<<endl;/return 1;void main()Relation w;w.inp