Help Tomisu UVA - 11440（phi函数的定义 + 思维 + gcd的思想）

Help Tomisu UVA - 11440

题意：

给你一个n和m。统计 2–n！，有多少个x的所有素因子都大于m。

思路：

1. x的所有素因子都大于m。等价于 x 和 $m!$ 互素。
2. 假如 x > $m!$ 。根据gcd里辗转相除法的定义，可以知道 x 和 $m!$ 互素 当且仅当 x% $m!$ 与 $m!$ 互素。这样只需要求出不超过 $m!$ 且与 $m!$ 互素的正整数个数即可。最后ans再乘 $\frac{n!}{m!}$即可。（因为 $n!$ 是 $m!$ 的倍数）
3. 现在问题就是怎么求 phifac【m】。

$\varphi (n)=n(1-\frac{1}{p1})(1-\frac{1}{p2})...(1-\frac{1}{pn})$

根据上面欧拉函数的定义有

1. 假如n不是质数，那么phifac【n】= n*phifac【n-1】。($n!$ 和 $（n-1）!$的质因子相同)
2. 假如n是质数，那么phifac【n】= n*phifac【n-1】$(1-\frac{1}{n})$ = phifac【n-1】（n-1）

反思：

本题的模数是1e8.。。。。。wa了几发。。。。。

AC

#include <iostream>
#include <algorithm>
#include <cmath>
#define For(i,x,y) for(int i=(x); i<=(y); i++)
#define fori(i,x,y) for(int i=(x); i<(y); i++)
using namespace std;
typedef long long ll;
const int maxn = 1e7+10;
const int mod = 1e8+7;
int vis[maxn], prime[maxn];
bool judge(ll x){
    if(x == 1)return false;
    for(ll i = 2; i*i<=x; i++){
        if(x%i==0)return false;
    }
    return true;
}
ll phifac[maxn];
void table(){
    fori(i,0,maxn) vis[i] = 0;
    fori(i,2,maxn){
        if(!vis[i]){
            prime[i] = 1;
            for(int j = i; j < maxn; j += i) vis[j] = 1;
        }
    }
    phifac[1] = phifac[2] = 1;
    fori(i,3,maxn) phifac[i] = phifac[i-1]*(prime[i] ? (i-1) : (i)) % mod;
}
int main()
{
    ios::sync_with_stdio(0);cin.tie(0);cout.tie(0);
    table();
    //for(int i=(int)1e7; i<=maxn; i++)if(judge(i))cout<<i<<' '<<endl;
    int n, m;
    while(cin>>n>>m,n||m){
        int ans = phifac[m];
        For(i,m+1,n) ans = (ll)ans * i % mod;
        cout<<((ans - 1 + mod) % mod)<<endl;
    }
    return 0;
}