就是推出循环矩阵乘积

然后一次操作后得到的c矩阵第一行第i列就是i的情况（b矩阵下标是a矩阵下标的转置）

两个循环矩阵乘积还是循环矩阵

以此推式子，发现c矩阵的第一行可以用a，b的第一行用循环卷积的形式表示

循环卷积也有结合律，可以快速幂

得到的多项式就是最终c矩阵第一行，直接输出即可。

emmm。。。

其实不用循环矩阵（这样只是说得方便严谨）

把b多项式reverse

画图可以得到

c=a*b *是循环卷积。

 

红色插入补充一点：

对于每个ai，质因数分解，枚举出约数。

开2e7个vector记录每个数是哪些数的约数，然后枚举每个数作为约数，遍历vector进行dfs

总复杂度是约数个数的5e7左右

但是2e7个vector会爆炸

开大数组然后分配内存：

cnt[2e7]表示是i的倍数的点多少个。i|d,cnt[i]++

通过cnt大小分配位置

再进行一次，把每个数填进去（这里可以开1e5个vector存约数，就不用再分解了）

 

第一步这个反演类似莫比乌斯反演

不同的是这是枚举倍数

分析法可以证明（把结果的g用条件换掉，然后大力反演，最后用f*[n==1]得到f）