POJ 1091(数论)

题目大意是给定两个整数n和m,求出长度为n+1满足条件的数列data的个数,数列的要求下:
1)1<=data[i]<=m,for1<=i<=n
2)data[n+1]=m;
3)这个n+1个数满足:存在x1,x2,...,xn,xn+1,满足x1*data[1]+x2*data[2]+...+x(n+1)*data[n+1]=1;
根据数论的知识,若存在这样的x1,x2...xn+1,则data[1],data[2]...data[n+1]的最大公约数为1

证明:若data[1],data[2]...data[n+1]满足题意,并且存在最大公约数d(为整数);则x1*data[1]+x2*data[2]+...+x(n+1)*data[n+1]的和是d的整数倍,必不等于1

其实举个例子就明白了,例如:n=2,m=360
360=3^2*2^3*5  所有不满足条件的数列,最大公约数是360质因子的乘积,只要将这些组合去掉,就是要求的答案

具体解题步骤如下:
1、求出满m的所有质因子,存入数组num
2、求出总的序列个数吗m^n
3、设t(k)表示数列最大公约数为(k个质因子乘积)的数列的个数

f=m^n-t(1)+t(2)-t(3)+..(-1)^k*t(k);
答案 = (m ^ n) - (有公因数2的n元组)- (有公因数3的n元组)- (有公因数5的n元组)+ (有公因数2,3的n元组) +(有公因数2,5的n元组) + (有公因数3,5的n元组)- (有公因数2,3,5的n元组)。这个比公式形象些
有公因数d的n元组,每个位置上有 (m/d)个选择(1 ~ m里面有m/d个d的倍数),根据乘法原理,可以得出有公因数d的n元组有 (m/d)^n 个。
ExpandedBlockStart.gifView Code
 1 #include<iostream>
 2 #include<cstdio>
 3 #include<cmath>
 4 using namespace std;
 5 __int64 n,m,per,total;
 6 __int64 s[130000],num[130000];
 7 
 8 void totalnum(__int64 x)//求质因子
 9 {
10     __int64 i;
11     total=0;
12     for(i=2;i*i<=x;i++)
13     {
14         if(x%i==0)
15         {
16             while(x%i==0) x=x/i;
17             num[total++]=i;
18         }
19     }
20     if(x!=1) num[total++]=x;
21 }
22 
23 __int64 por(__int64 x,__int64 y)//总的序列个数
24 {
25     __int64 i,k;
26     k=x;
27     for(i=1;i<y;i++)
28         x=k*x;
29     return x;
30 }
31 
32 void get(__int64 a,__int64 b,__int64 c)
33 {//a:序列起始位置 b个质因子乘积,c:公共质因子个数
34     __int64 i;
35     if(b==c)
36     {
37         __int64 t=m;
38         for(i=0;i<c;i++)
39             t=t/s[i];
40         per+=por(t,n);
41     }
42     else 
43     {
44         for(i=a;i<total;i++)
45         {
46             s[b]=num[i];
47             get(i+1,b+1,c);
48         }
49     }
50 }
51 
52 int main()
53 {
54     __int64 sum,i;
55     while(scanf("%I64d%I64d",&n,&m)!=EOF)
56     {
57         totalnum(m);
58         sum=por(m,n);
59         for(i=1;i<=total;i++)
60         {
61             per=0;
62             get(0,0,i);
63             if(i%2==1)
64                 sum-=per;
65             else sum+=per;
66         }
67         printf("%I64d\n",sum);
68     }
69     return 0;
70 }

转载于:https://www.cnblogs.com/yueshuqiao/archive/2012/02/29/2372783.html

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/262875.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

你真的了解.NET中的String吗?

你真的了解.NET中的String吗? Terrylee&#xff0c;2005年12月25日 概述 String在任何语言中&#xff0c;都有它的特殊性&#xff0c;在.NET中也是如此。它属于基本数据类型&#xff0c;也是基本数据类型中唯一的引用类型。字符串可以声明为常量&#xff0c;但是它却放在了堆中…

没有动任何配置文件,今天就出现了修改的JSP内容在页面不体现。依然是老的页面内容...

2019独角兽企业重金招聘Python工程师标准>>> 现象&#xff1a; 没有动任何配置文件&#xff0c;今天就出现了修改的JSP内容在页面不体现。依然是老的页面内容 问题分析: 既然可以用过&#xff0c;所以系统应该是检测文件更新的&#xff0c;之所以不重新编译JSP&…

DIV中class和id的区别

1.class是设置标签的类&#xff0c;id是设置标签的标识;2. 使用方法不同&#xff1a; class属性用于指定元素属于何种样式的类,如样式表可以加入.content1 { color: red; background: #ff80c0 } 使用方法&#xff1a;class"content1" id属性用于定义一个元素的独…

新浪微博api(.net)时间格式问题

Tue May 03 09:30:23 0800 2011 以前研究sinaapi的时候被这个时间格式纠结了好久&#xff0c;还自己写了一个函数去处理这种格式。最近在网上看到一个代码&#xff0c;很方便就解决了这个问题&#xff08;我就说.net怎么会没有这方法呢&#xff09;。 System.Globalization.Cul…

php代码上线,实现版本切换

以下为现有php业务&#xff0c;代码上线方式。实现&#xff1a;4套环境版本切换。[roothz-web-01 htdocs]# cat release.sh #!/usr/bin/env bashwork_dir/mnt/var/www/htdocs release_dir/var/www/htdocs# app environment app_env$1 # build release file path build_file_pat…

