【GDOI2016模拟3.16】幂(容斥 + 模型复杂转化)
【GDOI2016模拟3.16】幂
\(X\in[1,A],Y\in[1,B]\),问:\(x^y\)的不用取值个数.
\(A,B\)都是\(10^9\)级别.
然后我们开搞.
首先,假设一个合法的\(x\)可以表示为\(x=\prod p_i^{q_i}\),那么令\(d=gcd(q_1,q_2...q_k)\)
假设\(d>1\),显然我们不需要单独考虑,因为它可以继续化简,我们找到最简的那个数然后去一次性处理.
那么此时所有情况都变成了\(d=1\).
此时再分两种情况讨论,因为我们现在实际上把问题转化为了:
\(x\in[1,A]\),\(y\in [1,B]\),求\((c^x)^y\)的不同个数.
然后这个问题显然当\(x>\sqrt{A}\)时贡献就是\(B\)。
所以题目转化为:
在\(x\in[1,\sqrt{A}],y\in[1,B]\)时,\((c^x)^y\)的不同个数.
于是我们可以正难则反,计算一下不合法的个数,假设\(k={log_x}^n\),那么问题又可以转化为:
\(x\in[1,k],y\in[1,B]\),\(xy\)的相同个数.
于是这个问题值得思考一下.
咱有一个经典套路:分块计算;
即把总范围\([1,KB]\)变成\(K\)块,每一块都长的像这样的形式:
\]
那么现在单独考虑每一块的贡献.
然后假设我们现在处理第\(i\)个块,那么如果其中某一个数可以被表示为\(x*y\)的形式.
则必定存在一个\(d\in [i,k]\),满足\(d|x*y\).
因为这样子就必定能满足\(x\in [1,K],y\in [1,B]\)这两个条件.
所以现在题目再次被转化:
等价于我们要求在区间\([(i-1)B+1,iB]\)范围内有多少个数是\([i,k]\)区间中某一个数的倍数.
这次转化以后我们发现问题就变得熟悉了...
因为\(B\)比较大,\(k\)貌似小,所以,,,, 这TM容斥一下不就好了??
但是,我们发现,\(k\)最大可能为\(30\),这意味着,\(2^{30}\)是接受不了的...
但是我们很容易想到,如果当乘上一个数的lcm没有变化之后,就不需要去弄.
这样子的话时间一下子就降了下来,反正怎么打都能过~
#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#define ll long long
#define F(i, a, b) for (ll i = a; i <= b; i ++)
#define sqr(x) ((x) * (x))
const ll N = 5e5;
ll ans, sum;
ll A, B, up, P[N];
bool vis[N];
using namespace std;
void Init() {
#ifndef ONLINE_JUDGE
freopen("data.in", "r", stdin);
#endif
scanf("%lld%lld", &A, &B);
up = int(sqrt(A)) + 1;
ans = 1LL * (A - 1) * B + 1;
}
void GetPrime() {
F(i, 2, N - 1) {
if (!vis[i])
P[++ P[0]] = i;
F(j, 1, P[0]) {
if (P[j] * i >= N) break;
vis[P[j] * i] = 1;
if (i % P[j] == 0) break;
}
}
}
ll gcd(ll x, ll y){
return y == 0 ? x : gcd(y, x % y);
}
void dfs(ll x, ll v, ll st, ll en, ll g, ll k) {
if (x > k) {
sum += 1LL * (en / v - st / v) * g;
return;
}
if (x != gcd(v, x)) {
dfs(x + 1, v * x / gcd(v, x), st, en, g * (- 1), k);
dfs(x + 1, v, st, en, g, k);
}
}
ll Solve(ll x) {
ll cnt = -1;
for (ll tmp = 1; tmp <= A; tmp *= x, cnt ++);
sum = 0;
F(i, 2, cnt)
dfs(i, 1, (i - 1) * B, i * B, 1, cnt);
return sum;
}
void Doit() {
F(i, 2, up) {
ll tmp = i, yes = 0;
F(j, 1, P[0]) {
if (sqr(P[j]) > i) break;
ll cnt = 0;
while (tmp % P[j] == 0) tmp /= P[j], cnt ++;
yes = gcd(yes, cnt);
if (yes == 1) break;
}
if (tmp > 1 || yes == 1)
ans -= Solve(i);
}
printf("%lld\n", ans);
}
int main() {
Init();
GetPrime();
Doit();
}
【GDOI2016模拟3.16】幂(容斥 + 模型复杂转化)的更多相关文章
- bzoj 3930: [CQOI2015]选数【快速幂+容斥】
参考:https://www.cnblogs.com/iwtwiioi/p/4986316.html 注意区间长度为1e5级别. 则假设n个数不全相同,那么他们的gcd小于最大数-最小数,证明:则gc ...
- [CSP-S模拟测试]:阴阳(容斥+计数+递推)
题目传送门(内部题16) 输入格式 第一行两个整数$n$和$m$,代表网格的大小.接下来$n$行每行一个长度为$m$的字符串,每个字符若为$W$代表这个格子必须为阳,若为$B$代表必须为阴,若为$?$ ...
- CF285E Positions in Permutations(dp+容斥)
题意,给定n,k,求有多少排列是的 | p[i]-i |=1 的数量为k. Solution 直接dp会有很大的后效性. 所以我们考虑固定k个数字使得它们是合法的,所以我们设dp[i][j][0/1] ...
- [BZOJ3622]已经没有什么好害怕的了(容斥DP)
给定两个数组a[n]与b[n](数全不相等),两两配对,求“a比b大”的数对比“b比a大”的数对个数多k的配对方案数. 据说做了这题就没什么题好害怕的了,但感觉实际上这是一个套路题,只是很难想到. 首 ...
- hdu5072(鞍山regional problem C):容斥,同色三角形模型
现场过的第四多的题..当时没什么想法,回来学了下容斥,又听学长讲了一讲,终于把它过了 题目大意:给定n个数,求全部互质或者全部不互质的三元组的个数 先说一下同色三角形模型 n个点 每两个点连一条边(可 ...
- HUST 1569(Burnside定理+容斥+数位dp+矩阵快速幂)
传送门:Gift 题意:由n(n<=1e9)个珍珠构成的项链,珍珠包含幸运数字(有且仅由4或7组成),取区间[L,R]内的数字,相邻的数字不能相同,且旋转得到的相同的数列为一种,为最终能构成多少 ...
- 2017ACM暑期多校联合训练 - Team 8 1011 HDU 6143 Killer Names (容斥+排列组合,dp+整数快速幂)
题目链接 Problem Description Galen Marek, codenamed Starkiller, was a male Human apprentice of the Sith ...
- [JZOJ 5894] [NOIP2018模拟10.5] 同余方程 解题报告(容斥)
题目链接: http://172.16.0.132/senior/#contest/show/2523/0 题目: 题解:(部分内容来自https://blog.csdn.net/gmh77/arti ...
- [CSP-S模拟测试]:连连看(图论+容斥)
题目传送门(内部题74) 输入格式 输入文件$link.in$ 第一行三个整数$n,m,k$,之间用空格隔开,$n,m$表示地图行数和列数,$k$表示每个方块周围相邻的位置(至多有$4$个,至少有$2 ...
随机推荐
- 前后端分离djangorestframework—— 在线视频平台接入第三方加密防盗录视频
加密视频 在以后的开发项目中,很可能有做在线视频的,而在线视频就有个问题,因为在线播放,就很有可能视频数据被抓包,如果这个在线视频平台有付费视频的话,这样就会有人做点倒卖视频的生意了,针对这个问题,目 ...
- Oracle——DQL、DML、DDL、DCL
1.DQL:数据查询语言 基本结构:由select.from.where组成 子句组成的查询块: SELECT <字段名表> FROM <表或视图名> WHE ...
- c/c++ 多线程 等待一次性事件 std::promise用法
多线程 等待一次性事件 std::promise用法 背景:不是很明白,不知道为了解决什么业务场景,感觉std::async可以优雅的搞定一切的一次等待性事件,为什么还有个std::promise. ...
- IDF-简单题目writeup
1. 被改错的密码 原题: 从前有一个熊孩子入侵了一个网站的[数据库],找到了管理员密码,手一抖在[数据库]中修改了一下,现在的密码变成了ca9cc444e64c8116a30la00559c042b ...
- 【Teradata】安装SQL Assistant和Administrator 16.20(含查看.net版本)
1.安装介质获取: 获取的路径:connections==>Gateways==>Customer Services==>TOOLS & APPLICATIONS(点击Mor ...
- Python开发【第四篇】函数
函数的作用 函数可以让编程逻辑结构化以及模块化 无论是C.C++,Java还是Python,函数是必不可少的知识点,也是很重要的知识点,函数是完成一个功能的代码块,使用函数可以使逻辑结构变得更加清晰以 ...
- Spark1.6之后为何使用Netty通信框架替代Akka
解决方案: 一直以来,基于Akka实现的RPC通信框架是Spark引以为豪的主要特性,也是与Hadoop等分布式计算框架对比过程中一大亮点. 但是时代和技术都在演化,从Spark1.3.1版本开始,为 ...
- (四)esp8266 MDNS域名服务
(实例一)ESP8266 TFT(ST7735)彩屏-web刷图 https://www.arduino.cn/thread-42247-1-1.html (实例二) 自己当AP时建立MDNS域名 h ...
- (十四)Exploring Your Data
Sample Dataset Now that we’ve gotten a glimpse of the basics, let’s try to work on a more realistic ...
- 秒杀系统HTML倒计时设置
倒计时就是在秒杀表中设定好要秒杀的时间 再减去现在的时间 把时间差显示在页面上 由于第一次搞秒杀倒计时 辛苦搞了半天 留此微博 以留备用 <!DOCTYPE html PUBLIC &qu ...