uva 11440 Help Mr. Tomisu 欧拉函数变种

Problem:
给定正整数N和M,统计2和N!之间有多少个整数满足所有素因子都大于M,输出答案除以100000007的余数。
Solution:
所有素因子都大于M意味着这个整数x和M!互素,互素?想到了欧拉函数,而x和M!互素也就是说x%(M!)与M!互素,那我们最终要求的也就是(N!/M!)*phi(M!),然后在运算过程中作相应的取余运算即可。
note:
1. 要能够提炼出来模型转换为和M!互素。
2. 互素要知道x和m互素就是x%m和m互素。
3. 互素要想到欧拉函数。

#include<cstdio>
#include<iostream>
#include<sstream>
#include<cstdlib>
#include<cmath>
#include<cctype>
#include<string>
#include<cstring>
#include<algorithm>
#include<stack>
#include<queue>
#include<set>
#include<map>
#include<ctime>
#include<vector>
#include<fstream>
#include<list>
using namespace std;

#define ms(s) memset(s,0,sizeof(s))
typedef unsigned long long ULL;
typedef long long LL;

const int INF = 0x3fffffff;

const int N = 10000000;
bool primeTable[N+5];

void make_primeTable(){
    fill(primeTable,primeTable+N,1);
    primeTable[0] = false;
    primeTable[1] = false;

    int maxed = sqrt(N);
    for(int i = 2; i <= maxed; ++i){
        if(primeTable[i]){
            for(int j = i*i; j <= N; j += i)
                primeTable[j] = false;
        }
    }
}
const int maxn = 10000000;
const int mod = 100000007;
long long phi[maxn+5];
void phi_table(){
    memset(phi,0,sizeof(phi));
    phi[1] = phi[2] = 1;
    for(int i = 3; i <= maxn; ++i){
        if(primeTable[i] == true)
            phi[i] = (i-1)*phi[i-1]%mod;
        else
            phi[i] = i*phi[i-1]%mod;
    }
}

int main(){
//    freopen("F:\\input.txt","r",stdin);
//    freopen("F:\\output.txt","w",stdout);
//    ios::sync_with_stdio(false);

    make_primeTable();
    phi_table();
    int n,m;
    LL ans;
    while((~scanf("%d%d",&n,&m)) && n){
        ans = phi[m];
        for(int i = n; i > m; --i){
            ans = ans*(LL)i%100000007LL;
        }
        printf("%lld\n",ans-1);
    }

    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值