【ACM-ICPC 2018 沈阳赛区网络预赛 G】Spare Tire
【链接】 我是链接,点我呀:)
【题意】
在这里输入题意
【题解】
让你求出1..n中和m互质的位置i.
让你输出∑ai
这个ai可以oeis一波。
发现是ai = i*(i+1)
1..n中和m互质的数字的个数之前有做过一题。
然后发现是逆着做的。
删掉不互质的。剩下的就是互质的了。
是用容斥原理搞的。
其中有ans+=n/temp和ass-=n/temp
表示的是加上temp,2*temp,3*temp...t*temp这些数字
以及减去temp,2*temp,3*temp...t*temp这些数字
放在这一题的话其实就是
ans+=a[temp]+a[2*temp]...+a[t*temp]
然后发现是等差的下标。那么就去推推公式吧?
然后发现真的能推出来
在get_ans2里。自己看吧。
那么ans就是下标和m不互质的ai加起来的和
然后ans=$∑_1^na_i$-ans.
【代码】
#include <bits/stdc++.h>
#define LL long long
#define rep1(i,a,b) for (int i = a;i <= b;i++)
#define rep2(i,a,b) for (int i = a;i >= b;i--)
#define all(x) x.begin(),x.end()
#define pb push_back
#define lson l,mid,rt<<1
#define ri(x) scanf("%d",&x)
#define rl(x) scanf("%lld",&x)
#define rs(x) scanf("%s",x)
#define rson mid+1,r,rt<<1|1
using namespace std;
const int N = 10000;
const LL MOD = 1e9+7;
LL n,m,pri[N+10];
LL sixni;
int num;
LL Pow(LL x,LL y){ //求x^y
LL a = 1;
while (y){
if (y&1) a = (a*x)%MOD;
x=(x*x)%MOD;
y>>=1;
}
return a;
}
LL get_ans2(LL t,LL x){
LL temp1 = t*(t+1)%MOD*(2*t+1)%MOD;
temp1 = temp1*sixni%MOD*x%MOD*x%MOD;
temp1 = temp1 + (1+t)*t/2%MOD*x%MOD;
return temp1;
}
LL get_ans(LL x)
{
LL ans = 0;
rep1(i,1,(1<<num)-1){
LL y = 1,f = 0;
rep1(j,1,num)
if (i & (1<<(j-1))){
y = y*pri[j];
f++;
}
if (f&1)
ans += get_ans2(x/y,y);
else
ans -= get_ans2(x/y,y);
}
ans = 2*(n+2)*(n+1)%MOD*n%MOD*sixni%MOD-ans;
ans = ans%MOD;
ans=(ans+MOD)%MOD;
return ans;
}
void _init(){
num = 0;
for (LL i = 2;i*i<=m;i++)
if ((m%i)==0){
pri[++num] = i;
while ((m%i)==0) m/=i;
}
if (m > 1) pri[++num] = m;
}
int main()
{
sixni=Pow(6,MOD-2);
while (~scanf("%lld%lld",&n,&m)){
_init();
printf("%lld\n",get_ans(n));
}
return 0;
}
【ACM-ICPC 2018 沈阳赛区网络预赛 G】Spare Tire的更多相关文章
- ACM-ICPC 2018 沈阳赛区网络预赛 G Spare Tire (素因子分解+容斥)
. 样例输入复制 4 4 样例输出复制 14 #include<bits/stdc++.h> using namespace std; typedef long long ll; cons ...
- ACM-ICPC 2018 沈阳赛区网络预赛 G Spare Tire(容斥)
https://nanti.jisuanke.com/t/31448 题意 已知a序列,给你一个n和m求小于n与m互质的数作为a序列的下标的和 分析 打表发现ai=i*(i+1). 易得前n项和为 S ...
- ACM-ICPC 2018 沈阳赛区网络预赛 G. Spare Tire
这题很好啊,好在我没做出来...大概分析了一下,题目大概意思就是求 问所有满足1<=i<=n且i与m互素的ai之和 最开始我们队的做法是类似线性筛的方法去筛所有数,把数筛出来后剩下数即可, ...
- ACM-ICPC 2018 沈阳赛区网络预赛 G. Spare Tire (容斥原理)
可推出$a_n = n^2+n, $ 设\(S_n = \sum_{i=1}^{n} a_i\) 则 \(S_n = \frac{n(n+1)(2n+1)}{6} + \frac{n(n+1)}{2} ...
- ACM-ICPC 2018 沈阳赛区网络预赛 G 容斥原理
https://nanti.jisuanke.com/t/31448 解析 易得an=n*n+n O(1)得到前n项和 再删除与m不互素的数 我们用欧拉函数求出m的质因数 枚举其集合的子集 进行 ...
- ACM-ICPC 2018 沈阳赛区网络预赛 J树分块
J. Ka Chang Given a rooted tree ( the root is node 11 ) of NN nodes. Initially, each node has zero p ...
- ACM-ICPC 2018 徐州赛区网络预赛 G. Trace (思维,贪心)
ACM-ICPC 2018 徐州赛区网络预赛 G. Trace (思维,贪心) Trace 问答问题反馈 只看题面 35.78% 1000ms 262144K There's a beach in t ...
- ACM-ICPC 2018 沈阳赛区网络预赛 K Supreme Number(规律)
https://nanti.jisuanke.com/t/31452 题意 给出一个n (2 ≤ N ≤ 10100 ),找到最接近且小于n的一个数,这个数需要满足每位上的数字构成的集合的每个非空子集 ...
- ACM-ICPC 2018 沈阳赛区网络预赛-K:Supreme Number
Supreme Number A prime number (or a prime) is a natural number greater than 11 that cannot be formed ...
随机推荐
- _DataStructure_C_Impl:链串
//_DataStructure_C_Impl:链串 #include<stdio.h> #include<stdlib.h> #include<string.h> ...
- Android应用开发进阶篇-场景文字识别
因为研究生毕业项目须要完毕一个基于移动终端的场景文字识别系统.尽管离毕业尚早,但出于兴趣的缘故,近一段抽时间完毕了这样一套系统. 主要的架构例如以下: client:Android应用实现拍摄场景图片 ...
- 优雅的App全然退出方案(没有不论什么内存泄漏隐患)
在Android开发过程中,特别是界面比較多的情况下,用寻常的退出方式往往是不能全然退出这个应用,网络上也好多各种退出方案.当中一种应该是被广大开发人员採纳使用,也很的清晰方便.就是在Applicat ...
- Cocos2d-X直接使用OpenGL接口
Cocos2d-X是基于基于OpenGL ES的2D游戏引擎,所以Cocos2d-X能够直接使用OpenGL接口 首先建立一个Draw类,用于处理OpenGL接口 在Draw.h中加入以下的代码 #i ...
- 坚向的ViewPager,上下滑动的组件,android上下滑动 VerticalPager
package com.zhulin.android.atools; import android.content.Context; import android.util.AttributeSet; ...
- CodeForces 651A(水题)
Friends are going to play console. They have two joysticks and only one charger for them. Initially ...
- poj1700--贪心--Crossing River
Crossing River Time Limit: 1000MS Memory Limit: 10000K Total Submissions: 12260 Accepted: 4641 D ...
- 自定义数据类型 C++ 结构体类型 共同体类型 枚举类型 类类型{}
一.结构体类型 结构体类型,共用体类型,枚举类型,类类型等统称为自定义类型(user-defined-type,UDT). 结构体相当于其他高级语言中的记录(record);例如: struct St ...
- yii1 session
在 Yii框架中使用session 的笔记: 首先,在Yii框架中,你不需要像标准PHP代码那样使用session_start(),在Yii框架中,autoStart 属性缺省被设置为true,所以, ...
- shp系列(四)——利用C++进行Shx文件的读(打开)
1.shx文件的基本情况 shx文件又叫索引文件,主要包含坐标文件的索引信息,文件中每个记录包含对应的坐标文件记录距离坐标文件的初始位置的偏移量.通过索引文件可以很方便地在坐标文件中定位到指定目标的坐 ...