纪中9日T4 2298. 异或
2298. 异或
(File IO): input:gcdxor.in output:gcdxor.out
时间限制: 1000 ms 空间限制: 262144 KB 具体限制
题目描述
SarvaTathagata是个神仙,一天他在研究数论时,书上有这么一个问题:求不超过n两两的数的gcd。
SarvaTathagata这么神仙的人当然觉得这个是sb题啦。学习之余,他还发现gcd的某一个特别好的性质:如果有两个数i,j满足gcd(i,j)=i^j(这里的^为c++中的异或)的话,那么这两个数组成的数对(i,j)就是一个nb的数对(这里认为(i,j)和(j,i)为相同的,并不需要计算2次)。
当然,SarvaTathagata并不会只满足于判断一个数对是否nb,他还想知道满足两个数都是不超过n并且nb的数对有多少个。
由于SarvaTathagata实在是太神仙了,他认为这种题实在是太简单了。于是他找到了你,看看你是否能解决这个问题。
输入
共一行一个整数n,含义如题所述。
输出
一行一个整数,表示nb的数对的个数。
样例输入
样例输入1
12
样例输入2
123456
样例输出
样例输出1
8
样例输出2
214394
数据范围限制
提示
样例1中共有八对,分别是:
{1,3},{1,5},{1,7},{1,9},{2,6},{1,11},{2,10},{4,12}。
提示有误,特此隐藏
Solution
这是一道数论题,涉及gcd,xor(^),二进制减法。
以上三者中,若你对任何一样过敏,请谨慎食用。
Way one(40分)
用个程序找规律
//gcdxor table
#include<bits/stdc++.h>
using namespace std;
int gcd(int ta,int tb)
{
if(tb==) return ta;
if(ta%!=&&tb%!=) return gcd(tb,ta%tb);
if(ta%==&&tb%!=) return gcd(ta/,tb);
if(ta%!=&&tb%==) return gcd(ta,tb/);
return *gcd(ta/,tb/);
} int ans[];
int n;
bool vis[],memory[][];//memory[i][j]
int search(int num)
{
if(vis[num]) return ans[num];
if(num==)
return ; vis[num]=;
int now=;
for(int i=;i<num;i++)
{
if(gcd(num,i)==(num xor i))
now++;
// if(gcd(i,num)==(i xor num))
// now++; }
// if(gcd(num,num)==(num xor num))
// now++;
ans[num]=ans[num-]+now;
return ans[num];
}
int main()
{
freopen("table4.txt","w",stdout);
cin>>n;
// int s=clock();
int k;
/*
for(k=0;k<=n;k++)
{
if(k==0) {cout<<0<<endl;continue;}
for(int i=1;i<=k;i++)
{
for(int j=1;j<=k;j++)
{
if(gcd(i,j)==(i xor j)){
memory[i][j]=true;
ans[k]++; }
// cout<<memory[i][j]<<" ";
}
// cout<<endl;
} cout<<k<<" "<<ans[k]/2<<endl;
}
*/
// /*
for(int k=;k<=n;k++)
{
if(k==) {cout<<<<endl;continue;}
cout<<k<<" "<<search(k)<<","<<endl;
}
// */
// int e=clock();
// cout<<e-s; return ;
}
可以无视我
找不到……
Code(40分)
#pragma GCC optimize(2)
#include<bits/stdc++.h>
#define IL inline
using namespace std;
short int diff[]{,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,};
int n,ans;
IL int gcd(int ta,int tb)
{
if(tb==) return ta;
if(ta%!=&&tb%!=) return gcd(tb,ta%tb);
if(ta%==&&tb%!=) return gcd(ta/,tb);
if(ta%!=&&tb%==) return gcd(ta,tb/);
return *gcd(ta/,tb/);
}
IL int search(int num)
{
if(num==) {
int t=;
for(int i=;i<=;i++) t+=diff[i];
return t;
}
int prev;
if(num>)
prev=search(num-);
int now=;
for(int i=;i<num;i++) if(gcd(num,i)==(num xor i)) now++;
if(gcd(num,num)==(num xor num)) now++;
return prev+now;
}
int main()
{
freopen("gcdxor.in","r",stdin);
freopen("gcdxor.out","w",stdout);
cin>>n;
for(int i=;i<=min(n,);i++)
ans+=diff[i];
if(n<=)
cout<<ans<<endl;
else
cout<<search(n)<<endl;
return ;
}
40分
其实,要是我愿意我可以把表全打出来,但是这个OJ有代码长度限制(5kb)。有毒……
Way two
从大神视角理解这道题:
看到gcd就去想数论嘛
求证:若gcd(a,b)=a xor b,则gcd(a,b)=a-b (a≥b)
转到完整证明
证明:
由更相减损法,可知gcd(a,b)=gcd(b,a-b)
这个的原理与辗转相除法类似,只是把求模换成了减法
若a-b=0,则gcd(a,b)=gcd(b,a-b)=a=b
否则a-b>0,此时gcd(a,b)=gcd(b,a-b)>a-b
所以,gcd(a,b)≥a-b
再分析xor运算和二进制减法
xor:对与每一位,若是相同则为0,若是不同则得1
1^1=0 1^0=1 0^1=1 0^0=0
减法:对于每一位,若是相同则得0,若是1 0则得1,若是0 1则得1并使上一位-1
1-1=0 1-0=1 0-1=-1 0-0=0
那么,如果两数a,b中出现了0->1的情况,则减法会使那一位上出现退位,异或(xor)则不会
所以,a xor b≥a-b
又因为(前面证明的)gcd(a,b)≥a-b
所以gcd(a,b)≥a-b≥a xor b
所以当gcd(a,b)=a xor b时,gcd(a,b)=a-b=a xor b
所以若gcd(a,b)=a xor b,则gcd(a,b)=a-b (a≥b)
证毕
回到此题
现在我已经证出来
若gcd(a,b)=a xor b,则gcd(a,b)=a-b (a≥b)
一看题目,就是当gcd(a,b)=a xor b时
那么我们就可以直接用结论了
设c=a-b,则b=a-c
若数对(a,b)符合条件
则c=a^b
故我只要枚举a和b即可!
Code(90分)
#pragma GCC optimize(2)//这个程序不开O2会超时(90分)
#include<bits/stdc++.h>
using namespace std;
int n,cnt[],sum[];
int main()
{
// freopen("gcdxor.in","r",stdin);
// freopen("gcdxor.out","w",stdout);
cin>>n;
memset(cnt, , sizeof(cnt));
for(int c = ; c <= n; c++)
for(int a = c*; a <= n; a += c) //因为a>=b,所以需要从2*c开始枚举
{
int b = a - c;
if(c == (a ^ b)) cnt[a]++;//统计每个a对应的b的数量
}//^的优先级低于==,所以要打上括号
sum[] = ;
for(int i = ; i <= n; i++)
sum[i] = sum[i-] + cnt[i];
cout<<sum[n];
return ;
}
Code(100分)
#include<bits/stdc++.h>
using namespace std;
int n,sum;
int main()
{
// freopen("gcdxor.in","r",stdin);
// freopen("gcdxor.out","w",stdout);
scanf("%d",&n);
// int s=clock();
for(int c=;c<=n;c++)//c=a-b
for(int a=c*;a<=n;a+=c) //因为a>=b,所以需要从2*c开始枚举
if(c==(a ^ (a-c))) sum++;
//统计每个a对应的b的数量
//^的优先级低于==,所以要打上括号
printf("%d",sum);
// int e=clock();
// cout<<endl<<e-s;
return ;
}
TLE?
细心的你应该发现这个:
不要怕TLE啦
如果你会算复杂度,会一点微积分,会一点金融学,你就要知道:
所以整个程序的时间复杂度为O(nlog(n))
Attention
一、^的优先级低于==,所以要打上括号
二、因为a>=b,所以需要从2*c开始枚举,使得a的最小值为2*c,b=a-c,b的初值值为c且越来越小
三、用前缀和的思想,递归sum[i] = sum[i-1] + cnt[i],即n的答案是在n-1的基础上加上这次枚举的答案
END
纪中9日T4 2298. 异或的更多相关文章
- 纪中21日T3 2118. 【2016-12-30普及组模拟】最大公约数
纪中21日T3 2118. 最大公约数 (File IO): input:gcd.in output:gcd.out 时间限制: 1000 ms 空间限制: 262144 KB 具体限制 Goto ...
- 纪中17日T1 2321. 方程
纪中17日T1 2321. 方程 (File IO): input:cti.in output:cti.out 时间限制: 1000 ms 空间限制: 262144 KB 具体限制 Goto ...
- 纪中10日T1 2313. 动态仙人掌
纪中10日 2313. 动态仙人掌 (File IO): input:dinosaur.in output:dinosaur.out 时间限制: 1500 ms 空间限制: 524288 KB 具 ...
- 洛谷P1880 [NOI1995]石子合并 纪中21日c组T4 2119. 【2016-12-30普及组模拟】环状石子归并
洛谷P1880 石子合并 纪中2119. 环状石子归并 洛谷传送门 题目描述1 在一个圆形操场的四周摆放N堆石子,现要将石子有次序地合并成一堆.规定每次只能选相邻的2堆合并成新的一堆,并将新的一堆的石 ...
- 纪中23日c组T2 2159. 【2017.7.11普及】max 洛谷P1249 最大乘积
纪中2159. max 洛谷P1249 最大乘积 说明:这两题基本完全相同,故放在一起写题解 纪中2159. max (File IO): input:max.in output:max.out 时间 ...
- 纪中23日c组T3 2161. 【2017.7.11普及】围攻 斐波那契数列
2161. 围攻 (File IO): input:siege.in output:siege.out 时间限制: 1000 ms 空间限制: 262144 KB 具体限制 Goto Prob ...
- 纪中18日c组模拟赛
T2 GMOJ2127. 电子表格 (File IO): input:excel.in output:excel.out 时间限制: 1000 ms 空间限制: 262144 KB 具体限制 ...
- 纪中21日c组T1 1575. 二叉树
1575. 二叉树 (File IO): input:tree.in output:tree.out 时间限制: 1000 ms 空间限制: 262144 KB 具体限制 Goto Probl ...
- 纪中21日c组模拟赛
AWSL AWSL AWSL AWSL AWSL AWSL AWSL AWSL AWSL AWSL AWSL AWSL AWSL AWSL AWSL AWSL 题解传送 T1 ...
随机推荐
- C#系列之占位符的使用方法(二)
今天,我将简单记录下占位符的使用方法 首先,我们来看不使用占位符的方法来代码输出 int number = 10; int number_1 = 20; int number_2 = 30; Cons ...
- phpstrom激活码
今天PHPstorm又到期了,从网上找到一个激活码的网址,很好用,说是会时时更新的,所以特意记录一下 获取地址:https://www.php.cn/tool/phpstorm/408348.html
- VS Code 1.42 发布!2020 年首个大更新
近日(北京时间 2020 年 2 月 7 日),微软发布了 Visual Studio Code 1.42 版本,这也是 2020 年 VS Code 首次大更新.让我们来看看有哪些主要的更新. 支持 ...
- 应用层vc实现三种文件监视方法
http://hi.baidu.com/sadusaga/item/daa0d4b764c6dd76254b09cc http://bbs.csdn.net/topics/280032788 http ...
- GORM入门指南
gorm是一个使用Go语言编写的ORM框架.它文档齐全,对开发者友好,支持主流数据库. gorm介绍 Github GORM 中文官方网站内含十分齐全的中文文档,有了它你甚至不需要再继续向下阅读本文. ...
- CCF_201604-3_路径解析
没有用stack来写,直接用了字符串和指针,过程有点复杂的.首先对读入的每一个路径,判断它是绝对路径或者相对路径,然后确定起始的路径,然后继续一位位的判断,"//","/ ...
- MBMD(MobileNet-based tracking by detection algorithm)作者答疑
If you fail to install and run this tracker, please email me (zhangyunhua@mail.dlut.edu.cn) Introduc ...
- java9String类简单了解
public class jh_01_String类简单了解 { public static void main(String[] args) { /* * 函数:完成特定功能的代码块. * next ...
- RTEMS进程同步机制
互斥量 好像没有互斥量,信号量接收那儿有个图,互斥量似乎术语一类特殊的信号量. 信号量 12. Semaphore Manager 12.1. Introduction The semaphore m ...
- Cacti被监控机器 配置 snmp协议
SNMP(Simple Network Management Protocol,简单网络管理协议)的前身是简单网关监控协议(SGMP),用来对通信线路进行管理. snmpd.conf的 ...