[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;
}

版权声明:博客已搬迁至http://bill.moe

poj1091(容斥原理+质因子分解)跳蚤

跳蚤 Time Limit: 1000MS   Memory Limit: 10000K Total Submissions: 10520   Accepted: ...
  • update7
  • update7
  • 2017年07月24日 09:55
  • 16191

POJ1091 跳蚤 素因子分解+容斥原理

题目链接:http://poj.org/problem?id=1091 分析:我们假设卡片上的标号分别是a1,a2,...,an,M,跳蚤跳对应标号的卡片的次数分别为x1,x2,...,xn...
  • AC_Gibson
  • AC_Gibson
  • 2015年06月17日 10:38
  • 1318

基于Android的校园跳蚤市场(二手)的设计与实现

基于Android的校园跳蚤市场(二手)的设计与实现一、摘要当今移动互联网高速发展,本文详细介绍了基于Android的跳蚤市场app从可行性研究,需求分析,总体设计,详细设计,实现、编码、测试与维护的...
  • qq137722697
  • qq137722697
  • 2016年08月28日 23:21
  • 17046

HNOI2002跳蚤--容斥原理

题目描述:读入n,m两个数,假设一个合法的数列是n+1位,且前n位不大于m,第n+1位为m。对于每一个数列,跳蚤可以选取任意一个数列中的数k,往左或右走k步(可以走多次),若使用这一个数列跳蚤可以到达...
  • henyouyisiba123
  • henyouyisiba123
  • 2016年10月16日 21:25
  • 302

poj 1091 跳蚤(n元一次不定方程+斥容原理)

题意: 给一张卡片,上面写着N + 1个自然数,其中最后一个数是M,并且前N个数不超过M,并允许数字相同。 一只跳蚤每次可以选择从卡片上任意选择一个自然数S,向左或向右跳S个单位长度。 求最后能...
  • u013508213
  • u013508213
  • 2015年01月05日 13:29
  • 358

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

题目:http://poj.org/problem?id=1091 跳蚤 Time Limit: 1000MS   Memory Limit: 10000K Total Submiss...
  • theArcticOcean
  • theArcticOcean
  • 2015年08月29日 20:14
  • 718

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

跳蚤 Time Limit: 1000MS   Memory Limit: 10000K Total Submissions: 8129   Accepted: 2...
  • u011836218
  • u011836218
  • 2014年05月14日 15:49
  • 742

poj 1091-跳蚤(数论+容斥)

题目链接 题意: 可转化为已知N+1个数,x1, x2, .. xn, m;使得gcd(x1, x2, ..., xn) = 1的解的个数(xi 思路:m^n为所有不同序列个数,然后减去公因数>1...
  • u014492306
  • u014492306
  • 2015年09月28日 16:39
  • 209

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

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

POJ 1091 跳蚤 数论-容斥原理、扩展欧几里得

题目地址: http://poj.org/problem?id=1091 【题目大意】 输入整数n和m,找n个数字(都小于等于m),加上m共n+1个数字,使得这个n+1个数字满足存在x1…xn+1...
  • ilovexiaohao
  • ilovexiaohao
  • 2013年07月16日 10:10
  • 853
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:[HNOI2002][poj1091]跳蚤
举报原因:
原因补充:

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