题意简述

\(n\)是正整数,其四个最小的因子分别为\(d_1,d_2,d_3,d_4\)。

求对于所有的\(n \le m\)满足

\[d_1^2+d_2^2+d_3^2+d_4^2=n
\]

的\(n\)的个数。

分析

这是一道数竞题,下面来求\(n\)。

显然\(d_1=1\),则\(n=1+d_2^2+d_3^2+d_4^2\)。

若\(n\)是奇数,则其因子都是奇数,其平方和为四个奇数相加为偶数,矛盾。故n为偶数,\(d_2=2\),则\(n-5=d_3^2+d_4^2\)。

此时,\(d_3\)只可能是\(3,4\)或者大于\(4\)的质数

  • 若\(d_3=3\),则\(d_4=4\)或\(d_4=6\),验算得不存在这样的\(n\)
  • 若\(d_3=4\),则由奇偶性可得,\(d_4\)一定为奇数,所以\(d_4^2 \equiv 1(\mbox{mod } 4)\),但\(n-5-d_3^2=n-21 \equiv 3(\mbox{mod }4)\),所以矛盾。

所以\(n\)只能是一个大于\(4\)的质数,且\(d_4=2d_3\),故\(n-5=5d_3^2 \Longrightarrow n=5(d_3^2+1)\)。故\(5 \nmid n\),即\(d_3=5,d_4=10\),得\(n=130\)。故\(\forall n \in \mathbb{N^+}\)只有\(n=130\)符合题意。

Code

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<map>
#include<set>
#include<queue>
#include<vector>
#define IL inline
#define re register
#define LL long long
#define ULL unsigned long long
#define re register
#define debug printf("Now is %d\n",__LINE__);
using namespace std; template<class T>inline void read(T&x)
{
char ch=getchar();
while(!isdigit(ch))ch=getchar();
x=ch-'0';ch=getchar();
while(isdigit(ch)){x=x*10+ch-'0';ch=getchar();}
}
inline LL read()
{
int x=0;
char ch=getchar();
while(!isdigit(ch))ch=getchar();
x=ch-'0';ch=getchar();
while(isdigit(ch)){x=x*10+ch-'0';ch=getchar();}
return x;
}
int G[55];
template<class T>inline void write(T x)
{
int g=0;
if(x<0) x=-x,putchar('-');
do{G[++g]=x%10;x/=10;}while(x);
for(re int i=g;i>=1;--i)putchar('0'+G[i]);putchar('\n');
}
LL n;
int main()
{
n=read();
if(n>=130) cout<<1<<endl;
else cout<<0<<endl;
return 0;
}

最新文章

  1. js判断本机是否已安装app
  2. em与px换算关系以及常用列表
  3. 获取浏览器信息 c#
  4. CentOS6.3修复模式/单用户模式修改fstab文件
  5. 搭建DAO层和Service层代码
  6. 修改eclipse启动时eclipse使用的jre
  7. 启动C:\Windows\System32\logiLDA.DLL时出现问题,找不到指定模块
  8. [wikioi]能量项链
  9. CSS3绘制环形进度条
  10. 三 ICE开发初级研究
  11. 【转】AB实验设计思路及实验落地
  12. SparkStreaming+Kafka整合
  13. lua 5.3最简单plugin编写
  14. Jmeter 结构、原理介绍
  15. ELK篇---------elasticsearch集群安装配置
  16. Windows环境下,从零开始搭建Nodejs+Express+Ejs框架(二)---安装Express,ejs
  17. PHP array与Json学习
  18. zoj3696(泊松分布)
  19. opencv画图
  20. PIVOT 和 UNPIVOT 命令的SQL Server版本

热门文章

  1. 日志分析平台ELK之日志收集器logstash常用插件配置
  2. 在IDEA创建类时自动创建作者日期等信息设定
  3. 07 Sublime Text3常用快捷键
  4. Unicode和多字节字符集
  5. 【题解】小Z的袜子
  6. ElasticSearch 索引 VS MySQL 索引
  7. golang拾遗:指针和接口
  8. Elasticsearch(4):映射
  9. go 继承
  10. linux(centos8):用grep命令查找文件内容