HDU1299--Diophantus of Alexandria HDU(109)
2013-08-15 21:26
357 查看
看到这题的时候想到了刘汝佳里面的原题,就直接暴力枚举一下。
完全的TLE了。好吧,其实测试了一下99999999,发现果然过不了。虽然只是遍历了a~2*a但是好像超了范围。
然后试了试定义i为long long或者__int64,居然直接蹦了,你妹了。
索性看了看题解,发现大家代码基本上都一模一样(天下文章一大抄)。
看了看思路,发现原来不用这么做,直接用数论的素因子定理就好。于是明白了什么..............直接开工。
但是做到最后发现为什么q*3没明白,索性去查查
分解素因子时,如果n没法完全分解(n包含大于sqrt(n)素因子)sum*=3最后在sum=(sum+1)/2
在一个大牛的博客发现了这句话。
这就明白了,所以就A了...............
还是数论好玩。
#pragma warning(disable:4786)
#include<iostream>
#include<algorithm>
#include<cmath>
#include<stdio.h>
#include<time.h>
#include<stdlib.h>
#include<queue>
#include<set>
#include<vector>
#include<string>
#include<ctime>
#include<string.h>
using namespace std;
#define pi acos(-1.0)
//#define LL __int64
typedef long long LL;
#define INF 0x7fffffffffffffff
#define bug puts("hear!")
#define inf 0x7fffffff
#define eps 1e-10
#define FRE freopen("in.txt","r",stdin)
#define E exp(1.0)
#define mod 1000000007
const int N=40000;
int s
,m
,l;
void is_prime()
{
int i,j;
for (i=2;i<=N;i++)
{
for (j=2;i*j<N;j++)
{
s[i*j]=0;
}
}
l=0;
for (i=2;i<N;i++)
{
if(s[i])
m[l++]=i;
}
}
int main()
{
int t,n,l,num;
memset(s,1,sizeof(s));
is_prime();
while (scanf("%d",&t)!=EOF)
{
int z=1;
while (t--)
{
scanf("%d",&n);
int q=1;
for(int i=0;i<l;i++)
{
int o=(int)sqrt(n*1.0)+1;
if(m[i]>o)
break;
num=0;
while (n%m[i]==0)
{
num++;
n/=m[i];
}
q*=(2*num+1);
}
if(n>1)
q*=3;
printf("Scenario #%d:\n",z++);
printf("%d\n\n",(q+1)/2);
}
}
return 0;
}
#pragma warning(disable:4786) #include<iostream> #include<algorithm> #include<math.h> #include<stdio.h> #include<time.h> #include<stdlib.h> #include<queue> #include<set> #include<vector> #include<string> #include<ctime> #include<string.h> using namespace std; #define pi acos(-1.0) //#define LL __int64 typedef long long LL; #define INF 0x7fffffffffffffff #define bug puts("hear!") #define inf 0x7fffffff #define eps 1e-10 #define FRE freopen("in.txt","r",stdin) #define E exp(1.0) #define mod 1000000007 int gcd(int a,int b) { return b==0?a:gcd(b,a%b); } int main(){ int n,l; l=1; scanf("%d",&n); while(n--) { int a,q; q=0; scanf("%d",&a); for(int i=a;i<=2*a;i++) { int z,f,g; z=i/gcd(i,a)*a; f=z/a-z/i; g=z%f; if(f==1) q++; } printf("Scenario #%d:\n",l); printf("%d\n\n",q); l++; } return 0; }
完全的TLE了。好吧,其实测试了一下99999999,发现果然过不了。虽然只是遍历了a~2*a但是好像超了范围。
然后试了试定义i为long long或者__int64,居然直接蹦了,你妹了。
索性看了看题解,发现大家代码基本上都一模一样(天下文章一大抄)。
看了看思路,发现原来不用这么做,直接用数论的素因子定理就好。于是明白了什么..............直接开工。
但是做到最后发现为什么q*3没明白,索性去查查
分解素因子时,如果n没法完全分解(n包含大于sqrt(n)素因子)sum*=3最后在sum=(sum+1)/2
在一个大牛的博客发现了这句话。
这就明白了,所以就A了...............
还是数论好玩。
#pragma warning(disable:4786)
#include<iostream>
#include<algorithm>
#include<cmath>
#include<stdio.h>
#include<time.h>
#include<stdlib.h>
#include<queue>
#include<set>
#include<vector>
#include<string>
#include<ctime>
#include<string.h>
using namespace std;
#define pi acos(-1.0)
//#define LL __int64
typedef long long LL;
#define INF 0x7fffffffffffffff
#define bug puts("hear!")
#define inf 0x7fffffff
#define eps 1e-10
#define FRE freopen("in.txt","r",stdin)
#define E exp(1.0)
#define mod 1000000007
const int N=40000;
int s
,m
,l;
void is_prime()
{
int i,j;
for (i=2;i<=N;i++)
{
for (j=2;i*j<N;j++)
{
s[i*j]=0;
}
}
l=0;
for (i=2;i<N;i++)
{
if(s[i])
m[l++]=i;
}
}
int main()
{
int t,n,l,num;
memset(s,1,sizeof(s));
is_prime();
while (scanf("%d",&t)!=EOF)
{
int z=1;
while (t--)
{
scanf("%d",&n);
int q=1;
for(int i=0;i<l;i++)
{
int o=(int)sqrt(n*1.0)+1;
if(m[i]>o)
break;
num=0;
while (n%m[i]==0)
{
num++;
n/=m[i];
}
q*=(2*num+1);
}
if(n>1)
q*=3;
printf("Scenario #%d:\n",z++);
printf("%d\n\n",(q+1)/2);
}
}
return 0;
}
相关文章推荐
- Hdu 1299 Diophantus of Alexandria
- HDU 1299 Diophantus of Alexandria (数论)
- HDU 1299 Diophantus of Alexandria(数论)
- HDU 1299 Diophantus of Alexandria (公式变形 分解质因数)
- hdu-1299 Diophantus of Alexandria(分解素因子)
- HDU 1299 Diophantus of Alexandria
- 【HDU】 1299 Diophantus of Alexandria
- hdu 1299 Diophantus of Alexandria(数学题)
- hdu 1299 Diophantus of Alexandria (素数筛选)
- Hdu 1299 (Diophantus of Alexandria)
- hdu1299 Diophantus of Alexandria
- hdu 1299 Diophantus of Alexandria (数论)
- hdu 1299 Diophantus of Alexandria
- HDU-1299 Diophantus of Alexandria(数论:因子数)
- HDU 1299 Diophantus of Alexandria
- Diophantus of Alexandria HDU - 1299(数论)
- HDU1299-Diophantus of Alexandria
- HDU_1299 Diophantus of Alexandria
- HDU 1299 Diophantus of Alexandria 求因子个数
- HDU 1299 Diophantus of Alexandria