码迷,mamicode.com
首页 > 其他好文 > 详细

【vijos】1543 极值问题(数论+fib数)

时间:2014-10-07 13:30:53      阅读:208      评论:0      收藏:0      [点我收藏+]

标签:blog   http   io   os   ar   for   sp   div   问题   

https://vijos.org/p/1543

好神奇的一题。。

首先我竟然忘记n可以求根求出来,sad。

然后我打了表也发现n和m是fib数。。

严格证明(鬼知道为什么这样就能对啊,能代换怎么就能保证最大呢?):

(n^2-mn-m^2)^2=1

(m^2+mn-n^2)^2=1

(m(m+n)-n^2)^2=1

((m+n-n)(m+n)-n^2)^2=1

((m+n)^2-n(m+n)-n^2)^2=1

因为(n‘^2-m‘n‘-m‘^2)^2=1

取n‘=m+n, m‘=n使式子成立。。

然后就是暴力就行了

#include <cstdio>
#include <cstring>
#include <cmath>
#include <string>
#include <iostream>
#include <algorithm>
#include <queue>
using namespace std;
typedef long long ll;
#define rep(i, n) for(int i=0; i<(n); ++i)
#define for1(i,a,n) for(int i=(a);i<=(n);++i)
#define for2(i,a,n) for(int i=(a);i<(n);++i)
#define for3(i,a,n) for(int i=(a);i>=(n);--i)
#define for4(i,a,n) for(int i=(a);i>(n);--i)
#define CC(i,a) memset(i,a,sizeof(i))
#define read(a) a=getint()
#define print(a) printf("%d", a)
#define dbg(x) cout << (#x) << " = " << (x) << endl
#define printarr2(a, b, c) for1(_, 1, b) { for1(__, 1, c) cout << a[_][__]; cout << endl; }
#define printarr1(a, b) for1(_, 1, b) cout << a[_] << ‘\t‘; cout << endl
inline const int getint() { int r=0, k=1; char c=getchar(); for(; c<‘0‘||c>‘9‘; c=getchar()) if(c==‘-‘) k=-1; for(; c>=‘0‘&&c<=‘9‘; c=getchar()) r=r*10+c-‘0‘; return k*r; }
inline const ll max(const ll &a, const ll &b) { return a>b?a:b; }
inline const int min(const int &a, const int &b) { return a<b?a:b; }

ll k, a, b, ans;

int main() {
	read(k);
	ll n=1, m=0;
	while(n+m<=k) {
		ll c=m;
		m=n;
		n+=c;
	}
	printf("%lld\n", m*m+n*n);
	return 0;
}

 

 


 

 

背景

小铭的数学之旅2。

描述

已知m、n为整数,且满足下列两个条件:
① m、n∈1,2,…,K
② (n^ 2-mn-m^2)^2=1
编一程序,对给定K,求一组满足上述两个条件的m、n,并且使m^2+n^2的值最大。例如,若K=1995,则m=987,n=1597,则m、n满足条件,且可使m^2+n^2的值最大。

格式

输入格式

输入仅一行,K的值。

输出格式

输出仅一行,m^2+n^2的值。

样例1

样例输入1[复制]

1995

样例输出1[复制]

3524578

限制

每个测试点1秒。

提示

Source:
汕头市FXOI组
Phoeagon
ThanX2 Sivon
For TripleY

【vijos】1543 极值问题(数论+fib数)

标签:blog   http   io   os   ar   for   sp   div   问题   

原文地址:http://www.cnblogs.com/iwtwiioi/p/4009100.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!