Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- In[59]=
- ClearAll[a];
- a[n_]:=Plus[1/n c[n],Plus @@ (Function[m,(-m/n) a[m]^(n/m)] /@ Drop[Divisors[n],-1])] ;
- a[30]
- 2 3 5
- -c[1] c[2] 15 -c[1] c[3] 10 -c[1] c[5] 6
- 30 (------ + ----) (------ + ----) (------ + ----)
- -c[1] 2 2 3 3 5 5
- Out[59]= ------- - ----------------- - ----------------- - ---------------- -
- 30 15 10 6
- 2 3
- -c[1] c[2] 3 -c[1] c[3] 2
- 6 (------ + ----) (------ + ----)
- -c[1] 2 2 3 3 c[6] 5
- (------ - ---------------- - ---------------- + ----)
- 6 3 2 6
- > ------------------------------------------------------ -
- 5
- 2 5
- -c[1] c[2] 5 -c[1] c[5] 2
- 10 (------ + ----) (------ + ----)
- -c[1] 2 2 5 5 c[10] 3
- (------- - ---------------- - ---------------- + -----)
- 10 5 2 10
- > -------------------------------------------------------- -
- 3
- 3 5
- -c[1] c[3] 5 -c[1] c[5] 3
- 15 (------ + ----) (------ + ----)
- -c[1] 3 3 5 5 c[15] 2
- (------- - ---------------- - ---------------- + -----)
- 15 5 3 15 c[30]
- > -------------------------------------------------------- + -----
- 2 30
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement