[考试反思]1111csp-s模拟测试110：三思

题目名是为了照应3天的倒计时（我才不会说是因为我考场又摸鱼了）

在OJ上得到了295的好成绩，但是本地评测没有O2掉了10分。

总体来说还可以。T1全场切，T2半场切，T3纯暴力不卡常都有95。。。

又是AK场，做的不够好。

如上所述，有点摸鱼了。但是T3说是原题。。。正解的思想的确也想不到。。。

啊，总之，还有3天，加油吧！

这一场是Day1，和明天一起记分，和前面比起来还是落后了一些，Day2加油吧

T1：最大或

半结论题，很好想。

首先对于lr相同的部分最后的答案也一定与之相同。

最高的不同位，两个数分别取1/0，然后后面的位数都取0/1，这就是最优决策。

 #include<cstdio>
 int main(){
     freopen("maxor.in","r",stdin);freopen("maxor.out","w",stdout);
     int t;long long l,r,ans;
     scanf("%d",&t);
     while(t--){ans=;
         scanf("%lld%lld",&l,&r);
         for(int i=;~i;--i)
             if((l&1ll<<i)==(r&1ll<<i))ans|=l&1ll<<i;
             else {ans+=(1ll<<i+)-;break;}
         printf("%lld\n",ans);
     }
 }

这题有代码？

T2：答题

看到数据范围：40？肯定折半搜索meet in middle

meet in middle+two pointers板子题。

暴力思路就是分别处理前后20位，然后枚举前面的$2^{20}$个数lower_bound。

复杂度为$O(2^{20} \times 20 \times log 40000000)$

考场上在这里卡了一会，想卡常，计算发现处理前18位/后22位的总复杂度更低。

预处理：$O(22 \times 2^{22})$(状压枚举+sort)

二分答案：$O(26 \times 22 \times 2^{18})$(二分答案+lower_bound)

据说starsing用这种方式AC了。

正解：可以发现lower_bound挺蠢的，如果你把前几位也排序的话，那么后几遍的指针移动就是单调的了。

一个two pointers解决问题。

 #include<cstdio>
 #include<iostream>
 #include<algorithm>
 #include<cmath>
 using namespace std;
 int n,w[],pre,suf,W[],V[],tot;long double p;
 bool chk(int s){
     int ptr=<<suf;long long totp=;
     for(int i=;i<<<pre;++i){
         while(ptr&&V[ptr-]+W[i]>s)ptr--;
         totp+=ptr;
     }
     return totp*pow(0.5,n)>=p;
 }
 int main(){
     freopen("answer.in","r",stdin);freopen("answer.out","w",stdout);
     cin>>n>>p;
     for(int i=;i<=n;++i)cin>>w[i],tot+=w[i];
     pre=n>>;suf=n-pre;
     for(int i=;i<<<pre;++i)for(int j=;j<pre;++j)if(i&<<j)W[i]+=w[j+];
     for(int i=;i<<<suf;++i)for(int j=;j<suf;++j)if(i&<<j)V[i]+=w[j+pre+];
     sort(W,W+(<<pre));sort(V,V+(<<suf));
     int l=,r=tot,ans;
     while(l<=r)if(chk(l+r>>))r=ans=l+r>>,r--;else l=(l+r>>)+;
     cout<<ans<<endl;
 }

T3：联合权值·改

 #include<cstdio>
 #include<bitset>
 #include<ctime>
 using namespace std;
 int max(int a,int b){return a>b?a:b;}
 bitset<>B[];
 int w[],fir[],l[],to[],ec,n,m,t,mx;long long tot;
 void link(int a,int b){l[++ec]=fir[a];fir[a]=ec;to[ec]=b;}
 int main(){
     freopen("link.in","r",stdin);freopen("link.out","w",stdout);
     scanf("%d%d%d",&n,&m,&t);
     for(int i=,a,b;i<=m;++i)scanf("%d%d",&a,&b),B[a][b]=B[b][a]=,link(a,b),link(b,a);
     for(int i=;i<=n;++i)scanf("%d",&w[i]),B[i][i]=;
     for(int i=;i<=n;++i){
         if(t==&&clock()>)return printf("%d\n0\n",mx),;
         for(int j=fir[i];j;j=l[j])for(int k=fir[to[j]];k;k=l[k])
             if(!B[i][to[k]])tot+=w[i]*w[to[k]],mx=max(mx,w[i]*w[to[k]]);
     }
     printf("%d\n%lld\n",t==?:mx,t==?:tot);
 }

