题意: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; }