弹药科技

题目描述
经过精灵族全力抵挡，精灵终于坚持到了联络系统的重建，于是精灵向人类求助，
大魔法师伊扎洛决定弓}用博士的最新科技来抗敌。
伊扎洛:“博士，还没好吗?”
博士:“只差一步了!只需要在正确的位置装上弹药就可以了!”
博士的最新科技是全新的炸弹，但是现在还需要一步装弹药的操作。博士的炸弹有
N!个位置可以装弹药(>.<)，但是只有在正确的位置装上弹药才能启动，博士将
装弹药的位置编号为1到N!，一个位置i需要装弹药，当且仅当gcd(i, N!) ≠ 1且
N! mod i ≠ 0，现在博士不要求你求出所有装弹药位置，只需要你求出满足要求
的位置的数量就可以了。

输入
一个数N

输出
表示满足要求的位置数量，答案对mod1000000007输出

样例输入
4
样例输出
9
提示
 

【数据范围】

N <= 1000000

 

这就是我说的板子题，再加上一点小小的容斥思想，具体公式为：

其中suma表示n!以内的与n互质的数，sumb表示n!的约数个数，由于1被减了两次，要加上一次

其实就是：既然求不满足，就把总的减去满足的就是解

suma就是欧拉函数，sumb可以直接用公式求

 

至于怎么求p数组与w数组，用勒让德定理即可（有道以前的有道题用过）

推荐求n的欧拉函数用第二个，不用涉及到除法

记得取模

代码:

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <cstdlib>
using namespace std;
#define ll long long
const int mod = 1e9 + 7;
const int MAXN = 1000003;
int n;
ll prime[MAXN] , pn , phi[MAXN];
bool vis[MAXN];
ll psum[MAXN];
void pr(){
    for( int i =2 ; i <= n ; i ++ ){
        if( vis[i] == 0 ){
            phi[i] = i - 1;
            prime[++pn] = i;
        }
        for( int j =1 ; j <= pn && 1ll * i * prime[j] <= n ; j ++ ){
            vis[i*prime[j]] = 1;
            if( i % prime[j] == 0 ){
                phi[i*prime[j]] = phi[i] * prime[j];
                break;
            }
            phi[i*prime[j]] = phi[i] * ( prime[j] - 1 );
        }
    }
}
ll qpow( ll x , ll y , ll mod ){
    ll sum = 1;
    while( y ){
        if( y % 2 ) sum = sum * x % mod;
        x = x * x % mod;
        y /= 2;
    }
    return sum;
}
int main(){
    scanf( "%d" , &n );
    pr();
    for( int i = 1 ; i <= pn ; i ++ ){
        ll k = 1;
        for( int j = 1 ; k * prime[i] <= n ; j ++ ){
            k = k * prime[i];
            psum[i] += n / k;
        }
    }
    ll sum = 1;
    for( int i = 2 ; i <= n ; i ++ ){
        sum = sum * i % mod;
    }
    ll ans = 1 , ans2 = 1;
    for( int i = 1 ; i <= pn ; i ++ ){
        ans = ans * ( qpow( prime[i] , psum[i] - 1 , mod ) * ( prime[i] - 1 ) % mod ) % mod;
        ans2 = ans2 * ( psum[i] + 1 ) % mod;
    }
    sum = sum - ans - ans2 + 1;
    while( sum < 0 ){
        sum += mod;
    }
    printf( "%lld" , sum );
    return 0;
}
 

--------------------- 
作者：BIT_jzx 
原文：https://blog.csdn.net/weixin_43823476/article/details/89077146