【欧拉函数】【HDU1286】 找新朋友
找新朋友
Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others)
Total Submission(s): 7651 Accepted Submission(s): 4033
2
25608
24027
7680
16016
求约数 直接筛的裸算法过了。。
查题解发现有欧拉函数这个东西
在数论中,对正整数n,欧拉函数是小于或等于n的正整数中与n互质的数的数目。此函数以其首名研究者欧拉命名,它又称为φ函数、欧拉商数等。
p2……pn为x的所有质因数
因此 先打个素数表 求他是那些素数 然后直接用欧拉公式即可;
1.用约数筛
#include<stdio.h>
#include<math.h>
#include<string.h>
int yueshu[200];
int yuenum;
int ok[40000];
int ans=0;
int main()
{
freopen("a.in","r",stdin);
freopen("a.out","w",stdout);
int CN,i,N,mem,j,t,k;
while(scanf("%d",&CN)!=EOF)
{
for(i=1;i<=CN;i++)
{
memset(yueshu,0,sizeof(yueshu));
memset(ok,0,sizeof(ok));
yuenum=1;
ans=0;
scanf("%d",&mem);
t=sqrt(mem);
for(j=1;j<=t;j++)
if(mem%j==0)
{
yueshu[yuenum++]=j;
yueshu[yuenum++]=mem/j;
}
yuenum--;
for(j=2;j<=yuenum;j++)
{
for(k=1;k*yueshu[j]<=mem;k++)
ok[k*yueshu[j]]=1;
}
for(j=1;j<=mem;j++)
if(ok[j]==0) ans++;
printf("%d\n",ans);
}
}
return 0;
}
2.O(n)级别素数筛选+欧拉
#include<stdio.h>
int YNprime[40001];
int prime[40000];
int totprime=1;
int get_prime(int maxn)
{
int i,j;
for(i=2;i<=maxn;i++)
{
if(YNprime[i]==0) prime[totprime++]=i;
for(j=1;i*prime[j]<=maxn&&j<totprime;j++)
{
YNprime[i*prime[j]]=1;
if(i%prime[j]==0) break;
}
}
}
int main()
{
freopen("a.in","r",stdin);
freopen("a.out","w",stdout);
get_prime(40000);
int CN,N,i,j,ans;
while (scanf("%d",&CN)!=EOF)
{
for(i=1;i<=CN;i++)
{
scanf("%d",&N);
ans=N;
if(N==1) { printf("0\n");continue;}
for(j=1;prime[j]<=N&&j<totprime;j++)
{
if(N%prime[j]==0)
ans=(ans/prime[j])*(prime[j]-1);
}
printf("%d\n",ans);
} }
return 0;
}
但有个更优美的代码:
用最小的素因子筛掉每个数
int prime[N],phi[N],cnt;// prime:记录质数,phi记录欧拉函数
int Min_factor[N];// i的最小素因子
bool vis[N];
void Init()
{
cnt=0;
phi[1]=1;
int x;
for(int i=2;i<N;i++)
{
if(!vis[i])
{
prime[++cnt]=i;
phi[i]=i-1;
Min_factor[i]=i;
}
for(int k=1;k<=cnt&&prime[k]*i<N;k++)
{
x=prime[k]*i;
vis[x]=true;
Min_factor[x]=prime[k];
if(i%prime[k]==0)
{
phi[x]=phi[i]*prime[k];
break;
}
else phi[x]=phi[i]*(prime[k]-1);
}
}
}
【欧拉函数】【HDU1286】 找新朋友的更多相关文章
- 欧拉函数之HDU1286找新朋友
找新朋友 Time Limit : 2000/1000ms (Java/Other) Memory Limit : 65536/32768K (Java/Other) Total Submissi ...
- hdu1286 找新朋友 欧拉函数模板
首先这一题用的是欧拉函数!!函数!!不是什么欧拉公式!! 欧拉函数求的就是题目要求的数. 关于欧拉函数的模板网上百度一下到处都是,原理也容易找,这里要介绍一下另一个强势模板. 在这一题的讨论里看到的. ...
- hdu1286(找新朋友)&&POJ2407Relatives(欧拉函数模版题)
http://acm.hdu.edu.cn/showproblem.php?pid=1286 没什么好说的,模板题,主要是弄懂欧拉函数的思想. #include <iostream> #i ...
- hdu-1286 找新朋友(欧拉函数,水题)
题目链接: 找新朋友 Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others)Total ...
- hdu1286 找新朋友
找新朋友 http://acm.hdu.edu.cn/showproblem.php?pid=1286 Time Limit: 2000/1000 MS (Java/Others) Memory ...
- 欧拉函数 - HDU1286
欧拉函数的作用: 有[1,2.....n]这样一个集合,f(n)=这个集合中与n互质的元素的个数.欧拉函数描述了一些列与这个f(n)有关的一些性质,如下: 1.令p为一个素数,n = p ^ k,则 ...
- 找新朋友---hdu1286(欧拉函数)
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1286 欧拉函数:对正整数n,欧拉函数是求少于n的数中与n互质的数的数目: 素数(质数)指在一个大于1的 ...
- hdu 1286:找新朋友(数论,欧拉函数)
找新朋友 Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others)Total Submis ...
- hdoj 1286 找新朋友【欧拉函数】
找新朋友 Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others)Total Submis ...
随机推荐
- java集合总结【转】
Map.Set.Iterator迭代详解 Map接口定义了四种类型的方法,每个Map都包含这些方法. equals(Object o)比较指定对象与此Map的等价性. hashCode()返回此Map ...
- Oracle11g x64使用Oracle SQL Developer报错:Unable to find a Java Virtual Machine
原因oracle 11g中安装的Oracle SQL Developer是32位的,而我们现在给他指定的java.exe却是64位的,所以会出现这种错误.解决方法1)从网上下载Oracle SQL D ...
- oracle11gr2 netca 无法启动 报错
安装oracle软件后,必须要先配置listener才能dbca建库,但是netca却报下面的错误. Oracle Net Services Configuration: # # An unexpec ...
- ios 动态修改UILabel字体大小
- (IBAction)sliderChange:(id)sender { NSLog(@"sliderChange"); UISlider *slider = (UISl ...
- Gengxin讲STL系列——Queue和Stack
第三篇. 感觉队列和栈是必须的……所以决定加上这两个…… 我发现我已经买域名买上隐了……今天又买了个.top……真是智障…… Queue(队列FIFO)和Statk(栈LIFO). 那么为什么要这两个 ...
- (转) 学习C++ -> 指针初步
学习C++ -> 指针初步 一.指针 1. 什么是指针? 我们知道, 计算机的内存是由一个个独立的存储单元组成, 并且系统会对每一个存储单元分配一个唯一的号码, 称为这个存储 ...
- POJ 1698 (二分图的多重匹配)
转载:http://www.cppblog.com/MatoNo1/archive/2011/03/26/142766.aspx 我们知道在一个图中,每个点最多只能匹配一条边的情况,是二分图的最大匹配 ...
- C++程序设计实践指导1.10二维数组元素换位改写要求实现
改写要求1:改写为以单链表和双向链表存储二维数组 改写要求2:添加函数SingleLinkProcess()实现互换单链表中最大结点和头结点位置,最小结点和尾结点位置 改写要求3:添加函数Double ...
- c#一个分页控件的例子
一.首先下载一个dll,地址:http://pan.baidu.com/share/link?shareid=1628211605&uk=1342867987 二.添加到项目中 三.添加引用 ...
- 工作中小知识点汇总(c#)
1.OOP 实体与数据库字段转换(注意 此时实体字段必须和数据库中查询的字段列名相同) list = ModelExtend.GetByDataTablePart<EZRate>(ds.T ...