盖住关系的验证

//先将n所有因子求出保存在a数组中，如果a[i]是a[j]的整数倍，那么如果j<k<i并且a[k]是a[j]整数倍，
//a[i]是a[k]的整数倍,那么这组a[i],a[j]就可以忽略。所谓有补格就是任意一个因子都能找到另一个因子使得
//他们最大公约数是1，最小公倍数是n即可。
#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <cmath>
using namespace std;
int ok;
void caculate(int a,int b,int n)
{
int r,m,k1,k2;
k1=a;
k2=b;
do
{
r=a%b;
a=b;
b=r;
}
while(r);
if ((a==1)&&(k1*k2/a==n))
ok=1;
}
int main()
{
int i,j,m,k,n,sum;//sum记录该数有几个因子
int a[1000];
printf("请输入整数n\n");
while(scanf("%d",&n)==1&&n)
{
memset(a,0,sizeof(a));
sum=0;
k=(int)sqrt(n);
for (i=1;i<=k;i++)
if (n%i==0)
{
a[sum++]=i;
a[sum++]=n/i;
}
if (k*k==n)
sum--;//如果根号正好为整数
sort(a,a+sum);
for (i=0;i<sum-1;i++)
for (j=i+1;j<sum;j++)
{
ok=0;
if (a[j]%a[i]==0)
{
ok=1;
for (k=i+1;k<j;k++)
if ((a[k]%a[i]==0)&&(a[j]%a[k]==0))
ok=0;
}
if (ok)
printf("<%d %d>\n",a[i],a[j]);
}
for (i=0;i<sum;i++)
{
ok=0;
for (j=0;j<sum;j++)
caculate(a[i],a[j],n);
if (!ok)
{
printf("不是有补格\n");
break;
}
}
if (ok)
printf("是有补格\n");
}
return 0;
}