CF 241B
来源:互联网 发布:禁止软件安装的软件 编辑:程序博客网 时间:2024/06/05 05:02
题目很好可惜我不会做。。。
看了别人代码半天才懂。。。
看完感觉还是没学到啥。。。
于是写点东西记一下吧。。。
一堆数异或的问题,一个常用的思路是枚举每一位然后瞎搞
这里的做法是,先确定一个前缀,然后想办法求所有数对能得到多少个这个前缀
然后就可以根据求得的个数和剩下还需要求的个数的关系来计算答案跟确定下一个前缀
于是现在问题变成了有一个前缀,怎么确定这些数能搞出多少这个前缀
于是这里可以枚举所有数,对每个数求出,它异或这个前缀得出来的前缀在其他数中有多少。。。这样每对会被计算两遍,答案除个2就行了,值得一提的是不取模的情况下答案也不会爆longlong,那就可以求出最终答案再取模,当然也可以一边算一边模,就是要多写个快速幂求逆
那么显然这个子问题可以用trie解决,复杂度O(nlogn)
至此问题完美解决!总复杂度O(n(logn)^2)
(如果用静态分配内存记得算好要用多少空间。。。要卡得刚刚好。。。
代码
#include<iostream>#include<cstdio>#include<cstdlib>#include<cstring>#include<cmath>#include<algorithm>#include<string>#include<iomanip>#include<vector>#include<set>#include<map>#include<queue>using namespace std;typedef long long LL;typedef unsigned long long ULL;#define rep(i,k,n) for(int i=(k);i<=(n);i++)#define rep0(i,n) for(int i=0;i<(n);i++)#define red(i,k,n) for(int i=(k);i>=(n);i--)#define sqr(x) ((x)*(x))#define clr(x,y) memset((x),(y),sizeof(x))#define pb push_back#define mod 1000000007const int maxn=50003;const int maxnode=maxn*17;const int maxl=30;int n,a[maxn];LL m;struct trie{ int next[maxnode][2],sz[maxnode],c[maxnode][maxl+1][2]; int root,L; int newnode() { for(int i=0;i<2;i++)next[L][i]=-1; sz[L]=0; clr(c[L],0); return L++; } void init() { L=0; root=newnode(); } void insert(char str[]) { int len=strlen(str); int now=root; for(int i=0;i<len;i++) { sz[now]++; char s=str[i]-'0'; if(next[now][s]==-1) next[now][s]=newnode(); now=next[now][s]; } sz[now]++; } void dfs(int x,int k) { rep(i,0,1) { if(next[x][i]==-1)continue; int y=next[x][i]; dfs(y,k-1); rep(j,0,maxl) { c[x][j][0]+=c[y][j][0]; c[x][j][1]+=c[y][j][1]; } c[x][k][i]+=sz[y]; } } int f(int x,int k,int dep,int low) { if(dep==low)return x; int nxt=(k>>dep)&1; if(next[x][nxt]==-1)return -1; return f(next[x][nxt],k,dep-1,low); } LL g1(int pre,int k,LL &tmp) { LL ret=0;tmp=0; rep(i,1,n) { int p=f(root,a[i]^pre,maxl,k-1); if(p==-1)continue; ret+=sz[p]; tmp+=(LL)pre*sz[p]; rep(j,0,k-1)tmp+=(LL)c[p][j][((a[i]>>j)&1) ^1]<<j; } return ret; } LL gao(LL k) { LL ret=0,cnt=0,tmp=0; dfs(root,maxl); int pre=0; red(i,maxl,0) { pre|=1<<i; cnt=g1(pre,i,tmp); if(cnt<=k) { k-=cnt; ret+=tmp; pre^=1<<i; } } cnt=g1(pre,0,tmp); ret+=tmp/cnt*k; return ret; }}tr;inline void convert(int x,char s[]){ for(int i=30;i>=0;i--) { s[i]='0'+x%2; x>>=1; } if(x)cout<<"re"<<endl; s[31]=0;}int main(){ scanf("%d%I64d",&n,&m); char s[40]; tr.init(); rep(i,1,n) { scanf("%d",&a[i]); convert(a[i],s); tr.insert(s); } cout<<(tr.gao(m<<1)>>1)%mod<<endl; return 0;}
0 0
- CF 241B
- cf-B
- cf B
- CF 1B Spreadsheet
- CF 3B Lorry
- CF 176B
- CF 91B
- cf 131 DIV2 B
- CF 113B || CF196D
- CF 126B
- CF-236B
- CF 254B(日期)
- cf 79B
- CF - 158B - Taxi
- CF 296B
- CF 173(div2) B
- cf:B-QR code
- cf 1B. Spreadsheets
- hdu4614 Vases and Flowers(线段树+二分)
- JS canvas 表情
- 面向对象的一些心得
- eclipse中server中无法加载tomcat了
- 输出乘法口诀表
- CF 241B
- SDWebImage与UIScrollView或UICollectionView一起使用多次收到“Received memory warning.”,最终程序闪退的原因分析和解决办法。
- 更改host的重定向
- poj 1698 Just a Hook -- (线段树,区间值修改,区间求和)
- Http 请求
- 购买腾讯云服务器后的配置问题分析
- LeetCode---Pascal Triangle、Pascal Triangle II解题分析
- C Tricks(二)—— 按块访问二维数组
- 傅里叶变换