这个题最显眼的地方是10^9的前缀和和已知因子和的结果了,显然需要杜教筛。。
按照杜教筛的套路,我们设结果为F(n),那么
如果直接算空间和时间都承受不住。。需要预处理。。然后又发现一个大麻烦。。
首先可以发现无论是这个函数还是其和函数都不满足积性性质(f(1)!=1),所以需要转换思路。。。
运用莫比乌斯反演可以得到
然后发现f(n)就是三个积性函数的和。。而其中第二项其实就是欧拉函数 (推了半天才发现),第三项其实就是元函数e。。
然后就只剩下第一项了,看来只能用普通积性函数的分析方法分析一波。。令其为g(n),那么f(n)=g(n)-3φ(n)+2e(n)
而
然后就可以线性筛g(n)了,和φ(n)挺类似的其实。。
然后题目又提示了超过10^6的数据不超过5个,所以干脆直接预处理到10^6,也刚好能形成的复杂度。。
预处理完上杜教筛,解决
/**
* ┏┓ ┏┓
* ┏┛┗━━━━━━━┛┗━━━┓
* ┃ ┃
* ┃ ━ ┃
* ┃ > < ┃
* ┃ ┃
* ┃... ⌒ ... ┃
* ┃ ┃
* ┗━┓ ┏━┛
* ┃ ┃ Code is far away from bug with the animal protecting
* ┃ ┃ 神兽保佑,代码无bug
* ┃ ┃
* ┃ ┃
* ┃ ┃
* ┃ ┃
* ┃ ┗━━━┓
* ┃ ┣┓
* ┃ ┏┛
* ┗┓┓┏━━━━━━━━┳┓┏┛
* ┃┫┫ ┃┫┫
* ┗┻┛ ┗┻┛
*/
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<queue>
#include<cmath>
#include<map>
#include<stack>
#include<set>
#include<bitset>
#define inc(i,l,r) for(int i=l;i<=r;i++)
#define dec(i,l,r) for(int i=l;i>=r;i--)
#define link(x) for(edge *j=h[x];j;j=j->next)
#define mem(a) memset(a,0,sizeof(a))
#define ll long long
#define eps 1e-8
#define succ(x) (1LL<<(x))
#define lowbit(x) (x&(-x))
#define sqr(x) ((x)*(x))
#define mid ((x+y)/2)
#define NM 1000005
#define nm 10005
#define N 1000005
#define M(x,y) x=max(x,y)
const double pi=acos(-1);
const int inf=1e9+7;
using namespace std;
ll read(){
ll x=0,f=1;char ch=getchar();
while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
while(isdigit(ch))x=x*10+ch-'0',ch=getchar();
return f*x;
}
int m,prime[NM],tot,phi[NM];
ll _x,f[NM],inv[10],g[NM];
void init(){
bool v[NM];
mem(v);
m=1e6;inv[1]=1;phi[1]=g[1]=1;
inc(i,2,3)inv[i]=inv[inf%i]*(inf-inf/i)%inf;
inc(i,2,m){
if(!v[i])prime[++tot]=i,phi[i]=i-1,g[i]=((ll)i*i-1)%inf;
inc(j,1,tot){
if(i*prime[j]>m)break;
v[i*prime[j]]++;
if(i%prime[j]==0){
phi[i*prime[j]]=phi[i]*prime[j];
g[i*prime[j]]=g[i]*prime[j]%inf*prime[j]%inf;
break;
}
phi[i*prime[j]]=phi[i]*(prime[j]-1);
g[i*prime[j]]=g[i]*((ll)prime[j]*prime[j]-1)%inf;
}
f[i]=g[i]-3*phi[i];f[i]=(f[i-1]+f[i]+inf)%inf;
}
}
ll v[1005];
ll solve(ll n){
if(n<=m)return f[n];
int t=_x/n;
if(v[t])return v[t];
v[t]=(ll)n*(n-1)%inf*(n-2)%inf*inv[3]%inf;
for(ll i=2,j;i<=n;i=j+1){
j=n/(n/i);
v[t]=(v[t]-(j-i+1)*solve(n/i)%inf+inf)%inf;
}
return v[t];
}
int main(){
init();
int _=read();while(_--){
_x=read();if(_x>m)mem(v);
printf("%lld\n",solve(_x));
}
return 0;
}
functionTime Limit: 6000/3000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others)Total Submission(s): 954 Accepted Submission(s): 367 Problem Description There is a function f(x),which is defined on the natural numbers set N,satisfies the following eqaution
Input the first line contains a positive integer T,means the number of the test cases.
Output Tlines,each line contains a number,means the answer to the i-th test case.
Sample Input 1 3
Sample Output 2 $1^2-3*1+2=f(1)=0$ $2^2-3*2+2=f(2)+f(1)=0->f(2)=0$ $3^2-3*3+2=f(3)+f(1)=2->f(3)=2$ $f(1)+f(2)+f(3)=2$
Source
Recommend hujie | We have carefully selected several similar problems for you: 6447 6446 6445 6444 6443
| ||
|