sed 'a[digital_root(]Pdn[)=]P1-9%1+p'|dc