Pollux (./70) :
// mettre la somme de x^0 à x^(n-1) dans p (a priori l'expand() ne fait rien puisque n est pas connu, mais bon)
expand((x^n-1)/(x-1)) -> p

mod(subs(n,42,p),x-1,x) = 42
subs(n,42,mod(p,x-1,x)) = polynôme en NAN
subs(un,1,subs(n,42,mod(p,x-un,x))) = 42

Si déjà, ca marche comme çà direct, ca sera mieux que bien des CAS
(Et le résultat, c'est mod(42,x+1,x)

Pollux (./70) :
Comment est-ce que je peux me débrouiller pour que p soit vu comme un polynôme et pas comme une fraction rationnelle sans pour autant remplacer n par une constante ?

Je vois deux solutions:
+ créer une fonction 'id' identitée qui renvoie 'x' si 'x' est un entier/rationel/... ou polynome d'une variable exactement, ou id(x,flags) autrement (flags déterminant ce que l'on désire comme type). Une sorte de cast passif.
+ ajouter un parametre global qui fait que si on travaille mod 2, lorsqu'on tombe sur 1/2, on laisse 1/2, et on ne remplace pas par NAN..
+ même chose, mais au lieu d'être global, il est local et compris dans l'extension Mod.