POJ3090 Visible Lattice Points

题目

输入n，问从坐标(0,0)点能看到多少点，点的坐标范围是0<=(x,y)<=n
只有一个点没有被另一个点挡住才能看到。
如上图，n=5的时候，能看到21个点
Input
第一行输入t，表示样例个数。1<=t<=1000 接下来t行，每行输入一个正整数n 。 1<=n<=1000
Output
输出n行三列。 第一列从1到n表示第i个样例，第二列输出n，第三列输出能看到的点的个数
Sample Input
4
2
4
5
231
Sample Output
1 2 5
2 4 13
3 5 21
4 231 32549

我是这么想的，连接原点和可以看到的点，并延长，延长线上的点是看不到的，打表处理。

只画了一半，这个图关于主对角线对称
A点放在哪里都可以，为了打表方便，我把A点放在左上角，也就是把上图顺时针旋转90度

#include <stdio.h>
int vis[1010][1010]={0};
int main()
{
	for(int i=0;i<1001;i++)  //按行遍历
		for(int j=i+1;j<1001;j++)  //只遍历一半，所以从j=i+1开始（不包含主对角线，若想把主对角线包括进去从j=i开始，最后输出结果也要改）
			if(vis[i][j]==0)  //可到达的标为0
				for(int k=j+j,l=i+i;k<1001&&l<1001;k+=j,l+=i)  //找延长线上的点
				    vis[l][k]=1;  //不可到达的标为1
	int t,n,sum=0;
	scanf("%d",&t);
	for(int k=1;k<=t;k++)
	{
		sum=0;
		scanf("%d",&n);
		for(int i=0;i<=n;i++)
		    for(int j=i+1;j<=n;j++)
		        if(vis[i][j]==0)
		            sum++;
		printf("%d %d %d\n",k,n,sum*2+1);  //+1加的是主对角线上的点
	}
	return 0;
}