UVA 12716 GCDXOR 数论
题目链接: https://vjudge.net/problem/UVA-12716
题目描述: 给你一个N, 让你求1~N内所有的GCD(A,B) == A XOR B的个数, 其中1 <= B <= A <= N , N <= 3e7
解题思路: gcd(a, b) = a ^ b = c, 异或的性质是a ^ b = c 则 a ^ c = b 所以我们就枚举a, c得出b , 此时如果gcd(a, b) == c则++ans, 枚举a是O(n), 枚举因数模仿素数筛法枚举a, c是O(nlogn), 判断gcd是logn所以是O(nlognlogn), 这里N是3e7照理说是不会T了, 但是确实是T了, 可能是这个方法卡常, 这里我们还有O(nlogn)的做法, 这样是不会T的, 由于异或是不进位的加法, 所以a - b <= a ^ b 相等的时候是a == b , 同时a - b >= c 由于a ^ b == c 所以说a-b>=c && a-b<=所以a-b==c , 所以这里我们就能O(1)得出b, 再O(1)判断a^b是否==c这样复杂度就能降下来了, 题目就可以A了
代码:
#include <iostream> #include <cstdio> #include <string> #include <vector> #include <cstring> #include <iterator> #include <cmath> #include <algorithm> #include <stack> #include <deque> #include <map> #define lson l, m, rt<<1 #define rson m+1, r, rt<<1|1 #define mem0(a) memset(a,0,sizeof(a)) #define meminf(a) memset(a,0x3f,sizeof(a)) #define fi(n) for(i=0;i<n;i++) #define fj(m) for(j=0;j<m;j++) #define sca(x) scanf("%d",&x) #define scalld(x) scanf("%I64d",&x) #define print(x) printf("%d\n",x) #define printlld(x) printf("%I64d\n",x) #define d printf("=======\n") typedef long long ll; using namespace std; const int maxn = 3e7; ll ans[maxn+10]; ll gcd( ll a, ll b ) { return( b == 0 ? a : gcd( b, a%b ) ); } void build() { mem0(ans); int n = maxn >> 1; for( int c = 1; c <= n; c++ ) { for( int a = 2*c; a <= maxn; a+=c ) { ll b = a ^ c; // TLE!!!!!!!!!!!!!!!! if( a >= b && gcd(a, b) == c ) ++ans[a]; // ll b = a-c; // if( (a ^ b) == c ) ++ans[a]; } } // for( int i = 1; i <= 10; i++ ) { // cout << ans[i] << " "; // } // cout << endl; for( int i = 2; i <= maxn; i++ ) { ans[i] += ans[i-1]; } } int main() { int t; // cout << gcd(9,15) << endl; build(); sca(t); int cases = 1; while( t-- ) { ll res = 0; int n; sca(n); res = ans[n]; printf( "Case %d: %lld\n", cases++, res ); } return 0; }
思考 : 数论真的是好神奇啊, 对数字的敏感性要求是真的高, 而且这个东西得自己动脑筋去想才能够解题, 同时还要掌握更多的定理来支持自己, 一会儿hihocoder啦, 激动ing......
posted on 2017-08-20 15:19 FriskyPuppy 阅读(153) 评论(0) 编辑 收藏 举报