PWA

转载于:https://www.cnblogs.com/QianDingwei/p/11359160.html

Windows XP SP3?转自作者: 孤单的鸽子,  出处:天极yesky, 责任编辑: 袁绍龙

经过测试和检查之后&#xff0c;我们认为&#xff0c;SP3并未为Windows XP系统引入任何新的功能&#xff0c;它仅仅是一款安全升级和BUG修正补丁程序集……原文&#xff1a;http://winsystem.ctocio.com.cn/windowsxp/321/7492821.shtml目前&#xff0c;互联网上围绕Windows XP…

asp.net与JAVASCRIPT函数的相互调用

1.如何在JavaScript访问C#函数? 方法一&#xff1a; &#xff08;1&#xff09;、首先建立一个按钮&#xff0c;在后台将调用或处理的内容写入button_click中; &#xff08;2&#xff09;、在前台写一个js函数&#xff0c;内容为document.getElementById("btn1").…

汇编推荐书籍

1.《深入理解计算机系统》&#xff0c;必看&#xff01; 2.《see mips run linux》 &#xff0c; 千万不要认真看&#xff0c;看这本书的目录和小标题&#xff0c;了解mips需要学习些什么&#xff0c;而不是通过这本书来学基础知识&#xff0c;不懂的概念上网查就可以了&#x…

创建表头固定,表体可滚动的GridView(转)

本方法只需要设置一个GridView的宽度&#xff0c;其它宽度不需要设置。测试环境&#xff1a;IE6&#xff0c;Firefox通过。 C# <% Page Language"C#" AutoEventWireup"true" %> <!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional/…

3G与4G到底有何区别?

对每一个手机用户来说&#xff0c;3G和4G都是手机术语词典中最神秘的词汇&#xff0c;然而在任何卖手机和平板电脑的地方你都能发现它们的身影。当你想买一部手机的时候&#xff0c;通常不清楚该选择哪一种&#xff0c;但有一点要明确——不要随便追求高的数字。这篇文章将帮助…

海尔5D净水洗热水器引领中国制造未来

在过去的一段时间里&#xff0c;智能家居与物联网相继成为热潮&#xff0c;云计算与大数据开始被应用于不同领域的情况下&#xff0c;越来越多的智能家居产品也开始走入普通消费者生活之中……在这些独特的产业创新风起云涌的情况下&#xff0c;越来越多的行业开始推动着中国制…

Lambda表达式——注重过程的编程思想

一、使用匿名内部类的匿名对象创建线程和Lambda表达式写法 Lambda表达式写法不用去定义一个Runable接口的实现类&#xff1b; 二、方法入参是一个接口或者接口的实现类 三、对某个类的一些对象实例进行排序 其中&#xff0c;数组创建时候使用的是静态初始化&#xff0c;里面存放…

[转] STL中map用法详解

一&#xff0e;Map概述 Map是STL的一个关联容器&#xff0c;它提供一对一&#xff08;其中第一个可以称为关键字&#xff0c;每个关键字只能在map中出现一次&#xff0c;第二个可能称为该关键字的值&#xff09;的数据处理能力&#xff0c;由于这个特性&#xff0c;它…

HDU 1874 畅通工程续

畅通工程续 Time Limit: 3000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 9566 Accepted Submission(s): 3200 Problem Description某省自从实行了很多年的畅通工程计划后&#xff0c;终于修建了很多路。不过路多了也不好&am…

uboot中变量env(收集)

Env在u-boot中通常有两种存在方式&#xff0c;在永久性存储介质中&#xff08;flash、NVRAM等&#xff09;&#xff0c;在SDRAM中。可配置不适用env的永久存储方式&#xff0c;但不常用。U-boot在启动时会将存储在永久性存储介质中的env重新定位到RAM中&#xff0c;这样可以快速…

转换流

一、编码和解码 二、编码带来的问题 三、转换流InputStreamReader和OutputStreamWriter编码和解码的原理 看懂的字符变成看不懂的字节编码 和 把看不懂的字节编码变成看懂的字符&#xff1b; 字节和字符转换的桥梁&#xff0c;不写默认utf-8&#xff08;国际规定&#xff09;&a…

此情何时休

此情何时休 ——代腾飞 2008年5月6日 于成都 昨天你说今生要和我相依相守携手到白头 奈何今天你就突然说我们的缘分到了尽头 没有任何理由你就说要分手 你可知我的心里真的好难受 拉着你的手真的不想让你走 只是千言万语不知如何开口 只能轻轻松开你的手放你走 望着你远去的背…

BZOJ 3170: [Tjoi 2013]松鼠聚会 切比雪夫距离

3170: [Tjoi 2013]松鼠聚会 Time Limit: 1 Sec Memory Limit: 256 MB 题目连接 http://www.lydsy.com/JudgeOnline/problem.php?id3170Description 有N个小松鼠&#xff0c;它们的家用一个点x,y表示&#xff0c;两个点的距离定义为&#xff1a;点(x,y)和它周围的8个点即上下…

序列化流

一、对象的序列化和反系列化 二、对象的序列化流ObjectOutputStream 标记型接口&#xff1a; 使用序列化流ObjectOutputStream将对象写入到文件中进行序列化&#xff1a; 执行过后的文件&#xff1a; 二进制字节存储&#xff0c;无法直接查看&#xff1b; 三、对象的反序列化流…