HDU-5528-Count a * b-2015长春B题(数学推导)

最新推荐文章于 2021-01-24 22:09:46 发布

lyc1635566ty

最新推荐文章于 2021-01-24 22:09:46 发布

阅读量579

点赞数

分类专栏：数学推导数论文章标签：数学推导

本文链接：https://blog.csdn.net/lyc1635566ty/article/details/52781513

版权

数学推导同时被 2 个专栏收录

25 篇文章 0 订阅

订阅专栏

数论

12 篇文章 0 订阅

订阅专栏

题目链接

题意:给定一个n，求 $g(n)$ $g(n)=f(x_1)+f(x_2)...+f(x_m)$ ,其中 $x_i$ 是 $n$ 的约数， $f(n)是取(a，b)，a<n,b<n$ ，n不能整除 $a*b$ 的对数。

思路：进行公式推导把。

$f(x)=\sum_{a=0}^{x-1} \sum_{b=0}^{x-1}[x\nmid ab]\\=x^2-\sum_{a=0}^{x-1} \sum_{b=0}^{x-1}[x\mid ab]\\=x^2-\sum_{a=0}^{x-1} \sum_{b=0}^{x-1}[\frac{x}{gcd(x,a)}\mid \frac{a}{gcd(x,a)}*b]$

这里b只能取x的倍数，所以公式可以化简为

$=x^2-\sum_{a=0}^{x-1}\frac{x}{\frac{x}{gcd(x,a)}}\\=x^2-\sum_{a=0}^{x-1}{gcd(x,a)}=\\=x^2-\sum_{d|x}{d*euler(x/d)}$

$g(n)=\sum_{x|n}{f(x)}\\=\sum_{x|n}{x^2-\sum_{d|x}{deuler(x/d)}}\\=\sum_{x|n}{x^2}-\sum_{x|n}\sum_{d|x}{deuler(\frac{x}{d})}\\转化一下\\=\sum_{x|n}{x^2}-\sum_{d|n}\sum_{i|\frac{n}{d}}{deuler(i)}\\=\sum_{x|n}{x^2}-\sum_{d|n}d\sum_{i|\frac{n}{d}}{}euler(i)\\根据一个性质：一个数的约数的欧拉函数和还为这个数原式可化简为\\=\sum_{x|n}{x^2}-\sum_{\frac{d}{n}}{d\frac{n}{d}}\\=\sum_{x|n}{x^2}-\sum_{d|n}{n}\\=\sum_{x|n}{x^2}-\sum_{x|n}{x^2}-nh(n) （其中h(n)是n约数的个数）$

最后由于积性函数的性质

$g (n) = \prod i = 1 r ( p a i + 1 i ) 2 - 1 p 2 i - 1 - n * \prod i = 1 r a i + 1$ $g(n)=\prod_{i=1}^{r}{\frac{(p_{i}^{a_i+1})^2-1}{p_{i}^{2}-1}}-n*\prod_{i=1}^{r}{a_i+1}$

#include<bits/stdc++.h>
using namespace std;
#define maxn 100010
#define ull unsigned long long
bool vis[maxn];
int prime[maxn];
int tp=0;
void init(void) {
    memset(vis,0,sizeof vis);
    vis[1]=vis[0]=1;
    for(int i=1;i<maxn;i++) {
        if(!vis[i]) {
            prime[++tp]=i;
            for(int j=i+i;j<maxn;j+=i)
                vis[j]=1;
        }
    }
}
void sov(int n) {
    ull ans=1;
    ull ss=n;
    for(int i=1;prime[i]*prime[i]<=n;i++) {
        if(n%prime[i]==0) {
      //      cout << "n=" << n <<endl;
    //        cout << "prime=" << prime[i] << endl;
            int cnt=0;
            ull mul=1;
            while(n%prime[i]==0) {
                n/=prime[i];
                cnt++;
                mul*=prime[i];
            }
    //        cout << "mul" << mul << endl;
            mul*=prime[i];
            ss*=(cnt+1);
          //  ans*=((mul*mul)-1)/(prime[i]*prime[i]-1);
            ull a=(mul-1)/(prime[i]-1);
            ull b=mul+1;
            ull c=prime[i]+1;
            ans*=((a/c)*(b/c)*c)+a%c*(b/c)+b%c*(a/c);
        }
    }
  //  cout << ans <<endl;
    if(n>1) {
        ss*=2;
        ans*=(1+(ull)n*(ull)n);
    }
  //  cout << ans << " " << ss << endl;
    printf("%llu\n",ans-ss);
}
int main()
{
    init();
    int t,n;
    scanf("%d",&t);
    while(t--) {
        scanf("%d",&n);
        sov(n);
    }
    return 0;
}