escription
Lfx在复习离散的时候突然想到了一个算法题，毕竟是lfx，
算法题如下：
他想知道这样的问题，先定义1~n中即是3的倍数，又是11的倍数的那些数的和sum，
他想知道sum有多少个质因子，以及1~sum-1中有多少个数与sum互质？
1<= N <= 1e6
输入：
一个整数n
输出
两个整数，分别代表sum质因子的数量以及1~sum-1中与sum互质的数量。
思路：
先1~n扫一下求sum值，然后用唯一分解定理求质因子的数量，用欧拉函数求互质的数量。
唯一分解定理的步骤：
先打一个素数表，方法有很多种，然后用已知的素数去分解数值。
对于一个数x，小于x并与x互质的数的数量就是欧拉函数的定义，一个数论函数，很基础。
不知道的新名词应该去学习一下。
细节见代码：
#include <iostream> #include <cstdio> #include <cstring> #include <algorithm> #include <cmath> #include <queue> #include <stack> #include <map> #include <set> #include <vector> #include <iomanip> #define ALL(x) (x).begin(), (x).end() #define rt return #define dll(x) scanf("%I64d",&x) #define xll(x) printf("%I64d ",x) #define sz(a) int(a.size()) #define all(a) a.begin(), a.end() #define rep(i,x,n) for(int i=x;i<n;i++) #define repd(i,x,n) for(int i=x;i<=n;i++) #define pii pair<int,int> #define pll pair<long long ,long long> #define gbtb ios::sync_with_stdio(false),cin.tie(0),cout.tie(0) #define MS0(X) memset((X), 0, sizeof((X))) #define MSC0(X) memset((X), ' ', sizeof((X))) #define pb push_back #define mp make_pair #define fi first #define se second #define eps 1e-6 #define gg(x) getInt(&x) #define db(x) cout<<"== [ "<<x<<" ] =="<<endl; using namespace std; typedef long long ll; ll gcd(ll a,ll b){return b?gcd(b,a%b):a;} ll lcm(ll a,ll b){return a/gcd(a,b)*b;} ll powmod(ll a,ll b,ll MOD){ll ans=1;while(b){if(b%2)ans=ans*a%MOD;a=a*a%MOD;b/=2;}return ans;} inline void getInt(int* p); const int inf=0x3f3f3f3f; /*** TEMPLATE CODE * * STARTS HERE ***/ const int maxn = 1e6+50; bool noprime[maxn+50]; vector <int> p; int getPrime() { // 华丽的初始化 memset(noprime,false,sizeof(noprime)); p.clear(); int m=(int)sqrt(maxn+0.5); // 多姿的线性筛 for(int i=2;i<=m;i++) { if(!noprime[i]) { for(int j=i*i;j<=maxn;j+=i) { noprime[j] = true; } } } // 把素数加到vector里 for(int i=2;i<=maxn;i++) { if(!noprime[i]) { p.push_back(i); } } //返回vector的大小 return p.size(); } int pf[105][2];// 0 -> value 1->count int getPrifac( ll n,int len) { int pos = 0; for(int i=0; p[i]*p[i]<=n&&i<len;i++) { if( n% p[i] == 0) { pf[++pos][0]=p[i]; pf[pos][1]=0; // 算质因数的幂数 while(n%p[i]==0) { pf[pos][1]++; n/=p[i]; } } } if( n> 1) { pf[++pos][0] = n; pf[pos][1]=1; } return pos; // 优美的返回有多少个质因数 // 1~pos } ll euler(ll n) { //log(n)时间内求一个数的欧拉值 ll ans = n; for (ll i = 2; i*i <= n; i++) { if (n%i == 0) { ans -= ans / i; while (n%i == 0) n /= i; } } if (n>1) ans -= ans / n; return ans; } int main() { int len=getPrime(); int n; gg(n); ll cnt=0ll; repd(i,1,n) { if((i%3==0)&&(i%11==0)) { cnt+=i; } } // db(cnt); int num=getPrifac(cnt,len); printf("%d ",num); printf("%lld ",euler(cnt) ); return 0; } inline void getInt(int* p) { char ch; do { ch = getchar(); } while (ch == ' ' || ch == ' '); if (ch == '-') { *p = -(getchar() - '0'); while ((ch = getchar()) >= '0' && ch <= '9') { *p = *p * 10 - ch + '0'; } } else { *p = ch - '0'; while ((ch = getchar()) >= '0' && ch <= '9') { *p = *p * 10 + ch - '0'; } } }