結果
問題 |
No.1391 ±1 Abs Sum
|
ユーザー |
![]() |
提出日時 | 2021-02-12 23:07:52 |
言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
結果 |
TLE
|
実行時間 | - |
コード長 | 1,837 bytes |
コンパイル時間 | 121 ms |
コンパイル使用メモリ | 12,672 KB |
実行使用メモリ | 39,604 KB |
最終ジャッジ日時 | 2024-07-20 00:36:33 |
合計ジャッジ時間 | 4,789 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 8 TLE * 1 -- * 25 |
ソースコード
import sys from collections import defaultdict, Counter, deque from itertools import permutations, combinations, product, combinations_with_replacement, groupby, accumulate import operator from math import sqrt, gcd, factorial # from math import isqrt, prod,comb # python3.8用(notpypy) from bisect import bisect_left, bisect_right #from functools import lru_cache,reduce #from heapq import heappush,heappop,heapify,heappushpop,heapreplace #import numpy as np #import networkx as nx #from networkx.utils import UnionFind #from numba import njit, b1, i1, i4, i8, f8 #from scipy.sparse import csr_matrix #from scipy.sparse.csgraph import shortest_path, floyd_warshall, dijkstra, bellman_ford, johnson, NegativeCycleError # numba例 @njit(i1(i4[:], i8[:, :]),cache=True) 引数i4配列、i8 2次元配列,戻り値i1 def input(): return sys.stdin.readline().rstrip() def divceil(n, k): return 1+(n-1)//k # n/kの切り上げを返す def yn(hantei, yes='Yes', no='No'): print(yes if hantei else no) def main(): mod = 10**9+7 mod2 = 998244353 n, k = map(int, input().split()) A = list(map(int, input().split())) ans = 10**20 for i in range(n): ng, ok = -1, 2*10**9+10 # abs(ok)以下でk個揃えられるか? while (abs(ok - ng) > 1): mid = (ok + ng) // 2 r = bisect_right(A, A[i]+mid) l = bisect_left(A, A[i]-mid) if r-l >= k: ok = mid else: ng = mid r = bisect_right(A, A[i]+ok) l = bisect_left(A, A[i]-ok) plusone = sum(A[i:r])-A[i]*(r-i)-sum(A[l:i])+A[i]*(i-l)-ok*(k-(r-l)) minusone = sum(A[r:])-A[i]*(n-r)-sum(A[:l])+A[i]*l+ok*(k-(r-l)) ans = min(ans, plusone-minusone) #print(ans,l,r,ok,mid) print(ans) if __name__ == '__main__': main()