Sigma Function (因子和)

最新推荐文章于 2020-09-23 14:19:27 发布

Astronomical

最新推荐文章于 2020-09-23 14:19:27 发布

阅读量2.5k

点赞数 10

分类专栏：数论

本文链接：https://blog.csdn.net/strangedbly/article/details/50908522

版权

数论专栏收录该内容

18 篇文章 0 订阅

订阅专栏

Sigma Function

Description

Sigma function is an interesting function in Number Theory. It is denoted by the Greek letter Sigma (σ). This function actually denotes the sum of all divisors of a number. For example σ(24) = 1+2+3+4+6+8+12+24=60. Sigma of small numbers is easy to find but for large numbers it is very difficult to find in a straight forward way. But mathematicians have discovered a formula to find sigma. If the prime power decomposition of an integer is

Then we can write,

For some n the value of σ(n) is odd and for others it is even. Given a value n, you will have to find how many integers from 1 to n have even value ofσ.

Input

Input starts with an integer T (≤ 100), denoting the number of test cases.

Each case starts with a line containing an integer n (1 ≤ n ≤ 10¹²).

Output

For each case, print the case number and the result.

Sample Input

100

1000

Sample Output

Case 1: 1

Case 2: 5

Case 3: 83

Case 4: 947

题意：1~n (n:1~10¹²)中，因子和为偶数的有几个。
题解：My New Blog

因子和 Sum=(p1^0+p1^1….p1^e1)*(p2^0+p2^1…p2^e2)……(pn^0+…pn^en);

(p1^0+p1^1….p1^e1)，(p2^0+p2^1…p2^e2)，……(pn^0+…pn^en)中只要有一个是偶数，因子和sum就为偶数。所以只要找到一个是偶数就可以了。

若pi为偶数，则pi^x(x>0)为偶数，而pi^0=1（1+偶+偶….为奇数)。So,(pi^0+pi^1+…pi^ei)为奇数。所以pi只能是奇数,才能使(pi^0+pi^1+…pi^ei)为偶数。

再看pi^x (若x为奇数，pi^x为奇数(奇*奇*…为奇))，So,(pi^0+pi^1+…pi^ei)为偶数（1+奇+奇…)

所以，对m素因子分解，只要存在一个pi,ei都为奇数的pi^ei,就能使sum为偶数。

然而这么做必定TLE(-｡-;)

再想想，1~n中有多少个数的素因子分解中存在pi^ei(奇^奇）。

3^(2k-1) * (1,2,3,4..n/3) （k>1&&3^(2k-1)<=n)

5^(2k-1) * (1,2,3,4…n/5)

….

prime[i]^(2k-1) * (1,2,3,…n/prime[i]) 能包含所有的解，然而还有好多重复的解。看着有点像容斥定理，但多了个条件。（prime[i]^1,prime[i]^3…容斥有问题)。想了好久，感觉这题好难，不会了。。

不会咋办，换个思路呗!

不过之前还是忍不住暴力了一下果断TLE。

 
  #include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
const int N=1e6+5;

bool vis[N];
int prime[N],cnt;
void is_prime()
{
    cnt=0;
    memset(vis,0,sizeof(vis));
    for(int i=2;i<N;i++)
    {
        if(!vis[i])
        {
            prime[cnt++]=i;
            for(int j=i+i;j<N;j+=i)
                vis[j]=1;
        }
    }
}

bool is_even(long long n)
{
    for(int i=0;i<cnt&&prime[i]*prime[i]<=n;i++)
    {
        int count=0;
        if(n%prime[i]==0)
        {
            while(n%prime[i]==0)
            {
                n/=prime[i];
                count++;
            }
            if(prime[i]&1)
            {
                if(count&1)
                    return true;
            }
        }
    }
    if(n>1&&(n&1))
        return true;
    return false;
}

int main()
{
    int t;
    cin>>t;
    is_prime();
    for(int kase=1;kase<=t;kase++)
    {
        long long n;
        cin>>n;
        long long count=0;
        for(long long i=1;i<=n;i++)
        {
            if(is_even(i))
                cout<<i<<endl;
                count++;
        }
        printf("Case %d: %d\n",kase,count);
    }
} 
  
 
   
 

反过来想，什么时候因子和是奇数呢？
由前面的分析(只要存在一个pi,ei都为奇数的pi^ei,就能使sum为偶数)，素因子分解后，全为奇^偶，偶^偶，偶^奇，因子和就是奇数。

2是唯一一个偶素数。(特别的就得拿出来分开考虑。。)

((数^(偶/2)*(数^(偶/2))....)^2,这不是平方数吗！（数指的是奇数||偶数）（包含了数^偶）

所以这些数是平方数 || 2*平方数。（想不通的可以拿奇^偶，偶^偶，偶^奇组合，发现全被平方数 || 2*平方数包含了）

 
   
  #include<iostream>
#include<cmath>
#include<cstdio>
using namespace std;
typedef long long LL;

int main()
{
    int t;
    scanf("%d",&t);
    for(int kase=1;kase<=t;kase++)
    {
        LL n;
        scanf("%lld",&n);
        LL num1=(LL)sqrt((double)n);
        LL num2=(LL)sqrt((double)n/2.0);
        printf("Case %d: %lld\n",kase,n-num1-num2);
    }
    return 0;
} 
  
 
  
 
 
 
   
 