徐州赛区网络预赛 D Easy Math
比赛快结束的适合看了一下D题,发现跟前几天刚刚做过的HDU 5728 PowMod几乎一模一样,当时特兴奋,结果一直到比赛结束都一直WA。回来仔细一琢磨才发现,PowMod这道题保证了n不含平方因子,而Easy Math却没有。只要加一条特判,这道题就过了。
AC代码贴在下面。求莫比乌斯函数前缀和的那一部分是我当时在网上临时找的代码,向那位仁兄道谢。
#include<bits/stdc++.h>
using namespace std;
#define rep(i,a,b) for(int i=a;i<=b;++i)
#define ms(arr,a) memset(arr,a,sizeof arr)
#define debug(x) cout<<"< "#x" = "<<x<<" >"<<endl
typedef long long ll;
const int INF = 0x3f3f3f3f;
const int Maxn = 1e7 * 2;
const int mod = 2333333;
int pcnt = 0, prime[1300000]; // 质数
short mu[Maxn]; // 莫比乌斯函数值
bool vis[Maxn];
void init() {
mu[1] = 1;
for (int i = 2; i < Maxn; i++) {
if (vis[i] == 0) {
mu[i] = -1;
prime[++pcnt] = i;
}
for (int j = 1; j <= pcnt && i * prime[j] < Maxn; j++) {
vis[i * prime[j]] = 1;
if (i % prime[j] != 0)
mu[i * prime[j]] = -mu[i];
else {
mu[i * prime[j]] = 0;
break;
}
}
}
for (int i = 2; i < Maxn; i++) mu[i] += mu[i - 1]; // 函数前缀和
}
struct Hash {
long long key;
int value, next;
} node[mod];
int cnt = 0, Head[mod] = {0};
void Insert(long long N, int v) { // 记忆
int ha = N % mod;
node[++cnt].key = N;
node[cnt].value = v;
node[cnt].next = Head[ha];
Head[ha] = cnt;
}
int M(long long N) {
if (N < Maxn) return mu[N];
int ha = N % mod;
for (int i = Head[ha]; i != 0; i = node[i].next) {
if (node[i].key == N) // 如果已经计算过
{return node[i].value;}
}
int ans = 0;
for (long long i = 2, j; i <= N; i = j + 1) {
j = N / (N / i); // 分块加速
ans += (j - i + 1) * M(N / i);
}
Insert(N, 1 - ans); // 做记忆
return 1 - ans;
}
ll n,m;
long long fac[30];int cntfac;
bool init_fac(ll x)
{
cntfac=0;
for(int i=1;prime[i]<sqrt(x);++i)
if(x%prime[i]==0){fac[++cntfac]=prime[i];x/=prime[i];if(x%prime[i]==0)return false;}
if(x>1)fac[++cntfac]=x;
return true;
}
ll calc(ll n,ll m,int k)
{
//debug(n);
if(m==0)return 0;
if(n==1)
{
if(m==0)return 0;
else {return M(m);}
}
return -calc(n/fac[k],m,k-1)+calc(n,m/fac[k],k);
}
int main()
{
//freopen("Input.txt","r",stdin);
//freopen("1.txt","w",stdout);
init();
while(~scanf("%lld%lld",&m,&n))
{
if(init_fac(n))printf("%lld\n",calc(n,m,cntfac));
else printf("0\n");
}
//freopen("CON","w",stdout);
//system("start Output.txt");
}
徐州赛区网络预赛 D Easy Math的更多相关文章
- [ACM-ICPC 2018 徐州赛区网络预赛][D. Easy Math]
题目链接:Easy Math 题目大意:给定\(n(1\leqslant n\leqslant 10^{12}),m(1\leqslant m\leqslant 2*10^{9})\),求\(\sum ...
- ACM-ICPC 2018 徐州赛区网络预赛(8/11)
ACM-ICPC 2018 徐州赛区网络预赛 A.Hard to prepare 枚举第一个选的,接下来的那个不能取前一个的取反 \(DP[i][0]\)表示选和第一个相同的 \(DP[i][1]\) ...
- 计蒜客 1460.Ryuji doesn't want to study-树状数组 or 线段树 (ACM-ICPC 2018 徐州赛区网络预赛 H)
H.Ryuji doesn't want to study 27.34% 1000ms 262144K Ryuji is not a good student, and he doesn't wa ...
- ACM-ICPC 2018 徐州赛区网络预赛 I. Characters with Hash
Mur loves hash algorithm, and he sometimes encrypt another one's name, and call him with that encryp ...
- ACM-ICPC 2018 徐州赛区网络预赛 B(dp || 博弈(未完成)
传送门 题面: In a world where ordinary people cannot reach, a boy named "Koutarou" and a girl n ...
- ACM-ICPC 2018 徐州赛区网络预赛 B. BE, GE or NE
In a world where ordinary people cannot reach, a boy named "Koutarou" and a girl named &qu ...
- ACM-ICPC 2018 徐州赛区网络预赛 H. Ryuji doesn't want to study
262144K Ryuji is not a good student, and he doesn't want to study. But there are n books he should ...
- ACM-ICPC 2018 徐州赛区网络预赛 F. Features Track
262144K Morgana is learning computer vision, and he likes cats, too. One day he wants to find the ...
- ACM-ICPC 2018 徐州赛区网络预赛 D 杜教筛 前缀和
链接 https://nanti.jisuanke.com/t/31456 参考题解 https://blog.csdn.net/ftx456789/article/details/82590044 ...
随机推荐
- lly的数列询问(最小生成树 + 思维)
lly的数列询问 Description 这个问题很简单,lly lly lly给你一些提示,让你试着确定长度为n n n的数列A [1] A [2] ... A [n]的值,但是他想尽一切办法为大家 ...
- Java并发基础10:原子性操作类的使用
在 java5 以后,我们接触到了线程原子性操作,也就是在修改时我们只需要保证它的那个瞬间是安全的即可,经过相应的包装后可以再处理对象的并发修改,本文总结一下Atomic系列的类的使用方法,其中包含: ...
- vscode 保存自动 格式化eslint 代码
在网上找了很多种方法,大多都没有成功 一下是一种成功的 配置方法: 1) First, you need to install the ESLint extension in the VS code ...
- Go语言笔记 (2) 变量命名与多重赋值
变量命名 1.大小写 观摩以下代码: func main() { var m int = "你" var M int = "我" fmt.Println(m,M ...
- Jmeter 压力测试笔记(1)--服务器迁移失败
近期,公司服务器因技术架构升级等原因需要迁移,在经过开发,运维DBA,测试多部门进行联合讨论后,制定出了迁移方案.迁移前也对APP应用进行了各种测试,并没有发现问题. 凌晨2点开始迁移,5点完成迁移. ...
- python 面向对象反射以及内置方法
一.反射 什么是反射:可以用字符串的方式去访问对象的属性,调用对象的方法(但是不能去访问方法),python中一切皆对象,都可以使用放射. 反射的四种方法: hasattr:hasattr(objec ...
- Github上面拉取别人提交的PR
在github上面协同开发,避免不了拉取别的同学的PR,那么如何拉取呢? 1.首先我们看下upstream liz@liz-PC:~/jimeng/handle-api$ git remote -v ...
- iOS 头条一面 面试题
1.如何高效的切圆角? 切圆角共有以下三种方案: cornerRadius + masksToBounds:适用于单个视图或视图不在列表上且量级较小的情况,会导致离屏渲染. CAShapeLayer+ ...
- ClickOnce的安装路径
win 7下C:\Users\Administrator.U5G4L4PUY34SH5C\AppData\Local\Apps\2.0\KPVZOAYK.0JE\56B55RCH.A7A\winr.. ...
- Linux常用命令02(远程管理)
01 关机/重启 序号 命令 对应英文 作用 01 shutdown 选项 时间 shutdown 关机/重新启动 1.1 shutdown shutdown 命令可以 安全 关闭 或者 重新启动系统 ...