题目链接:
HDU 6053
题意:
给出你一个长度为
n
的
题解:
容斥原理+莫比乌斯函数。
答案就是
∑ min(A) i=2 −mu(i)∗(∏ n i=1 ⌊A i k ⌋)
。
对于
f(k)
那么任意
B
都应该含有
那么
B i
可选择的方案数就是
⌊A i k ⌋
.
根据乘法原理:
f(k)=∏ n i=1 ⌊A i k ⌋
。
因为
gcd(1,n)>=2
,那么我们可以考虑分别计算
gcd
是
2
的倍数,
设
cnt[a]
是
A
中
因为
A
中
cnt[(t+1)k−1]−cnt[tk−1]
。
那么,那么
A
中与
t cnt[(t+1)k−1]−cnt[tk−1]
。
总复杂度
O(nlognlogn)
。
AC代码:
#include<bits/stdc++.h>
using namespace std;
const int mod=1e9+7;
typedef long long ll;
int mu[100005];
void mobius(int n)
{
mu[1]=1;
for(int i=1;i<=n;i++){
for(int j=i+i;j<=n;j+=i){
mu[j] -= mu[i];
}
}
}
int cnt[200005];
ll ans;
ll q_mod(ll x,int y)
{
ll res=1;
while(y)
{
if(y&1){
res = res * x % mod;
}
y >>= 1;
x = x * x % mod;
}
return res;
}
int main()
{
int cas = 0;
int n,t,Min;
mobius(100000);
scanf("%d",&t);
while(t--)
{
scanf("%d",&n);
Min=100000;
for(int i=0;i<=200000;i++) cnt[i]=0;
while(n--)
{
int x;
scanf("%d",&x);
cnt[x]++;
Min=min(Min,x);
}
ll ans=0;
for(int i=1;i<=200000;i++)
{
cnt[i]+=cnt[i-1];
}
for(int i=2;i<=Min;i++)
{
ll temp=1;
for(int j=1;i*j<=100000;j++)
{
temp=(temp * q_mod(1LL*j,cnt[i*j+i-1]-cnt[i*j-1]))%mod;
}
ans=(ans-temp*mu[i]+mod)%mod;
}
printf("Case #%d: ",++cas);
printf("%lld\n",ans);
}
return 0;
}