【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 ...
随机推荐
- ListViewItem中的图片不能动态改变的解决方法
近期遇到了一个问题,就是我的listviewitem中有个图片,点击的时候须要变成还有一种图片.结果在getView()中设置了响应.可是能够运行.就是不起作用.在网上查了非常多资料也没有解决.最后发 ...
- 零基础学python-2.16 列表解析
这一节聊聊强大的列表解析 主要就是在一行里面赋值给列表 以下我们举两个样例: 上面的样例我们引入了range函数,他主要作用是在一定范围里面取整数值 我来解释一下中括号中面的那一句:x**2 for ...
- unity游戏开发之NGUI的UISprite染色
游戏的UI开发中常常会遇到染色问题.比如button失效变灰的效果,同一个道具通过策划表配的颜色值染上红绿蓝紫等颜色,效果例如以下 最笨最挫的方法当然是让美术多出几个资源图.这种一个缺点是浪费资源,在 ...
- C#之简易猜数字游戏
游戏规则说明: 由系统生成一个随机数,玩家有三次猜数字的机会,如果在三次内猜出数字反馈玩家猜对了,否则Game Over! 代码设计说明: 1.首先设计一个简易的欢迎界面,并提示玩家是否开始游戏: p ...
- oc10--练习
// // main.m // 练习 #import <Foundation/Foundation.h> @interface Car : NSObject { @public int w ...
- Node.js:工具模块
ylbtech-Node.js:工具模块 1.返回顶部 1. Node.js 工具模块 在 Node.js 模块库中有很多好用的模块.接下来我们为大家介绍几种常用模块的使用: 序号 模块名 & ...
- php添加文件到压缩文件夹
extension=php_zip.dll注释掉 <?php$filename = "test.zip"; ob_end_clean(); $zip = new ZipArc ...
- 关于一些UI的插件(杂)
1.时间插件 //日期框 $('.date-picker').datepicker(); 2.checkbox 保存checkbox的值 // 组装选择的标签 var check = $(" ...
- SqlServer备份恢复
备份 (Transact SQL) | Microsoft Docs https://docs.microsoft.com/zh-cn/sql/t-sql/statements/backup- ...
- python os 模块常用操作
python 2.7 os 常用操作 官方document链接 文件和目录 os.access(path, mode) 读写权限测试 应用: try: fp = open("myfile&q ...