print('Hello world!')
import datetime as dt
t0=733500
Dm=lambda y,m,d:(dt.date(y,m,d)-dt.date(2009,4,1)).days
fy=[(3496,4182),(3640,4355),(3735,4469),(3780,4523),(3710,4438),(3819,4569),(3867,4627),(4044,4838),(4090,4893),(4074,4874),(3945,4720)]
P=sum([[(t0+Dm(2009+i,4,1),n),(t0+Dm(2009+i,12,16),n),(t0+Dm(2010+i,1,15),s),(t0+Dm(2010+i,2,13),s),(t0+Dm(2010+i,3,15),n),(t0+Dm(2010+i,3,31),n)] for i,(n,s) in enumerate(fy)],[])
pub=[23046,23263,23286,23378,23686,23891,24161,24392,24621,25155,25661]
B=sum([[(t0+Dm(2009+i,4,1),v),(t0+Dm(2010+i,3,31),v)] for i,v in enumerate(pub)],[])
pri=[3946,4098,4033,3882,3906,4014,4226,4644,4657,5056,5050]
R=[(t0,4022),(t0+Dm(2009,12,31),4022)]+sum([[(t0+Dm(2010+i,1,1),v),(t0+Dm(2010+i,12,31),v)] for i,v in enumerate(pri)],[])
al=[7.5,7.5,7.2,7.5,7.4,7.3,7.2,7.1,7.2,7.2,7.5]
AA=sum([[(t0+Dm(2009+i,4,1),v),(t0+Dm(2010+i,3,31),v)] for i,v in enumerate(al)],[])
T=[(t0,0),(t0+Dm(2019,3,31),0),(t0+Dm(2019,4,1),8000)]
AC=[(t0+Dm(2008,3,31),0.2353),(t0+Dm(2009,12,31),0.2511),(t0+Dm(2011,12,1),0.2722),(t0+Dm(2014,6,1),0.3243),(t0+Dm(2016,12,1),0.3439),(t0+Dm(2018,12,31),0.3531),(t0+Dm(2019,3,31),0.3542),(t0+Dm(2019,12,31),0.3677)]
IV=lambda p,x:(p[0][1] if x<=p[0][0] else (p[-1][1] if x>=p[-1][0] else sum([aa[1]+(bb[1]-aa[1])*(x-aa[0])/(bb[0]-aa[0]) for aa,bb in zip(p,p[1:]) if aa[0]<=x<bb[0] and bb[0]>aa[0]])))
N=3928
DN=[IV(P,t0+k) for k in range(N)]
PB=[IV(B,t0+k) for k in range(N)]
RB=[IV(R,t0+k) for k in range(N)]
AL=[IV(AA,t0+k) for k in range(N)]
TX=[IV(T,t0+k) for k in range(N)]
AH=[IV(AC,t0+k) for k in range(N)]
print(DN[0],DN[259],DN[260],PB[0],PB[365],RB[274],RB[275],AL[0],TX[3652],AH[275])
BASE=0.222392
OD=0.5
res=365.0/3
UH=0.57
SP=0.06
SF=0.39
BD=0.055
TD=0.138
a1=-0.0275
a2=a1
a3=a1
b1=0.0
b2=0.0
b3=0.0
dgi=0.0
ipu=0.0
ipr=0.0
H=[]
for k in range(N): HI=BASE+a3+b3; H.append(round(HI*100,4)); qq=DN[k]; gsp=qq*HI*UH+qq*(1-HI)*SP; gsu=qq*HI*(1-UH)+qq*(1-HI)*(1-SP); la=AL[k]; dif=ipu/PB[k]-ipr/RB[k]; xx=(dif-BD)*OD; yy=(TX[k]*0.85*0.06/2876.0)*TD; c1=a1+(xx-a1)/res; c2=a2+(a1-a2)/res; c3=a3+(a2-a3)/res; e1=b1+(yy-b1)/res; e2=b2+(b1-b2)/res; e3=b3+(b2-b3)/res; dgi=dgi+(qq-gsp-gsu); ipu=ipu+(gsu-ipu/la); ipr=ipr+(gsp-ipr/(la*SF)); a1=c1; a2=c2; a3=c3; b1=e1; b2=e2; b3=e3
print(H[0],H[1],H[5],H[30],H[60],H[120],H[200],H[274],H[275],H[276])
print(H[275],H[1887],H[3652],H[3926])
print(AH[275],AH[1887],AH[3652],AH[3926])
print(H[1887]-32.96,H[3652]-35.08,H[3926]-36.61)