""" 作戦 0: 等差数列で写真を取る 作戦 1: 上位 L 枚を採用 作戦 2: 山登り 作戦 3: アニーリング 作戦 4: 初期解を頑張る. 作戦 5: 作戦 4 の下でアニーリング 作戦 6: ビームサーチ """ #================================================== # 関数部 def debug(*message): print(*message,file=sys.stderr) def General_Binary_Increase_Search_Integer(L, R, cond, default=None): """ 条件式が単調増加であるとき, 整数上で二部探索を行う. L: 解の下限 R: 解の上限 cond: 条件(1変数関数, 広義単調増加を満たす) default: Lで条件を満たさないときの返り値 """ if not(cond(R)): return default if cond(L): return L R+=1 while R-L>1: C=L+(R-L)//2 if cond(C): R=C else: L=C return R def move(start,direction,time): return start+direction*time def position(i,T): if T0: T-=(B[i]-H[i]*E[i])+M[i] q,r=divmod(T,2*M[i]) direction=pow(-1,H[i]+2+q) start=pow(-1,H[i]+1+q)*M[i] return move(start,direction,r) else: return 0 def PROBABILITY(delta,temp): if delta>=0: return 1 else: return exp(delta/temp) Memo={} def picture(t): if t in Memo: return Memo[t] X=[position(i,t) for i in range(N)] score=0 for i in range(N): for j in range(i+1,N): score+=abs(X[i]-X[j])/(B[i]+B[j]) score*=pow(10,9)*(2/(N*(N-1))) Memo[t]=round(score) return Memo[t] #================================================== # インポート from operator import itemgetter from time import * from heapq import * from random import * from math import * import sys #================================================== # 入力 N,L,Delta=map(int,input().split()) B=[0]*N; M=[0]*N; E=[0]*N for i in range(N): B[i],M[i],E[i]=map(int,input().split()) #================================================== H=[0]*N; F=[0]*N for i in range(N): H[i]=(B[i]-M[i])//E[i] F[i]=2*B[i]*H[i]-E[i]*H[i]*H[i] #================================================== BEGIN_TIME=time() TIME_LIMIT=1.800 TEMPERATURE=5 ALPHA=0.999 BEAM_SIZE=20 TIME_DELTA=5 #================================================== A=[] X=[(picture(t),t) for t in range(Delta, Delta+BEAM_SIZE)] while time()-BEGIN_TIME=Delta if flag: T.append(u) if len(T)==L: break #================================================== T.sort() print(len(T)) print(*T)