1243: CKJ老师爱数学

时间限制: 1 Sec  内存限制: 128 MB
提交: 56  解决: 13
[提交][状态][讨论版]

题目描述

众所周知,CKJ老师非常热爱数学,他对于方程组的有自己的独到的研究,today,他抛给了你一个too simple的方程组x^2+y^2=z^2,z是一个已给的正整数,然后他毫不客气地问你,这个方程组的整数解有几个?

输入

包含一个整数T表示输入数据组数,

接下来T行每行一个整数z(z<=2000 000 000)

输出

每行一个整数表示方程组的解的个数

样例输入

1
4

样例输出

4
---------------------------------------------------------------------------------------------

自认为这是xdoj上思考最多的一道题: 网上好像有这样题的公式,但做题最要的不是理解学习心得知识吗。

首先我贴出一些知识点:

1) 定理 1   x2 + y2 = z2 的互質整數解為 x=2uvy=u2-v2, 其中 uv 是任意互質整數,且 uv 不同時為奇數。

證明 1 令 xyz 是互質正整數,且 x2 + y2 = z2。若 x 與 y 都是奇數,則 z 是偶數。故 x2 + y2 是 4l+2 的型式,z2 是 4l' 的型式(l 與 l' 是整數)。矛盾。

因此,不妨假定 x 是偶數,y 與 z 是奇數。

由 x2 = z2 - y2 = (z+y) (z-y),令 x=2rz-y = 2sz+y = 2t。得 r2 = st

s 與 t 互質,因 y 與 z 互質。但 st 是完全平方。故 s=u2 , t=v2。得證。

网址链接:http://episte.math.ntu.edu.tw/articles/mm/mm_07_4_01/page3.html

2)定理二:费马平方和定理    奇质数能表示为两个平方数之和的充分必要条件是该素数被4除余1。

3)定理三:  (a^2+b^2)*(c^2+d^2)==(ac+bd)^2+(ad-bc)^2

网址链接: https://zh.wikipedia.org/wiki/%E8%B4%B9%E9%A9%AC%E5%B9%B3%E6%96%B9%E5%92%8C%E5%AE%9A%E7%90%86

问题解决:

 1 借助定理1 将x^2+y^2=z^2  转化为 x^2+y^2=n 有多少整数解

 2 每个整数解我们可以看出是由本原解组成(三个数互质)【 3,4,5】是,而【6,8,10】不是

3 z=p1^k1*p2^k2*....pn^kn      

根据定理二  对于每个素数pk且满足pk%4==1   一个本原解  a^2+b^2==pk

对于每一个素数要么可以充当本原解的组成,要么组成解的因子  

所以对于每一个因子可以有(2*sum+1) 种选择

