54 等差数列

最新推荐文章于 2024-03-03 10:50:03 发布

Dyanxier

最新推荐文章于 2024-03-03 10:50:03 发布

阅读量317

点赞数

分类专栏：东华大学OJ 文章标签：等差数列双平方数算法编程搜索

本文链接：https://blog.csdn.net/Dyanxier/article/details/123533295

版权

东华大学OJ 专栏收录该内容

11 篇文章

订阅专栏

这篇博客探讨了如何寻找双平方数集中长度为N的等差数列。作者通过编程实现了一个程序，该程序首先生成指定范围内的双平方数集，然后检查是否存在满足条件的等差数列。输出包括所有找到的等差数列，若未找到则输出`NONE`。示例输入输出展示了程序的运行效果。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

54 等差数列

作者: xxx 时间限制: 1S 章节: 一维数组

问题描述 :

一个等差数列是一个能表示成a, a+b, a+2b,..., a+nb (n=0,1,2,3,...) 在这个问题中a是一个非负的整数，b是正整数。

写一个程序来找出在双平方数集合S中长度为n的等差数列。双平方数集合是所有能表示成p2＋q2的数的集合。

输入说明 :

第一行: N(3<= N<=25),要找的等差数列的长度。第二行: M(1<= M<=250),搜索双平方数的上界0 <= p,q <= M。

输出说明 :

如果没有找到数列,输出`NONE'。

如果找到了，输出一行或多行, 每行由两个整数组成:a,b 这些行应该先按b排序再按a排序(均为升序)。

将不会有多于10,000个等差数列。

输入范例 :

5 7

输出范例 :

1 4
37 4
2 8
29 8
1 12
5 12
13 12
17 12
5 20
2 24

#include<stdio.h>
#include<string.h>

int N,M;
int set[125000];//双平方数集
int sup;//上界
int ans;//符合要求的a和b的对数

void initSet(){
    //获取0<=p,q<=M的双平方数集
	int pos;
	for(int i=0;i<=M;i++){
		for(int j=i;j<=M;j++){
			pos=i*i+j*j;
			set[pos]=1;//标记该下标为双平方数
		}
	}
	sup=M*M+M*M;
}

bool isOk(int n){
    //验证整数n是否为确定范围之内的双平方数
	if(set[n]==1)return true;
	return false;
}

void findOut(){
	ans=0;
    int a,b,i;
	for(b=1;b*(N-1)<=sup;b++){
		//b在外层，a在内层
		//使得结果先按b排序再按a排序(均为升序)
		for(a=0;a+b*(N-1)<=sup;a++){
			for(i=0;i<N;i++){
				if(!isOk(a+b*i))break;
			}
			if(i==N){
				printf("%d %d\n",a,b);
				ans++;
			}
		}
	}
}

int main(){

    memset(set,0,sizeof(set));//初始化set数组
	scanf("%d%d",&N,&M);
	initSet();
    findOut();
	if(ans==0)printf("NONE\n");
	return 0;

}