斐波那契首项 BAPC2014 I&&HUNNU11589:Interesting Integers

原创 2015年11月21日 16:09:21

题目大意:给定一个数n,按斐波那契数规则生成,求最小的b,并使a最小。G1 = a , G2 = b; b >= a.

思路分析:

因为n<10^9,所以斐波那契数列的项数不可能超过46.题目转换为求ax+by=n的最小的x,y。根据欧几里得

#include <iostream>
#include <cmath>
#include <string>
#include <cstring>
#include <algorithm>
#include <cstdio>

using namespace std;

#define maxn 99999999
#define mem(a) memset(a , 0 , sizeof(a))
#define LL __int64
LL fib[65];

void init()
{
    fib[0] = 1;
    fib[1] = 1;
    fib[2] = 2;
    for(int i = 3 ; i < 65 ; i ++)
    {
        fib[i] = fib[i-1] + fib[i-2];
    }
    return ;
}

void gcd2(LL a , LL b , int &d ,LL &x , LL &y)
{
    if(!b)
    {
        d = a;
        x = 1 ;
        y = 0;
    }
    else
    {
        gcd2(b , a%b , d , y , x);
        y -= x*(a/b);
    }
}


int main()
{
    int t;
    LL n;
    init();
    scanf("%d" , &t);
    while(t-- )
    {
        scanf("%I64d" , &n);
        int pos;
        pos = 65;
        LL ansx = (LL)maxn*1000, ansy = (LL)maxn*1000;
        LL x , y ;int d;
        while(pos > 2)
        {
            LL a = fib[pos - 2], b = fib[pos - 1];
            pos -- ;
            gcd2(a , b , d , x , y);
            x *= n/d ;
            y *= n/d;
            LL t1 = b / d;
            LL t2 = a / d;
            x = (x % t1 + t1) % t1;     //找到最小正整数x
            y = (n - a * x) / b;     //对应的y;
            
            ////////////并没有看懂这个
             ll   k = (y-x-1) / (ma+mb);//优化  
            x += ma*k, y -= mb*k;  
           /* int abc;
            if(x <= 0)
            {
                abc = abs(x / b);
                if(x % b != 0)  abc ++;
                x = x + abc * b / d;
                y = y - abc * a / d;
            }
            if(y <= 0)
            {
                abc = abs(y / a);
                if(y % a != 0)  abc ++;
                x = x - abc * b / d;
                y = y + abc * a / d;
            }
            if(x < 0 || y + x < 0) continue;
            if( x == 0) swap(x , y);
            else if(y != 0 && x > y) swap(x , y);
            if(x == 0 && y == 0) continue;*/
            while(y >= x && y < ansy)
            {
                if(y < ansy || (y == ansy && x < ansx) && x > 0) ansx = x , ansy = y;
                x += t1  , y -= t2;
            }
        }
        if(ansx == 0) ansx = ansy ;
        printf("%I64d %I64d\n" , ansx ,ansy);
    }
}


版权声明:本文为博主原创文章,未经博主允许不得转载。

相关文章推荐

BAPC2014 I&&HUNNU11589:Interesting Integers

题意: 对于我们已知的斐波那契数列,现在给出一个n,要我们求出一个新的斐波那契数列起始项,使得n能在新的斐波那契数列中,要求起始项y最小 思路: 我们知道 a3 = a1+a2 a4 = a1+2...

Interesting Integers(CF---BAPC 14 + hnoj11589)扩展欧几里得

I Interesting Integers source: xkcd.com/587 Undoubtedly you know of the Fibonacci numbers. Starting ...

BAPC2014 B&&HUNNU11582:Button Bashing(BFS)

题意: 给出n,m,代表微波炉有n个按钮,要求达到总时间为m 然后给出n个数,代表n个按钮能增加的时间,问最少几步,能够使得按出的总时间大于等于要求的时间,并且相差最小 输出最小的步数与相差的最小值 ...

BAPC2014 E&&HUNNU11585:Excellent Engineers(线段树)

题意: 有n个程序员,每个程序员有三个能力值,当一个程序员三个能力值都小于另外一个程序员的时候,那么就可以淘汰那个程序员,问经过淘汰之后,还剩下几个程序员。 思路: 我们首先按照x排序,那么我们就...

BAPC2014 C&&HUNNU11583:Citadel Construction(几何)

题意: 给出一系列的点,要求寻找最多4个点,使得组成一个面积最大的多边形 思路: 很显然只有两种情况,要么是三角形,要么是四边形 首先不难想到的是,先要把最外面的点都找出来,其实就是找凸包 但是并...

2014蓝桥杯本科C/C++组预赛第9题<斐波那契>解题报告

首先是题目链接 http://lx.lanqiao.org/problem.page?gpid=T121 如有错误,欢迎指出,讨论~ 蓝桥杯里面的数据相当水...很容易就能够水过去 现场思路 ...
  • ronnoc
  • ronnoc
  • 2014-03-26 20:38
  • 3486

du熊学斐波那契I

// 11号的入门题,手机收到邀请,但没时间参加,加班去了,现在补上!/* du熊学斐波那契I Time Limit : 2000/1000ms (C/Other) M...

【中途相遇法】【STL】BAPC2014 K Key to Knowledge (Codeforces GYM 100526)

题目链接:   http://codeforces.com/gym/100526   http://acm.hunnu.edu.cn/online/?action=problem&type=sho...

BAPC2014 K题 Key to Knowledge(中途相遇法,hash,状态压缩)

题意:n个学生和m道选择题,每个学生有关于m个选择题的答案,表示为一个字符串,1位√ 0位X,后面给一个老师给的对了几个答案 现在给你n个学生的答案,问,一共有多少种可能的答案,如果只有一个答案输出...

2014 Benelux Algorithm Programming Contest (BAPC 14) 部分题解

RT
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:深度学习:神经网络中的前向传播和反向传播算法推导
举报原因:
原因补充:

(最多只允许输入30个字)