[HNOI2002][poj1091]跳蚤

原创 2017年01月03日 17:00:54

题目描述

  Z城市居住着很多只跳蚤。在Z城市周六生活频道有一个娱乐节目。一只跳蚤将被请上一个高空钢丝的正中央。钢丝很长,可以看作是无限长。节目主持人会给该跳蚤发一张卡片。卡片上写有N+1个自然数。其中最后一个是M,而前N个数都不超过M,卡片上允许有相同的数字。跳蚤每次可以从卡片上任意选择一个自然数S,然后向左,或向右跳S个单位长度。而他最终的任务是跳到距离他左边一个单位长度的地方,并捡起位于那里的礼物。
  比如当N=2,M=18时,持有卡片(10, 15, 18)的跳蚤,就可以完成任务:他可以先向左跳10个单位长度,然后再连向左跳3次,每次15个单位长度,最后再向右连跳3次,每次18个单位长度。而持有卡片(12, 15, 18)的跳蚤,则怎么也不可能跳到距他左边一个单位长度的地方。
  当确定N和M后,显然一共有MN张不同的卡片。现在的问题是,在这所有的卡片中,有多少张可以完成任务。


输入格式

输入文件有且仅有一行,包括用空格分开的两个整数N和M。


输出格式

输出文件有且仅有一行,即可以完成任务的卡片数。
1≤M≤10^8,1≤N≤M,且MN≤10^16。


样例数据

样例输入

2 4

样例输出

12


样例说明

这12张卡片分别是:
(1, 1, 4), (1, 2, 4), (1, 3, 4), (1, 4, 4), (2, 1, 4), (2, 3, 4),
(3, 1, 4), (3, 2, 4), (3, 3, 4), (3, 4, 4), (4, 1, 4), (4, 3, 4)


题目分析

这题稍微有点难。。
题意很好理解,即为不定方程

x1a1+x2a2+...+xn+1an+1=1

要有解,须满足系数最大公约数为1,详细见这儿
于是只需要将小于M中不与M互质的数去掉即可,利用容斥原理可以得到设t(k)为数列gcd为k个质因子乘积的数列个数
故Ans=m^n-t(1)+t(2)-t(3)…(-1)^k*t(k)


源代码

#include<algorithm>
#include<iostream>
#include<iomanip>
#include<cstring>
#include<cstdlib>
#include<vector>
#include<cstdio>
#include<cmath>
#include<queue>
using namespace std;
typedef long long LL;
inline const LL Get_Int() {
    LL num=0,bj=1;
    char x=getchar();
    while(x<'0'||x>'9') {
        if(x=='-')bj=-1;
        x=getchar();
    }
    while(x>='0'&&x<='9') {
        num=num*10+x-'0';
        x=getchar();
    }
    return num*bj;
}
LL n,m,Ans=0,tot,cnt=0,num[105],ans[105];
LL Quick_Pow(LL a,LL b) {
    LL ans=1;
    while(b>0) {
        if(b&1)ans*=a;
        b>>=1;
        a*=a;
    }
    return ans;
}
void Divide(LL x) {
    cnt=0;
    for(int i=2; i<=sqrt(x); i++)
        if(x%i==0) {
            num[++cnt]=i;
            while(x%i==0)x/=i;
        }
    if(x!=1)num[++cnt]=x;
}
void Dfs(LL Last,LL Limit,LL Now) {
    if(Now>Limit) {
        LL sum=m;
        for(int i=1; i<=Limit; i++)sum/=ans[i];
        tot+=Quick_Pow(sum,n);
        return;
    }
    for(int i=Last+1; i<=cnt; i++) {
        ans[Now]=num[i];
        Dfs(i,Limit,Now+1);
    }
}
int main() { 
    n=Get_Int();
    m=Get_Int();
    Divide(m);
    Ans=Quick_Pow(m,n);
    for(int i=1; i<=cnt; i++) {
        tot=0;
        memset(ans,0,sizeof(ans));
        Dfs(0,i,1);
        if(i&1)Ans-=tot;
        else Ans+=tot;
    }
    printf("%lld\n",Ans);
    return 0;
}

版权声明:本文为博主原创文章,未经博主允许不得转载。 举报

相关文章推荐

poj 1091 跳蚤(组合数学-容斥原理)这题要求的数据范围根本就是胡说八道!

跳蚤 Time Limit: 1000MS   Memory Limit: 10000K Total Submissions: 8129   Accepted: 2...

poj 1091 跳蚤(最大公约数原理+容斥原理)

题目:http://poj.org/problem?id=1091 跳蚤 Time Limit: 1000MS   Memory Limit: 10000K Total Submiss...

我是如何成为一名python大咖的?

人生苦短,都说必须python,那么我分享下我是如何从小白成为Python资深开发者的吧。2014年我大学刚毕业..

POJ 1091 跳蚤 解题报告(2014携程第一场A题,聪明的猴子)(容斥原理)

/*  * 大数模版  * 完美支持负数及加减法操作  * 支持乘法,大小比较,赋值  * 支持与long long类型的赋值比较  * SF-_- 14.04.01  */ #incl...
  • kbdwo
  • kbdwo
  • 2014-04-11 00:04
  • 976

POJ - 1091 跳蚤

积性函数

容斥原理 + 大数模板(跳蚤 POJ - 1091)

欢迎大家来刷我的专题:https://vjudge.net/contest/177165#problem/I 这是一道数范围很大的题,据说数据很水.首先题目的意识也就是 一个方程有解,系数都是<=m...

POJ 1091 跳蚤(经过三天的研究,终于真的把它搞明白了)

题意:一个序列有n+1个数 最后一个为m ,前面的数都不大于m ,一个跳蚤按照这个序列使其能够跳到距离起始位置左边一个单位长度的的地方,求有多少种这种数列; 一开始读题,我 就想着 这n+1个数相加...

poj1091跳蚤---------容斥

跳蚤 Time Limit: 1000MS   Memory Limit: 10000K Total Submissions: 6454   Accepted: 1...

poj 1091 跳蚤

题目大意:给你两个数n与m,n代表跳蚤共可以走n+1步,规定最后一步的步长为m,前n步的步长可以为1,2,......,m之中的任意一个。对于每一步步长,跳蚤可以选择向左跳或向右跳相应的长度,如果经过...

POJ_1091 跳蚤

题目描述 输入 输出 思路 如何计算需要去掉的组合 代码实现 附 big_int 类 备注题目描述Z城市居住着很多只跳蚤。在Z城市周六生活频道有一个娱乐节目。一只跳蚤将被请上一个高空钢丝的正中央。钢丝...
返回顶部
收藏助手
不良信息举报
您举报文章:深度学习:神经网络中的前向传播和反向传播算法推导
举报原因:
原因补充:

(最多只允许输入30个字)