7.22训练总结

考场错误：
这次考试基础题，wa了两发，所以导致排名比较低，应该多加注意
注意环形图处理的一些细节！

HackerRank - geometric-trick

关于枚举$i\ast j=k^2$这样的形式，要注意$i,j$的次数为奇数的质因数应该相同才可以
这样我们可以给每个质数赋一个哈希值，然后把他们异或起来作为这个数的哈希值，放在一个map里面，进行配对即可

#include<bits/stdc++.h>
using namespace std;
typedef unsigned long long ull;
const int maxn=6e5+5;
const ull base=233333;
int pr[maxn],cnt,is[maxn],mn[maxn];
ull has[maxn];
void init()
{
    for(int i=2;i<maxn;i++)
    {
        if(!is[i])
        {
            pr[++cnt]=i;
            mn[i]=i;
        }
        for(int j=1;j<=cnt && 1ll*pr[j]*i<maxn;j++)
        {
            is[i*pr[j]]=1;
            mn[i*pr[j]]=pr[j];
            if(i%pr[j]==0) break;
        }
    }
}
int n;
char s[maxn];
map <ull,vector<int> > G;
int main()
{
    init();
    scanf("%d",&n);
    has[0]=1;
    for(int i=1;i<=n;i++) has[i]=has[i-1]*233333+rand();
    scanf("%s",s+1);
    ull res=0;
    long long ans=0;
    // cerr<<mn[16];
    for(int i=1;i<=n;i++)
    {
        res=0;
        for(int j=i;j>1;j/=mn[j]) res^=has[mn[j]];
        // cerr<<res<<endl;
        if(s[i]!='b' && G.count(res))
        {
            for(auto j:G[res])
            {
                int v=sqrt(1ll*i*j);
                if(s[v]=='b' && s[i]=='a' && s[j]=='c') ans++;
                if(s[v]=='b' && s[j]=='a' && s[i]=='c') ans++;
            }
        }
        G[res].push_back(i);
    }
    printf("%lld\n",ans);
    return 0;
}

HackerRank - colliding-circles

这是一个推式子的题，考场上不少人切掉了，应该注意提升这类问题

#include<bits/stdc++.h>
using namespace std;
const int maxn=1e5+5;
int n,k,r[maxn];
int main()
{
    scanf("%d%d",&n,&k);
    double sum=0,res=0;
    for(int i=1;i<=n;i++)
    {
        scanf("%d",&r[i]);
        sum+=r[i];
        res+=1.0*r[i]*r[i];
    }
    sum=sum*sum-res;
    double ans=0,now=1;
    for(int i=1;i<=k;i++)
    {
        double p=2.0/(n-i+1)/(n-i);
        ans+=p*now;
        now*=(1-p);
    }
    printf("%.10lf",acos(-1)*(ans*sum+res));
    return 0;
}

HackerRank - nominating-group-leaders

莫队裸题，为了平衡复杂度，使用分块记录出现次数为i的在j块里有几个为cnt[i][j]，从前往后，如果这个块里有值，那么就可以做到总复杂度$O(n\sqrt{n})$

#include<bits/stdc++.h>
using namespace std;
const int maxn=1e5+5;
int n,sz,v[maxn],m;
int cnt[maxn][320],Sz;
struct query
{
    int l,r,z,id;
    bool operator < (const query &oth) const{
        if(l/sz!=oth.l/sz) return l<oth.l;
        if((l/sz)&1) return r<oth.r;
        return r>oth.r;
    }
}q[maxn];
int gs[maxn],bl[maxn];
void add(int x)
{
    cnt[gs[x]][bl[x]]--;
    gs[x]++;
    cnt[gs[x]][bl[x]]++;
}
void del(int x)
{
    cnt[gs[x]][bl[x]]--;
    gs[x]--;
    cnt[gs[x]][bl[x]]++;
}
int ans[maxn];
int main()
{
    freopen("a.in","r",stdin);
    freopen("a.out","w",stdout);
    int T; scanf("%d",&T);
    while(T--)
    {
        memset(cnt,0,sizeof(cnt));
        int n; scanf("%d",&n); Sz=sqrt(n);
        // cerr<<n;
        for(int i=1;i<=n;i++) scanf("%d",&v[i]),v[i]++,gs[i]=0;
        scanf("%d",&m); sz=sqrt(m); 
        for(int i=1;i<=m;i++)
        {
            scanf("%d%d%d",&q[i].l,&q[i].r,&q[i].z),q[i].id=i;
            q[i].l++; q[i].r++;
            ans[i]=0;
        }
        for(int i=1;i<=n;i++) bl[i]=(i+Sz-1)/Sz;
        sort(q+1,q+m+1);
        int nowl=1,nowr=0;
        // cerr<<m;
        for(int i=1;i<=m;i++)
        {
            while(nowr<q[i].r) add(v[++nowr]);
            while(nowr>q[i].r) del(v[nowr--]);
            while(nowl<q[i].l) del(v[nowl++]);
            while(nowl>q[i].l) add(v[--nowl]);
            // cerr<<cnt[3][1]<<endl;
            for(int j=1;j<=bl[n];j++) if(cnt[q[i].z][j])
            {
                // cerr<<q[i].l<<" "<<q[i].r<<" "<<gs[1]<<endl;
                for(int k=(j-1)*Sz+1;k<=Sz*j;k++) if(gs[k]==q[i].z) {ans[q[i].id]=k; break;}   
                break;
            }
        }
        // cerr<<m;
        for(int i=1;i<=m;i++) printf("%d\n",ans[i]-1);
    }
    return 0;   
}