The 2019 ICPC China Nanchang National Invitational and International Silk-Road Programming Contest B、H

比赛链接https://www.jisuanke.com/contest/3098?view=challenges

B题 拉格朗日插值

题意  T组输入。一个n次多项式 f(x) ，每项的系数不知道，只知道f(0),f(1)..f(n) 的值,m个询问，L,R。计算$\sum_{i=L}^{R}f(i)\quad mod(9999991)$

$(1\leq T\leq 5) $

$(1\leq n\leq 1000) $

$(1\leq m\leq 2000) $

$(1\leq L\leq R \leq 9999990)$

解析 遇到这题我是崩溃的，听大家说是拉格朗日插值，找到了一个快速拉格朗日的板子，贴上去就过了。。。

#include<bits/stdc++.h>
using namespace std;
#define maxn 1010
typedef long long LL;
const LL mod = ;

LL powmod(LL aa, LL x) {
    LL res = ;
    for(; x > ; x >>= ) {
        if(x & )res = (res * aa) % mod;
        aa = (aa * aa) % mod;
    }
    return res;
}

struct lagrange {
#define ll long long
#define rep(i,a,n) for (int i=a;i<n;i++)
#define per(i,a,n) for (int i=n-1;i>=a;i--)
#define D 2010 //D比MAXN大100就行
    ll a[D], f[D], g[D], p[D], p1[D], p2[D], b[D], h[D][], C[D];
    void init(int M) {//初始化：参数填MAXN + 20
        f[] = f[] = g[] = g[] = ;
        rep(i, , M + ) f[i] = f[i - ] * i % mod;
        g[M + ] = powmod(f[M + ], mod - );
        per(i, , M + ) g[i] = g[i + ] * (i + ) % mod;
    }
    /*给定一组样本数据a[]，规模为0-d，计算出第n项*/
    ll calcn(int d, ll *a, ll n) {
        if (n <= d) return a[n];
        p1[] = p2[] = ;
        rep(i, , d + ) {
            ll t = (n - i + mod) % mod;
            p1[i + ] = p1[i] * t % mod;
        }
        rep(i, , d + ) {
            ll t = (n - d + i + mod) % mod;
            p2[i + ] = p2[i] * t % mod;
        }
        ll ans = ;
        rep(i, , d + ) {
            ll t = g[i] * g[d - i] % mod * p1[i] % mod * p2[d - i] % mod * a[i] % mod;
            if ((d - i) & ) ans = (ans - t + mod) % mod;
            else ans = (ans + t) % mod;
        }
        return ans;
    }
    /*
    给定一组观测点(0, a[0]), (1, a[1]), ...，(m, a[m])，、
    样本点的个数为a(x)的最高次+1。
    求在该函数模型下，a[0]+a[1]+...+a[n]的和。
    */
    ll ta[D];
    ll polysum(ll m, ll *a, ll n) { // 给定a[0].. a[m]，求\sum_{i=0}^{n}a[i]
        memcpy(ta, a, sizeof(a[]) * (m + ));
        ta[m + ] = calcn(m, ta, m + );
        rep(i, , m + )ta[i] = (ta[i - ] + ta[i]) % mod;
        return calcn(m + , ta, n);
    }
};
int main()
{
    int t;
    scanf("%d",&t);
    while(t--){
        ll a[maxn],n,m;
        scanf("%lld%lld",&n,&m);
        for(int i=;i<=n;i++){
            scanf("%lld",&a[i]);
        }
        lagrange ri;
        ri.init(maxn+);
        while(m--){
            ll l,r;
            scanf("%lld%lld",&l,&r);
            printf("%lld\n",(ri.polysum(n,a,r)-ri.polysum(n,a,l-)+mod)%mod);
        }
    }
}

H题  FWT+线段树

题意   一个n代表A，B数组的长度，A,B两个数组中的数两两或(二进制运算)一下 ，得到一个不去重C数组(显然C的长度为n*n)。接下来一个m代表操作次数，每次输入两个数L,R

