【玲珑杯Round17】xjb总结

zcy真是垃圾，啥都不会的那种。

菜的不行。

这场手速上了三题，然后各种E被卡……

日个吗居然E不开栈，傻逼吧

有毒吧

来看题：

A.sqc给的我的神奇公式，gtmd居然能A？

#include<bits/stdc++.h>
using namespace std;
int n;double ans=;
int main(){
    scanf("%d",&n);
    for(int i=;i<n-;i++){
        int j=(n-i)/,k=(n-i)-(n-i)/;
        ans=max(ans,sin(i)+sin(j)+sin(k));
    }
    printf("%.9lf\n",ans);
}

B.小学生容斥，2^m次系列。

#include <stdio.h>
#include <iostream>
#include <algorithm>
#include <memory.h>
using namespace std;
typedef long long ll;
#define N 2
const int maxn=;
const int L=<<;
int a[maxn],T,n,m;ll f[L];int bin[L];
ll gcd(ll n,ll m){
    return m?gcd(m,n%m):n;
}
inline ll lcm(ll n,ll m){return n/gcd(n,m)*m;}
void dfs(int stp,int s,ll w){
    if (w>n) return ;
    if (stp>m) {f[s]+=n/w;return ;}
    dfs(stp+,s<<,w);
    dfs(stp+,s<<|,lcm(w,a[stp]));
}
int main(){
 
    for (int i=;i<L;i++) bin[i]=bin[i>>]+(i&);
    for (scanf("%d",&T);T--;) {
        memset(f,,sizeof f);
        scanf("%d %d",&n,&m);
        for (int i=;i<=m;i++) scanf("%d",&a[i]);
        dfs(,,);ll ans=;
        for (int j=;j<m;j++)
        for (int i=;i<(<<m);i++)
        if (((i>>j)&)==) f[i]-=f[i+(<<j)];
        for (int i=;i<(<<m);i++)
        if (bin[i]&) ans+=f[i];
        printf("%lld\n",ans);
    }
    return ;
}

剩下就一堆不会……

Ctm烦得跟什么一样，Etm卡空间卡常数，还狗日的不开栈。

D裸FFT算下。

#include<bits/stdc++.h>
#define N 1000005
typedef long long ll;
using namespace std;
const double pi=acos(-);
struct E{
    double x,y;
    E (double x=,double y=):x(x),y(y){}
}a[N];
E operator+(const E &a,const E &b){return E(a.x+b.x,a.y+b.y);}
E operator-(const E &a,const E &b){return E(a.x-b.x,a.y-b.y);}
E operator*(const E &a,const E &b){return E(a.x*b.x-a.y*b.y,a.x*b.y+a.y*b.x);}
int n,m,l,q,L,nmax;
int r[N],c[N];ll ans[N];
void FFT(E *x,int k){
    for (int i=;i<n;i++) if (i<r[i]) swap(x[i],x[r[i]]);
    for (int i=;i<n;i<<=){
        E wn(cos(pi/i),sin(pi*k/i));
        for (int j=;j<n;j+=(i<<)){
            E w(,),l,r;
            for (int p=j;p<j+i;p++,w=w*wn){
                l=x[p],r=x[p+i];
                x[p]=l+w*r;
                x[p+i]=l-w*r;
            }
        }
    }
    if (k==-) for (int i=;i<=n;i++) x[i].x/=n;
}
inline int read(){
    int f=,x=;char ch;
    do{ch=getchar();if(ch=='-')f=-;}while(ch<''||ch>'');
    do{x=x*+ch-'';ch=getchar();}while(ch>=''&&ch<='');
    return f*x;
}
int main(){
    int T=read();
    while(T--){
        n=read();q=read();m=n;nmax=;L=;
        memset(ans,,sizeof(ans));memset(a,,sizeof(a));
        for(int i=;i<=m;i++)c[i]=read(),nmax=max(n,c[i]),a[c[i]].x++;
        nmax*=;for (n=;n<=nmax;n<<=,L++);
        for(int i=;i<n;i++)r[i]=(r[i>>]>>)|((i&)<<(L-));
        FFT(a,);
        for(int i=;i<=n;i++)a[i]=a[i]*a[i];
        FFT(a,-);
        for(int i=;i<=n;i++)ans[i]=(ll)(a[i].x+0.5);
        for(int i=;i<=m;i++)ans[c[i]*]--;
        for(int i=;i<=n;i++)ans[i]>>=;
        for(int i=n;i>=;i--)ans[i]+=ans[i+];
        for(int i=;i<=q;i++){
            int k=read();printf("%lld\n",ans[k]);
        }
    }
}

好气啊，要是知道开栈我就能rank1了

巨气。