eg (z==25   sum=2  有五种选择

     1+4=5;  解得x=5^2*(2*u*v)=4*25=100; y=25*(u^2-v*2)=75;【一个5是因子,一个5组成本原解】

     3^2+4^2=25  解得y=2uv=24; x=u*2-v*2=7;                         【两个5都是本原解】

    两个解调换位置 4个解

   25全部作为因子  1个解

4  根据定理三  两个本原解可以合并乘一个本原解  所以不同因子之间是相乘关系

一共有   ans=∏ (2*pk+1)-1  本原解  【 减一是排除所有素数作为因子】

考虑正负号   ans*=4

再加上四个点  sum=ans+4;         【端点我们不认为是本原解】

 #include<cstdio>
#include<algorithm>
#include<iostream>
using namespace std;
typedef long long LL;
const LL N=1e6+;
LL p[N+];
int main ()
{
for (LL i=;i<=N;i++) {
if (!p[i]) p[++p[]]=i;
for (LL j=;j<=p[]&&i*p[j]<=N;j++) {
p[i*p[j]]=;
if (i%p[j]==) break;
}
}
int T;
scanf ("%d",&T);
while (T--) {
LL x;
scanf("%lld",&x);
LL i=;
LL ans=;
while (p[i]*p[i]<=x) {
LL sum=;
if (x%p[i]==) {
while (x%p[i]==) {
sum++;
x=x/p[i];
}
if (p[i]%==)
ans=ans*(*sum+);
}
i++;
}
if(x!=&&(x%==)) ans=ans*;
printf("%lld\n",(ans-)*+);
}
return ;
}

这完全是我自己的理解,可能存在说的不清晰地方或者不严谨,请

多多指教;

xdoj-1243 (费马平方和问题)的更多相关文章

  1. BZOJ 1041 [HAOI2008]圆上的整点:数学【费马平方和定理】

    题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=1041 题意: 给定n(n <= 2*10^9),问你在圆x^2 + y^2 = n^ ...

  2. 费马平方和定理&&斐波那契恒等式&&欧拉四平方和恒等式&&拉格朗日四平方和定理

    费马平方和定理 费马平方和定理的表述是:奇素数能表示为两个平方数之和的充分必要条件是该素数被4除余1. 1. 如果两个整数都能表示为两个平方数之和的形式,则他们的积也能表示为两个平方数之和的形式. $ ...

  3. hdu 4704 Sum (整数和分解+快速幂+费马小定理降幂)

    题意: 给n(1<n<),求(s1+s2+s3+...+sn)mod(1e9+7).其中si表示n由i个数相加而成的种数,如n=4,则s1=1,s2=3.                  ...

  4. nyoj1000_快速幂_费马小定理

    又见斐波那契数列 时间限制:1000 ms  |  内存限制:65535 KB 难度:4   描述 斐波那契数列大家应该很熟悉了吧.下面给大家引入一种新的斐波那契数列:M斐波那契数列. M斐波那契数列 ...

  5. poj 3734 Blocks 快速幂+费马小定理+组合数学

    题目链接 题意:有一排砖,可以染红蓝绿黄四种不同的颜色,要求红和绿两种颜色砖的个数都是偶数,问一共有多少种方案,结果对10007取余. 题解:刚看这道题第一感觉是组合数学,正向推了一会还没等推出来队友 ...

  6. 数论初步(费马小定理) - Happy 2004

    Description Consider a positive integer X,and let S be the sum of all positive integer divisors of 2 ...

  7. 【BZOJ1951】【SDOI2010】古代猪文 Lucas定理、中国剩余定理、exgcd、费马小定理

    Description “在那山的那边海的那边有一群小肥猪.他们活泼又聪明,他们调皮又灵敏.他们自由自在生活在那绿色的大草坪,他们善良勇敢相互都关心……” ——选自猪王国民歌 很久很久以前,在山的那边 ...

  8. 数论 --- 费马小定理 + 快速幂 HDU 4704 Sum

    Sum Problem's Link:   http://acm.hdu.edu.cn/showproblem.php?pid=4704 Mean: 给定一个大整数N,求1到N中每个数的因式分解个数的 ...

  9. HDU 5667 Sequence 矩阵快速幂+费马小定理

    题目不难懂.式子是一个递推式,并且不难发现f[n]都是a的整数次幂.(f[1]=a0;f[2]=ab;f[3]=ab*f[2]c*f[1]...) 我们先只看指数部分,设h[n]. 则 h[1]=0; ...

随机推荐

  1. GCC 命令行详解 -L 指定库的路径 -l 指定需连接的库名 zhuan

    1.gcc包含的c/c++编译器gcc,cc,c++,g++,gcc和cc是一样的,c++和g++是一样的,(没有看太明白前面这半句是什么意思:))一般c程序就用gcc编译,c++程序就用g++编译 ...

  2. zabbix3.4.7监控linux进程

    利用zabbix proc.num方法监控Linux服务进程 proc.num[<name>,<user>,<state>,<cmdline>] 监控用 ...

  3. 普通01背包问题(dp)

    有n个物品,重量和价值分别为wi和vi,从这些物品中挑选出重量不超过W的物品,求所有挑选方案中物品价值总和的最大值 限制条件: 1 <= n <= 100; 1 <= wi,vi & ...

  4. 逆袭之旅DAY20.XIA.程序调试

    2018-07-16 20:25:50 F5:进入方法 F6:单步执行

  5. Win10系列:UWP界面布局基础12

    画刷 画刷(Brush)用于为图形元素填充颜色.在XAML中,画刷有许多属性,其中较常使用的是Fill属性和Stroke属性,Fill用于填充图形的背景色,而Stroke用于设置图形的线条颜色. 在实 ...

  6. POJ 1840 Eqs 解方程式, 水题 难度:0

    题目 http://poj.org/problem?id=1840 题意 给 与数组a[5],其中-50<=a[i]<=50,0<=i<5,求有多少组不同的x[5],使得a[0 ...

  7. HTML5 ②

    html代码中的各种标签和特殊符号: 1.<p>段落标签 默认会在段落前后都有空行</p>: 2.<hr/> 水平线标签     <br/>换行标签   ...

  8. Saiku缓存处理(七)

    Saiku缓存处理方案 Saiku默认是从缓存中读取数据的(如果缓存中有数据的话),所以用户看到的数据不一定是最新的,如果需要看到最新的的数据需要手动刷新数据或者更改配置信息. Saiku获取实时数据 ...

  9. day17 面向对象-成员

    今日主要学习了: 1 .成员 2. 变量 3. 方法 4.属性 5.私有 1. 成员 在类中你能写的所有内容都是成员 2 .变量 1) 实例变量: 昨天写的就是实例变量 ,由对象去访问的变量 . cl ...

  10. Python中os与sys模块的区别

    os与sys模块的官方解释如下: os: This module provides a portable way of using operating system dependent functio ...