class="markdown_views prism-atom-one-light">
Description
给出
n,K" role="presentation" style="position: relative;">n,K,令
si=∑j≤ijK" role="presentation" style="position: relative;">si=∑j≤ijK,求
∑1≤i≤nsi[gcd(i,n)=1]" role="presentation" style="position: relative;">∑1≤i≤nsi[gcd(i,n)=1]
Input
第一行一整数
T" role="presentation" style="position: relative;">T表示用例组数,每组用例首先输入一整数
K" role="presentation" style="position: relative;">K,之后输入一整数
m" role="presentation" style="position: relative;">m表示
n" role="presentation" style="position: relative;">n的素因子数量,最后输入
m" role="presentation" style="position: relative;">m对整数
pi,ai" role="presentation" style="position: relative;">pi,ai表示
n=∏i=1npiai" role="presentation" style="position: relative;">n=∏i=1npaii
(K≤105,∑K≤106,m≤20,pi,ai≤109,pi<pi+1)" role="presentation" style="position: relative;">(K≤105,∑K≤106,m≤20,pi,ai≤109,pi<pi+1)
Sample Input
2
1
2
2 1
3 1
233
1
23333 1
Sample Output
16
32727388
Solution
对所求表达式反演有
ans=∑i=1n[(i,n)=1]∑j=1ijK=∑i=1n∑d|(i,n)μ(d)S(i,K)=∑d|nμ(d)∑i=1ndS(id,K)" role="presentation">ans===∑i=1n[(i,n)=1]∑j=1ijK∑i=1n∑d|(i,n)μ(d)S(i,K)∑d|nμ(d)∑i=1ndS(id,K)
其中
S(n,k)=∑i=1nik=1k+1∑i=1k+1Ck+1iBk+1−i(n+1)i" role="presentation" style="position: relative;">S(n,k)=∑i=1nik=1k+1∑i=1k+1Ci