如果L等于0，表示询问C数组中第R个数，否则表示C数组中第L个数到第R个数 开根号。

解析 用FWT求出来 or 之后 每个数的个数，然后建立权值线段树，或者前缀和+二分，都可以log时间复杂度求出第k个数是几。只需要知道这个数开了几次根号，L，R会很大，

但是数量只有那么多，离散化一下就可以解决了,1e5开根号大于5次就是1了，小于5次暴力开就好了。

#include <bits/stdc++.h>
using namespace std;
#define rep(i,a,n) for (int i=a;i<n;i++)
#define per(i,a,n) for (int i=n-1;i>=a;i--)
#define pb push_back
#define mp make_pair
#define all(x) (x).begin(),(x).end()
#define fi first
#define se second
#define SZ(x) ((int)(x).size())
typedef long long ll;
typedef pair<int,int> Pii;
const int maxn = 3e5+;

ll a[maxn],b[maxn],c[maxn];
void FWT_or(ll *a,int N,int opt)
{
    for(int i=;i<=N;i<<=)
        for(int p=i<<,j=;j<=N;j+=p)
            for(int k=;k<i;++k)
                if(opt==)a[i+j+k]=a[j+k]+a[i+j+k];
                else a[i+j+k]=a[i+j+k]-a[j+k];
}
ll C[maxn];
int lowbit(int x)
{
    return x&(-x);
}
ll bitgetsum(int x)
{
    ll ans=;
    for(int i=x;i>;i-=lowbit(i))
        ans+=C[i];
    return ans;
}
void bitupdate(int x,int z)
{
    for(int i=x;i<=3e5;i+=lowbit(i))
        C[i]+=z;
}
struct ndoe
{
    ll l,r;
}q[maxn];
vector<ll> v;
int getid(ll x){
    return lower_bound(all(v),x)-v.begin()+;
}
ll sum[maxn*];
void pushUp(int rt)
{
    sum[rt]=sum[rt<<]+sum[rt<<|];
}
void build(int l,int r,int rt)
{
    if(l==r){
        sum[rt]=c[l];
        return;
    }
    int m=(l+r)>>;
    build(l,m,rt<<);
    build(m+,r,rt<<|);
    pushUp(rt);
}
int query(ll val,int l,int r,int rt)
{
    if(l==r){
        return l;
    }
    int mid=(l+r)>>;
    if(sum[rt<<]>=val)
        return query(val,l,mid,rt<<);
    else
        return query(val-sum[rt<<],mid+,r,rt<<|);
}
int main()
{
    int n,x,m;
    scanf("%d",&n);
    for(int i=;i<n;i++){
        scanf("%d",&x);a[x]++;
    }
    for(int i=;i<n;i++){
        scanf("%d",&x);b[x]++;
    }
    n=2e5;
    FWT_or(a,n,);FWT_or(b,n,);
    for(int i=;i<=n;i++) c[i]=1ll*a[i]*b[i];
    FWT_or(c,n,-);build(,n,);
    scanf("%d",&m);
    for(int i=;i<m;i++){
        scanf("%lld%lld",&q[i].l,&q[i].r);
        v.pb(q[i].r);
        if(q[i].l!=)
            v.pb(q[i].l);
    }
    sort(all(v));
    v.erase(unique(all(v)),v.end());
    for(int i=;i<m;i++){
        if(q[i].l!=){
            bitupdate(getid(q[i].l),);
            bitupdate(getid(q[i].r)+,-);
        }
        else{
            int ans,times = bitgetsum(getid(q[i].r));
            if(times>=){
                ans=;
            }
            else{
                ans = query(q[i].r,,n,);
                while(times--){
                    ans=floor(sqrt(ans));
                }
            }
            printf("%d\n",ans);
        }
    }
    return ;
}