from fractions import Fraction as F
from itertools import product
def run(load,check):
 reg=load('C286_DATA.json')['target_register'];bios=load('C322_DATA.json')['biographies'];rows=[];phases=[];summaries=[];endpoints=[]
 for c,t in product([0,1],repeat=2):
  context=f'{c}{t}';state=[]
  for bio in bios:
   g=bio['gear'];mt={}
   for role in bio['roles']:
    matches=[x for x in reg if x['tradition']=='MT' and x['gear']==g and x['role']==role and x['Cainan']==c and x['Terah60']==t and x['full430']==1]
    check(f'{context}/G{g}/{role}/unique',len(matches)==1);mt[role]=matches[0]
    delta=bio['SP_records'][role]['center']-mt[role]['center'];check(f'{context}/G{g}/{role}/delta',delta==650-130*c-60*t)
    endpoints.append({'context':context,'gear':g,'role':role,'SP':bio['SP_records'][role],'MT':mt[role],'delta':delta})
   for a,terminal in product(bio['roles'],repeat=2):
    A=bio['SP_records'][a]['center'];B=mt[a]['center'];C=bio['SP_records'][terminal]['center'];d=A-C;r=B-C;delta=A-B
    row={'context':context,'gear':g,'anchor_role':a,'terminal_role':terminal,'A_SP':A,'B_MT':B,'C_SP':C,'delta':delta,'source_duration':d,'residual':r,'sign':'positive' if r>0 else 'zero' if r==0 else 'negative'};rows.append(row);state.append(row)
    for p,q,s in product([F(-1,4),F(1,4)],repeat=3):
     vals=[F(A)+p,F(B)+q,F(C)+s];exact=(vals[0]-vals[1],vals[0]-vals[2],vals[1]-vals[2])==(delta,d,r)
     phases.append({'context':context,'gear':g,'anchor_role':a,'terminal_role':terminal,'phases':list(map(str,[p,q,s])),'coordinates':list(map(str,vals)),'exact':exact,'positive_residual':r>0})
   check(f'{context}/G{g}/same biography',[mt[x]['center']-mt[y]['center'] for x,y in zip(bio['roles'],bio['roles'][1:])]==[500,100,350,150])
  summaries.append({'context':context,'Cainan':c,'Terah60':t,'delta':650-130*c-60*t,'annual_positive':sum(x['residual']>0 for x in state),'annual_zero':sum(x['residual']==0 for x in state),'annual_negative':sum(x['residual']<0 for x in state),'positive_residual_values':sorted({x['residual'] for x in state if x['residual']>0}),'Noah_lifespan_residual':next(x['residual'] for x in state if x['anchor_role']=='Noah birth' and x['terminal_role']=='Noah death')})
 check('allphasecriterion',all(x['exact']==(len(set(x['phases']))==1) for x in phases))
 for x in load('PREVIOUS_DATA.json')['rows']:
  y=next(y for y in rows if y['context']=='11' and (y['gear'],y['anchor_role'],y['terminal_role'])==(x['gear'],x['anchor_role'],x['terminal_role']))
  check(f'11 replicates {x["gear"]}/{x["anchor_role"]}/{x["terminal_role"]}',all(y[k]==x[k] for k in ['source_duration','residual','sign','A_SP','B_MT','C_SP']))
 counts={'source_endpoint_pairs':len(endpoints),'annual_candidates':len(rows),'phase_candidates':len(phases),'exact_partitions':sum(x['exact'] for x in phases),'positive_annual':sum(x['residual']>0 for x in rows),'positive_exact':sum(x['exact'] and x['positive_residual'] for x in phases)}
 check('counts',counts=={'source_endpoint_pairs':60,'annual_candidates':300,'phase_candidates':2400,'exact_partitions':600,'positive_annual':48,'positive_exact':96})
 return {'step':'C328','endpoint_pairs':endpoints,'rows':rows,'phase_trials':phases,'context_summaries':summaries,'counts':counts,'formula':'delta=650−130c−60t; r=d−delta; SP OFF/T0 and bothfull430 held.','source_selector_established':False,'new_contexts_adopted':False}