poj 1320 Street Numbers(解pell方程)
2014-03-08 15:52
483 查看
题意:求满足1+2+3+……+(n-1)=(n+1)+(n+2)+……+m的前10项n和m。
思路:把公式化简一下,n(n-1)/2=(m-n)(m+n+1)/2,然后可以求出2n^2=m^2+m,给m配方可以得到(2m+1)^2-8n^2=1。设x=2m+1,y=n,这样就变成了x^2-dy^2=1的形式,显然这是个pell 方程。利用迭代公式就可以求解了:
Xn=Xn-1 * X1 +dYn-1 *Y1
Yn=Xn-1 * Y1 + Yn-1 * X1
比较令人无语的是我最开始是看书上的,书上翻译就错了。。。然后我公式死活推不出来,上网上翻了一下题解,发现好多题解的公式中间那一步都不对。。。最后又看了看题才发现问题,太坑了
代码:
思路:把公式化简一下,n(n-1)/2=(m-n)(m+n+1)/2,然后可以求出2n^2=m^2+m,给m配方可以得到(2m+1)^2-8n^2=1。设x=2m+1,y=n,这样就变成了x^2-dy^2=1的形式,显然这是个pell 方程。利用迭代公式就可以求解了:
Xn=Xn-1 * X1 +dYn-1 *Y1
Yn=Xn-1 * Y1 + Yn-1 * X1
比较令人无语的是我最开始是看书上的,书上翻译就错了。。。然后我公式死活推不出来,上网上翻了一下题解,发现好多题解的公式中间那一步都不对。。。最后又看了看题才发现问题,太坑了
代码:
#include<iostream> #include<cstdio> #include<cstring> #include<string> #include<algorithm> #include<map> #include<queue> #include<stack> #include<cmath> #include<vector> #define inf 0x3f3f3f3f #define Inf 0x3FFFFFFFFFFFFFFFLL #define eps 1e-9 #define pi acos(-1.0) using namespace std; typedef long long ll; int main() { //freopen("in.txt","r",stdin); //freopen("out.txt","w",stdout); int xn=3,yn=1; for(int i=0;i<10;++i) { int tmp1=xn*3+8*yn; int tmp2=xn+yn*3; printf("%10d%10d\n",tmp2,(tmp1-1)/2); xn=tmp1;yn=tmp2; } return 0; }
相关文章推荐
- POJ 1320-Street Numbers(解佩尔方程)
- POJ 1320 Street Numbers 佩尔方程
- poj 1320 Street Numbers 解pell方程
- POJ - 1320 - Street Numbers - (佩尔方程)
- POJ 1320:Street Numbers——佩尔方程
- POJ 1320 Street Numbers 笔记
- POJ 1320 Street Numbers [佩尔方程]
- POJ 1320Street Numbers(佩尔方程定理)
- POJ1320 Street Numbers【佩尔方程】
- poj 1320——Street Numbers
- POJ1320 Street Numbers(佩尔方程)
- POJ1320 UVA138 Street Numbers【Pell方程+数学】
- POJ 1320 Street Numbers 解佩尔方程
- UVa 138 Street Numbers(数论&Pell方程)
- POJ1320 Street Numbers【佩尔方程】
- POJ 2115 C Looooops (模线性方程)
- poj_2115C Looooops(模线性方程)
- POJ 1186 方程的解数
- POJ 2115 C Looooops 解题报告(模线性方程)
- POJ 1186 方程的解数 [解题报告] Java