# include<stdio.h>
# include<math.h>
# define N 12
int main(void)
{
float a[3][N] = {{0},{0},{0.50,1.41,3.00,3.37,3.40,4.72,4.72,5.43,7.57,8.38,8.98,9.50}};
float b[3][N] = {{0},{0},{2.00,3.50,1.50,3.49,5.50,2.00,6.24,4.10,2.01,4.50,3.41,0.80}};
float t;
int i,j,k,l,m,n,an=0,bn=0,maxnum=0;
int a_idx[N][N]={0},b_idx[N][N]={0};
int anum[N]={0},bnum[N]={0},max[N]={0},tmp[N]={0};
//步骤1：初始化坐标数组，并处理坐标值
for(i=0; i<N; i++)
{
a[0][i] = i+1;
t = a[2][i]-(int)a[2][i];
if(t > 0.5)
t -= 1;
a[1][i] = t;
b[0][i] = i+1;
t = b[2][i]-(int)b[2][i];
if(t > 0.5)
t -= 1;
b[1][i] = t;
}
//步骤2：排序
for(i=0; i<N; i++)
{
for(j=0; j<N-1-i; j++)
{
if(a[1][j] > a[1][j+1])
for(l=0; l<3; l++)
t=a[l][j], a[l][j]=a[l][j+1], a[l][j+1]=t;
if(b[1][j] > b[1][j+1])
for(l=0; l<3; l++)
t=b[l][j], b[l][j]=b[l][j+1], b[l][j+1]=t;
}
}
//步骤3：分子组
for(i=0,k=0; i<N; i++,k++)
{
t = a[1][i]+0.1;
for(j=i,l=0; j<N&&a[1][j]<=t; j++,l++)
{
a_idx[k][l] = a[0][j];
anum[k] = l;
}
a_idx[k][l] = 0;
an = k;
if(k!=0 && a_idx[k][anum[k]]==a_idx[k-1][anum[k-1]])
{
for(n=0; n<N; n++)
a_idx[k][n] = 0;
k -= 1;
}
}
for(i=0,k=0; i<N&&k<N; i++,k++)
{
t = b[1][i]+0.1;
for(j=i,l=0; j<N&&b[1][j]<=t; j++,l++)
{
b_idx[k][l] = b[0][j];
bnum[k] = l;
}
b_idx[k][l] = 0;
bn = k;
if(k!=0 && b_idx[k][bnum[k]]==b_idx[k-1][bnum[k-1]])
{
for(n=0; n<N; n++)
b_idx[k][n] = 0;
k -= 1;
}
}
//步骤4：求交集
for(i=0; i<an; i++)
for(k=0; k<bn; k++)
{
n=0;
for(j=0; j<=anum[i]; j++)
{
for(l=0; l<=bnum[k]; l++)
{
if(a_idx[i][j] == b_idx[k][l])
{
tmp[n]=a_idx[i][j];
n++;
break;
}
}
}
if(n>=maxnum)
{
maxnum=n;
for(m=0; m<maxnum; m++)
max[m]=tmp[m];
}
}
//步骤5：输出
//输出二维数组a,b
printf("\nA:\n");
for(i=0; i<3; i++)
{
for(j=0; j<N; j++)
if(i==0)
printf("%6.0f",a[i][j]);
else
printf("%6.2f",a[i][j]);
printf("\n");
}
printf("\nB:\n");
for(i=0; i<3; i++)
{
for(j=0; j<N; j++)
if(i==0)
printf("%6.0f",b[i][j]);
else
printf("%6.2f",b[i][j]);
printf("\n");
}
//输出子组
printf("\nAi:\n");
for(i=0; i<an; i++)
{
for(j=0; j<=anum[i]; j++)
printf("%6d", a_idx[i][j]);
printf("\n");
}
printf("\nBj:\n");
for(i=0; i<bn; i++)
{
for(j=0; j<=bnum[i]; j++)
printf("%6d", b_idx[i][j]);
printf("\n");
}
//输出最大交集及其对应坐标
printf("\nAi∩Bj:\n");
for(i=0; i<maxnum; i++)
printf("%8d\t", max[i]);
printf("\n");
for(i=0; i<maxnum; i++)
{
for(j=0; j<N; j++)
{
if(abs(a[0][j]-max[i]) < 1e-3)
printf("(%.2f,", a[1][j]);
if(abs(b[0][j]-max[i]) < 1e-3)
printf("%.2f)\t", b[1][j]);
}
}
printf("\n");
return 0;
}