NOIP2008普及组 第二题 排座椅
原题 试题描述上课的时候总有一些同学和前后左右的人交头接耳，这是令小学班主任十分头疼的一件事情。不过，班主任小雪发现了一些有趣的现象，当同学们的座次确定下来之后，只有有限的D对同学上课时会交头接耳。同学们在教室中坐成了M行N列，坐在第i行第j列的同学的位置是（i，j），为了方便同学们进出，在教室中设置了K条横向的通道，L条纵向的通道。于是，聪明的小雪想到了一个办法，或许可以减少上课时学生交头接耳的问题：她打算重新摆放桌椅，改变同学们桌椅间通道的位置，因为如果一条通道隔开了两个会交头接耳的同学，那么他们就不会交头接耳了。
请你帮忙给小雪编写一个程序，给出最好的通道划分方案。在该方案下，上课时交头接耳的学生对数最少。
输入
输入文件seat.in的第一行，有5各用空格隔开的整数，分别是M，N，K，L，D（2<=N，M<=1000，0<=K，0<=L，D<=2000）。
接下来D行，每行有4个用空格隔开的整数，第i行的4个整数Xi，Yi，Pi，Qi，表示坐在位置(Xi,Yi)与(Pi,Qi)的两个同学会交头接耳（输入保证他们前后相邻或者左右相邻）。
输入数据保证最优方案的唯一性。
输出
输出文件seat.out共两行。
第一行包含K个整数，a1a2……aK，表示第a1行和a1+1行之间、第a2行和第a2+1行之间、…、第aK行和第aK+1行之间要开辟通道，其中ai< ai+1，每两个整数之间用空格隔开（行尾没有空格）。
第二行包含L个整数，b1b2……bk，表示第b1列和b1+1列之间、第b2列和第b2+1列之间、…、第bL列和第bL+1列之间要开辟通道，其中bi< bi+1，每两个整数之间用空格隔开（行尾没有空格）。
输入样例
4 5 1 2 3
4 2 4 3
2 3 3 3
2 5 2 4
输出样例
2
2 4
输入/输出文件基本名
seat
我的代码是：
程序代码：#include <iostream>
#include <fstream>
#include <math.h>
using namespace std;
long int a[201][201]={0};
int b[2][201]={0},c[2][201]={0},e[6600]={0};
int main(){
ifstream fin("seat.in");
ofstream fout("seat.out");
int m,n,k,l,d,i,j,r=2,t,p,q;
fin >> m >> n >> k >> l >> d ;
e[1]=2;
for (i=3; i<65536; i++){
float s=sqrt(i);
for (j=1 ; e[j]<=s; j++)
if(i%e[j]==0)break;
if (e[j]>s){
e[r]=i;
r++;
}
}
for (r=1; r<=d; r++)
for (t=1; t<3; t++){
fin >> i >> j;
if (t==1) { p=i; q=j; }
if (a[i][j]!=0){
if (t==1) { a[i][j]=e[r]*a[i][j]; continue ;}
else { a[i][j]=a[p][q]*a[i][j]; continue ;}
}
a[i][j]=e[r];
}
for (i=1; i<m; i++){
r=0;
for (j=1; j<=n; j++) if (a[i][j]!=0&&a[i+1][j]!=0&&(a[i][j]%a[i+1][j]==0||a[i+1][j]%a[i][j]==0)) r++;
b[0][i]=i; b[1][i]=r;
}
for (i=1; i<m-1; i++)
for (j=1; j<m-i; j++)
if (b[1][j]<b[1][j+1])
{ p=b[0][j]; b[0][j]=b[0][j+1]; b[0][j+1]=p;
q=b[1][j]; b[1][j]=b[1][j+1]; b[1][j+1]=q; }
for (i=1; i<k; i++)
for (j=1; j<=k-i; j++)
if (b[0][j]>b[0][j+1])
{ p=b[0][j]; b[0][j]=b[0][j+1]; b[0][j+1]=p;
q=b[1][j]; b[1][j]=b[1][j+1]; b[1][j+1]=q; }
fout << b[0][1];
for (i=2; i<=k; i++) fout << " " << b[0][i];
fout << endl;
for (j=1; j<n; j++){
r=0;
for (i=1; i<=m; i++) if (a[i][j]!=0&&a[i][j+1]!=0&&(a[i][j]%a[i][j+1]==0||a[i][j+1]%a[i][j]==0)) r++;
c[0][j]=j; c[1][j]=r;
}
for (i=1; i<n-1; i++)
for (j=1; j<n-i; j++)
if (c[1][j]<c[1][j+1])
{ p=c[0][j]; c[0][j]=c[0][j+1]; c[0][j+1]=p;
q=c[1][j]; c[1][j]=c[1][j+1]; c[1][j+1]=q; }
for (i=1; i<l; i++)
for (j=1; j<=l-i; j++)
if (c[0][j]>c[0][j+1])
{ p=c[0][j]; c[0][j]=c[0][j+1]; c[0][j+1]=p;
q=c[1][j]; c[1][j]=c[1][j+1]; c[1][j+1]=q; }
fout << c[0][1];
for (i=2; i<=l; i++) fout << " " << c[0][i];
fout << endl;
fin.close();fout.close();
return 0;
}
有几组比较大的测试数据总是不能通过，但是我找了很多遍也找不到哪里错了。请大家帮帮忙，谢谢了！
