結果

問題 No.195 フィボナッチ数列の理解(2)
ユーザー shotoyooshotoyoo
提出日時 2021-06-28 23:28:08
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
AC  
実行時間 204 ms / 5,000 ms
コード長 2,452 bytes
コンパイル時間 101 ms
コンパイル使用メモリ 11,180 KB
実行使用メモリ 42,816 KB
最終ジャッジ日時 2023-09-08 00:05:15
合計ジャッジ時間 7,671 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge15
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 204 ms
42,460 KB
testcase_01 AC 202 ms
42,448 KB
testcase_02 AC 204 ms
42,672 KB
testcase_03 AC 201 ms
42,432 KB
testcase_04 AC 203 ms
42,480 KB
testcase_05 AC 194 ms
42,664 KB
testcase_06 AC 195 ms
42,508 KB
testcase_07 AC 198 ms
42,452 KB
testcase_08 AC 188 ms
42,480 KB
testcase_09 AC 186 ms
42,796 KB
testcase_10 AC 198 ms
42,704 KB
testcase_11 AC 198 ms
42,404 KB
testcase_12 AC 199 ms
42,816 KB
testcase_13 AC 198 ms
42,788 KB
testcase_14 AC 190 ms
42,712 KB
testcase_15 AC 192 ms
42,496 KB
testcase_16 AC 199 ms
42,636 KB
testcase_17 AC 198 ms
42,644 KB
testcase_18 AC 193 ms
42,624 KB
testcase_19 AC 191 ms
42,424 KB
testcase_20 AC 196 ms
42,436 KB
testcase_21 AC 193 ms
42,444 KB
testcase_22 AC 194 ms
42,684 KB
testcase_23 AC 192 ms
42,376 KB
testcase_24 AC 196 ms
42,624 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys
input = lambda : sys.stdin.readline().rstrip()

sys.setrecursionlimit(2*10**5+10)
write = lambda x: sys.stdout.write(x+"\n")
debug = lambda x: sys.stderr.write(x+"\n")
writef = lambda x: print("{:.12f}".format(x))

def solve(a,b,c,d,e,f):
    """
    ax + by + c = 0
    dx + ey + f = 0
    """
    v = (a*e-b*d)
    if v==0:
        # 縮退している
        if b!=0:
            return (0, -c/b)
        else:
            return (-c/a,0)
    x = - (c*e - f*b) / v
    y = (d*c - a*f) / v
    return (x,y)
def solve(A,B):
    (a,b),(c,d) = A
    v = a*d-b*c
    assert v!=0
    return (d*B[0]-b*B[1])/v, (-c*B[0]+a*B[1])/v
import scipy.linalg
import numpy as np
x,y,z = list(map(int, input().split()))
f = [1,0]
for i in range(301):
    f.append(f[-2]+f[-1])
s = set([x,y,z])
ans = None
if len(s)==1:
    a = 1
    for i in range(50)[::-1]:
        if (x - f[i]) % f[i+1] == 0 and (x-f[i])//f[i+1]>0:
            b = (x-f[i])//f[i+1]
            ans = (a,b)
            break
    else:
        if x==1:
            ans = (1,1)
        else:
            ans = (1, x-1)
elif len(s)==2:
    x,y = s
    ans = None
    for i in range(50):
        for j in range(50):
            if i==j:
                continue
            v1,v2 = solve(([[f[i], f[i+1]], [f[j], f[j+1]]]), ([x,y]))
            v1,v2 = int(v1+0.1), int(v2+0.1)
            if not (f[i]*v1 + f[i+1]*v2 == x and f[j]*v1 + f[j+1]*v2==y):
                continue
            if v1>0 and v2>0 and (ans is None or (v1,v2)<ans):
                ans = (v1,v2)
    if ans is None:
        aa = min(x,y)
        bb = max(x,y)
        ans = min((aa,bb), (bb-aa, aa), (aa, bb-aa))
else:
    x,y,z = s
    ans = None
    for i in range(50):
        for j in range(50):
            if i==j:
                continue
            try:
                v1,v2 = solve(([[f[i], f[i+1]], [f[j], f[j+1]]]), ([x,y]))
            except Exception:
                continue
            v1 = int(v1+0.1)
            v2 = int(v2+0.1)
            if not (f[i]*v1 + f[i+1]*v2 == x and f[j]*v1 + f[j+1]*v2==y):
                continue
            if not (v1>0 and v2>0):
                continue
            for k in range(50):
                if f[k]*v1 + f[k+1]*v2 == z:
                    if v1>0 and v2>0:
                        if (ans is None or (v1,v2)<ans):
                            ans = (int(v1+0.1), int(v2+0.1))
                    
    if ans is None:
        ans = (-1,)
print(*ans)
0