95分暴力

好题。难题。（呃，乱搞不算）

如果写正解的话数据范围开到100000，权值出到10000，内存卡到8MB都是可做的。

但是出题人比较仁慈数据又小又水，让很多人多拿了不少分数～（就像我的95分暴力）

但是暴力归暴力，得分归得分。还是要承认：正解的确很牛逼%%%skyh

正解可以做到在任何题目里以$O(m^{1.5})$的复杂度求出图里的所有三元环。（无自环/重边），思路值得积累。

以前都是bitset$O(\frac{mn}{32})$暴干的，没什么水平。

首先一个结论是三元环个数不会超过$O(m^{1.5})$。因为在完全图里会达到上届。

对于所有点，我们将其按照度数排序，度数相同的要内部严格有序。得到每个点的排名。

然后重新建边，边只从排名小的连向大的。

这样的话我们枚举每一个点i，找出所有出边j，再找出所有到达点的出边k，检查i，k两点之间是否有边即可。

这样的话，对于每个度数大于根号的点j，它们的总出边小于等于m，而每个点的入边小于根号，所以总数小于$m \sqrt{m}$

如果一个点的度数小于根号，那么它们的总入边小于等于m，而每个点的出边小于根号，所以总数小于$m \sqrt{m}$

这样，计算总值的时候，枚举每一个点计算其邻点的和的平方减去平方的和，再取出三元环贡献即可。

计算最大值时，从大到小枚举邻点（二重循环），遇到非三元环就跳出。因为三元环不超过$m\sqrt{m}$所以总复杂度也是$m \sqrt{m}$\

代码也好写。是道好题。

 #include<bits/stdc++.h>
 using namespace std;
 bitset<>B[];
 struct Q{int p,w;friend bool operator<(Q x,Q y){return x.w>y.w;}}q[];
 vector<Q>v[];
 struct P{int deg,ord;friend bool operator<(P x,P y){return x.deg>y.deg;}}p[];
 int w[],fir[],l[],to[],ec,n,m,t,mx,FIR[],L[];
 int deg[],TO[],rk[],EC;long long tot;
 void link(int a,int b){l[++ec]=fir[a];fir[a]=ec;to[ec]=b;deg[b]++;}
 void LINK(int a,int b){L[++EC]=FIR[a];FIR[a]=EC;TO[EC]=b;}
 int main(){
     freopen("link.in","r",stdin);freopen("link.out","w",stdout);
     scanf("%d%d%d",&n,&m,&t);
     for(int i=,a,b;i<=m;++i)scanf("%d%d",&a,&b),B[a][b]=B[b][a]=,link(a,b),link(b,a);
     for(int i=;i<=n;++i)scanf("%d",&w[i]),B[i][i]=,p[i]=(P){deg[i],i},q[i]=(Q){i,w[i]};
     sort(p+,p++n);
     for(int i=;i<=n;++i)rk[p[i].ord]=i;
     for(int i=;i<=n;++i)for(int j=fir[i];j;j=l[j])if(rk[to[j]]<rk[i])LINK(i,to[j]);
     for(int i=;i<=n;++i){
         long long sum=;
         for(int j=fir[i];j;j=l[j])tot-=w[to[j]]*w[to[j]],sum+=w[to[j]];
         tot+=sum*sum;
     }
     for(int i=;i<=n;++i)for(int j=FIR[i];j;j=L[j])for(int k=FIR[TO[j]];k;k=L[k])
         if(B[i][TO[k]])tot-=w[i]*w[TO[j]]+w[TO[j]]*w[TO[k]]+w[i]*w[TO[k]]<<;
     for(int i=;i<=n;++i)for(int j=fir[i];j;j=l[j])v[i].push_back(q[to[j]]);
     for(int i=;i<=n;++i)sort(v[i].begin(),v[i].end());
     for(int i=;i<=n;++i)for(int j=;j<v[i].size();++j)for(int k=j+;k<v[i].size();++k)
         if(!B[v[i][j].p][v[i][k].p]){mx=max(mx,v[i][j].w*v[i][k].w);break;}
         else if(v[i][j].w*v[i][k].w<=mx)break;
     printf("%d\n%lld\n",t==?:mx,t==?:tot);
 }