【计蒜客系列】挑战难题26：等差数列

题目来源：计蒜客

一个等差数列是一个能表示成a, a+b, a+2b,..., a+nb (n=0,1,2,3,...)的数列。
在这个问题中a是一个非负的整数，b是正整数。写一个程序来找出在双平方数集合(双平方数集合是所有能表示成p2＋q2的数的集合)S中长度为n的等差数列。
输入包括两行，第一行为N要找的等差数列的长度（3≤ N≤25）。第二行是找到的双平方数p和q的上界M（0≤p,q≤M）。
输出一行或者多行，如果没有找到数列,输出NONE。否则输出一个整数对a b（这些行应该先按b排序再按a排序）
样例1
输入：
5
7
输出：
1 4
37 4
2 8
29 8
1 12
5 12
13 12
17 12
5 20
2 24

注：本题确实没有看懂，代码参考网上的，在此贴出来，大家看看就好。不喜勿喷。

#include<iostream>
#include<cstring>
using namespace std;
bool is[125001];
int a[125001];
int main() {

	int n, m;
	scanf("%d%d", &n, &m);
	for (int i = 0; i <= m; ++i)
		for (int j = 0; j <= m; ++j)
			is[i * i + j * j] = 1;
	int k = 0;
	for (int i = 0, z = m * m * 2; i <= z; ++i)
		if (is[i])
			a[k++] = i;
	int l = 1, ma = (m * m * 2 - 0) / (n - 1);
	int nu = 0;
	for (; l <= ma; ++l) {
		for (int i = 0; a[i] + (n - 1) * l <= m * m * 2 && i < k; ++i) {
			bool f = 1;
			for (int j = 1; j < n && f; ++j)
				if (!is[a[i] + j * l]) {
					f = 0;
				}
			if (f) {
				printf("%d %d\n", a[i], l);
				++nu;
			}
		}
	}
	if (!nu)
		puts("NONE");
	return 0;
}