题意:f(n)为n所有约数的和,给你一个数n,让你求从1到n中f(n)为偶数的数有多少个
题解:
唯一分解定理:x=p1^a1*p2*a2*…*pn*an
数x的因子和 f(x)= (1+p1+p1^2+p1^3+…+p1^a1)*(1+p2+p2^2+…+p2^a2)*…*(1+pn+pn^2+…+pn^an);//穷尽所有因子的多项式相加
(1)偶数乘偶数还是偶数 //所以偶数因子必定因为加1而括号内为奇数
(2)奇数乘奇数还是奇数 //所以某个奇数因子若有奇数个则括号内为偶数,偶数个则括号内为奇数
(3)奇数乘偶数是偶数 //所以因子和为奇数当且仅当所有括号内都为奇数
由此可得,若使得因子和为奇数,则除了2以外(只有2是偶素数),所有的素因子的次数必须为偶数,所以结果必定是某个数的平方(比如3^4*5^2==45^2),或者是某个数平方的2倍(其实可以是2的n次方倍,但若2的次数为偶数,本身就在平方之内,若为奇数,则提出一个2,剩余放入平方内),所以1~n中f(x)%2==0的数的个数为n-sqrt(n)-sqrt(n/2);

#include<algorithm>
#include <iostream>
#include  <cstdlib>
#include  <cstring>
#include  <cassert>
#include   <cstdio>
#include   <vector>
#include   <string>
#include    <cmath>
#include    <queue>
#include    <stack>
#include      <set>
#include      
<map>
using namespace std;
#define P(a,b,c) make_pair(a,make_pair(b,c))
#define rep(i,a,n) for (int i=a;i<=n;i++) #define per(i,a,n) for (int i=n;i>=a;i--)
#define CLR(vis) memset(vis,0,sizeof(vis))
#define MST(vis,pos) memset(vis,pos,sizeof(vis))
#define pb push_back
#define mp make_pair
#define all(x) (x).begin(),(x).end()
#define fi first
#define se second
#define SZ(x) ((int)(x).size())
typedef pair<int,pair<int,int> >pii;
typedef long long ll;
const int INF = 0x3f3f3f3f;
const ll mod = 1000000007;
const int maxn = 1e6+5;

int main(){
    int T;
    scanf("%d", &T);
    for(int kase = 1; kase <= T; kase++) { ll n, sum; cin>>n;
        sum = n;
        sum -= (int)sqrt(n);
        sum -= (int)sqrt(n/2);//可以适当提高精度
        printf("Case %d: ", kase);
        cout<<sum;
        if(kase!=T)cout<<endl;
    }
	return 0;
} 

发表评论

邮箱地址不会被公开。 必填项已用*标注