sushu.rar_K.

#include<iostream.h> #include<math.h> int gcd(int a,int b) { int r; if(a<b) { int temp=a; a=b; b=temp;} while(b!=0) { r=a%b; a=b; b=r;} return a;} bool prime(int p) { int k=sqrt(p+1); for(int i=2;i<=k;i++) if(p%i==0) return false; return true;} void main() { int i,j,D,p,max,m[3][3]; cout<<"please input the elements of matrix"<<endl; for(i=0;i<3;i++) for(j=0;j<3;j++) cin>>m[i][j]; D=m[0][0]*(m[1][1]*m[2][2]-m[2][1]*m[1][2])-m[0][1]*(m[1][0]*m[2][2]-m[2][0]*m[1][2]) +m[0][2]*(m[1][0]*m[2][1]-m[2][0]*m[1][1]);cout<<D<<endl; cout<<"the maximum p :"; cin>>max; for(j=2;j<=max;j++) { p=gcd(D,j); if(prime(j)&&p!=1) cout<<"p="<<j<<","; }}