首页 > 代码库 > 【vijos】1543 极值问题(数论+fib数)
【vijos】1543 极值问题(数论+fib数)
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 << endlinline 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数)
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。