結果
問題 | No.1905 PURE PHRASE |
ユーザー | U S |
提出日時 | 2022-04-15 22:47:56 |
言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
結果 |
AC
|
実行時間 | 636 ms / 2,000 ms |
コード長 | 1,584 bytes |
コンパイル時間 | 123 ms |
コンパイル使用メモリ | 12,672 KB |
実行使用メモリ | 46,620 KB |
最終ジャッジ日時 | 2024-12-25 01:56:48 |
合計ジャッジ時間 | 28,363 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 38 |
ソースコード
# import sys # input = sys.stdin.readline # import pypyjit # pypyjit.set_param('max_unroll_recursion=-1') def mp():return map(int,input().split()) def lmp():return list(map(int,input().split())) def mps(A):return [tuple(map(int, input().split())) for _ in range(A)] def stoi(LIST):return list(map(int,LIST)) def itos(LIST):return list(map(str,LIST)) def bitA(X,A):return X & 1<<A == 1<<A import math import bisect import heapq import time from copy import copy as cc from copy import deepcopy as dc from itertools import accumulate, product from collections import Counter, defaultdict, deque def ceil(U,V):return (U+V-1)//V def modf1(N,MOD):return (N-1)%MOD+1 inf = int(1e18+20) mod = 998244353 import numpy as np n = int(input()) a = lmp() hz = [261.6,294.3,327,348.8,392.4,436,490.5] al = ["C","D","E","F","G","A","B"] F = np.fft.fft(a) amp = np.abs(F) ans = -1 h = -inf for i in range(510): if amp[i] > h: ans = i h = amp[i] er = inf now = -1 for i in range(7): ner = abs(ans-hz[i]) if ner < er: er = ner now = i print("{}4".format(al[now])) # def calc_var(lst): # mean = sum(lst)/len(lst) # var = 0 # for i in lst: # var += (i-mean)**2 # var /= len(lst) # return var # # # ans = -1 # er = inf # for i in range(7): # now = [[] for _ in range(hz[i])] # for j in range(n): # now[j%hz[i]].append(a[j]) # v = [] # for j in now: # v.append(calc_var(j)) # ter = max(v) # if ter < er: # er = ter # ans = i # print(ter) # print("{}4".format(al[ans])) #