2017.10.7 phi的反函数 失败总结
2017-10-07 16:56
405 查看
数论题好难啊,,分析也好,爆搜也好,剪枝也好,什么都不会啊
一开始能想到的就只有根据单点求phi来进行反推。。。
然后并不能保证最小解,,而且搜索意味不明。
其实需要注意到任意的phi值 全部都是可以写成素数的phi乘积 就好做了
就可以枚举素数,作为phi反函数的一个因子,对当前phi值的影响就是/ 素数-1
注意: phi(81)=phi(3)*3*3;;
码:
#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;
#define N 50000
#define ll long long
int he
,su
,tot,i,j,k;
ll daan,n;
void eular(int o)
{
for(i=2;i<=o;i++)
{
if(!he[i])
{
su[++tot]=i;
}
for(j=1;j<=tot&&i*su[j]<=o;j++)
{
he[i*su[j]]=1;
if(i%su[j]==0)break;
}
}
}
bool shisu(int o)
{
int i;
for(i=1;su[i]<=sqrt(o);i++)
{
if(o%su[i]==0)return 0;
}
return 1;
}
void dfs(int last,ll ans,int qi)
{
if(ans>=daan)return;
if(last==1){daan=ans;return;}
if(last>sqrt(n)&&shisu(last+1)){ daan=min(daan,ans*(last+1)); return;}
int i;
for(i=qi+1;su[i]-1<=sqrt(n);i++)
{
if(su[i]-1>last)break;
if(!(last%(su[i]-1)))
{
int x=last/(su[i]-1); ll lin=ans*(su[i]); dfs(x,lin,i);
while(!(x%(su[i])))
{
x/=(su[i]);
lin*=(su[i]);
dfs(x,lin,i);
}
}
}
}
int main()
{
scanf("%lld",&n);
eular(sqrt(2147483647));
daan=2147483648;
dfs(n,1,0);
if(daan==2147483648)printf("-1");
else printf("%lld",daan);
}
一开始能想到的就只有根据单点求phi来进行反推。。。
然后并不能保证最小解,,而且搜索意味不明。
其实需要注意到任意的phi值 全部都是可以写成素数的phi乘积 就好做了
就可以枚举素数,作为phi反函数的一个因子,对当前phi值的影响就是/ 素数-1
注意: phi(81)=phi(3)*3*3;;
码:
#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;
#define N 50000
#define ll long long
int he
,su
,tot,i,j,k;
ll daan,n;
void eular(int o)
{
for(i=2;i<=o;i++)
{
if(!he[i])
{
su[++tot]=i;
}
for(j=1;j<=tot&&i*su[j]<=o;j++)
{
he[i*su[j]]=1;
if(i%su[j]==0)break;
}
}
}
bool shisu(int o)
{
int i;
for(i=1;su[i]<=sqrt(o);i++)
{
if(o%su[i]==0)return 0;
}
return 1;
}
void dfs(int last,ll ans,int qi)
{
if(ans>=daan)return;
if(last==1){daan=ans;return;}
if(last>sqrt(n)&&shisu(last+1)){ daan=min(daan,ans*(last+1)); return;}
int i;
for(i=qi+1;su[i]-1<=sqrt(n);i++)
{
if(su[i]-1>last)break;
if(!(last%(su[i]-1)))
{
int x=last/(su[i]-1); ll lin=ans*(su[i]); dfs(x,lin,i);
while(!(x%(su[i])))
{
x/=(su[i]);
lin*=(su[i]);
dfs(x,lin,i);
}
}
}
}
int main()
{
scanf("%lld",&n);
eular(sqrt(2147483647));
daan=2147483648;
dfs(n,1,0);
if(daan==2147483648)printf("-1");
else printf("%lld",daan);
}
相关文章推荐
- Visual studio 2013安装失败总结
- Ajax前后台交互的失败总结
- 人人旗下风车网CEO的创业失败教训总结
- json_decode()解析json数据失败原因总结
- 人人旗下风车网CEO的创业失败教训总结
- 项目1失败总结
- 创业失败的10个教训总结
- Spring框架:跨域问题之使用@CrossOrigin注解解决失败的原因总结
- openstack---失败总结、、、
- 早起习惯养成失败总结
- 虚拟机克隆网卡启动失败总结
- Atitit.java c#这类编程语言的设计失败点attilax总结
- linux 3.4.103 内核移植到 S3C6410 开发板 移植失败 (问题总结,日本再战!)
- 2017.3.19 约数个数和 失败总结
- 2017.9.27 可持久化并查集 失败总结
- 失败——阿里校招内推电话面试总结
- 这个学期总结 失败啊
- 一个电子商务网站的失败总结
- 2017.9.2 王室联邦 失败总结
- 关于TransactionScope出错:“与基础事务管理器的通信失败”的解决方法总结