2019牛客多校第一场

A.Equivalent Prefixes

二分答案，check的时候每次找到区间最小值，如果位置不一样则返回false，否则两边继续递归

#include <bits/stdc++.h>
#define fo(i,a,b) for(int i=a;i<=b;i++)
using namespace std;
const int maxn=1e5+5;
int posa[maxn],posb[maxn],a[maxn],b[maxn];
int n,l,r,ans;
struct RMQ
{
    int n,d[maxn][20],bit[maxn];
    void init(int *a)
    {
        for (int i=1;i<=n;i++)
            d[i][0]=a[i];
        for (int j=1;(1<<j)<=n;j++)
            for (int i=1;(i+(1<<j))<=n+1;i++)
                d[i][j]=min(d[i][j-1],d[i+(1<<(j-1))][j-1]);
        bit[1]=0;
        for (int i=2;i<=n;i++)
            bit[i]=bit[i>>1]+1;
    }
    int query(int l,int r)
    {
        int k=bit[r-l+1];
        return min(d[l][k],d[r-(1<<k)+1][k]);
    }
};
RMQ u,v;
bool check(int l,int r)
{
    int mina,minb;
    if (l >= r) return true;
    mina = posa[u.query(l,r)];
    minb = posb[v.query(l,r)];
    //cout<<u.query(l,r)<<' '<<v.query(l,r)<<endl;
    if (mina != minb) return false;
    return check(l,mina-1) && check(mina+1,r);
}
int main()
{
    while (~scanf("%d",&n))
    {
        fo(i,1,n) scanf("%d",&a[i]);
        fo(i,1,n) scanf("%d",&b[i]);
        fo(i,1,n) posa[a[i]] = i,posb[b[i]] = i;
        u.n=n,v.n=n;
        u.init(a); v.init(b);
         
        l = 1; r = n;
        while (l <= r)
        {
            int mid = (l + r) >> 1;
            if (check(1,mid)) {ans = mid; l = mid + 1;} else r = mid - 1;
        }
        cout<<ans<<endl;
    }
}

B.Integration

把乘法拆成若干个加法即可

C.Euclidean Distance

一开始直接用拉格朗日发现答案可能小于零。。。
我的想法是，这个距离本质上是$N$维空间上两点间的距离，$\sum p_i=1$实际上是一个约束了一个$N-1$维的空间，而$p_i>0$实际上是这个空间的边界。
当$p_i<0$的时候，他的几何意义是这个最近点落在了边界外面，那么这时候显然$p_i=0$，即在边界上寻找答案，问题就转化为一个$N-1$维的子问题了。
而事实上题解给出的方法是$p_i$的取值范围都是可以独立计算的。

#include <bits/stdc++.h>
#define fo(i,a,b) for(i=a;i<=b;i++)
using namespace std;
typedef long long ll;
ll i,n,m,sum,nn,mark;
ll a[10005];
int main()
{
    while (~scanf("%lld%lld",&n,&m))
    {
        fo(i,1,n) cin>>a[i];
        sort(a+1,a+n+1);
        sum = 0; fo(i,1,n) sum += a[i];
        nn = n; mark = 0;
        fo(i,1,n)
        {
            if (nn * a[i] + m - sum < 0)
            {
                nn--; sum -= a[i];
                mark = i;
                continue;
            }
            break;
        }
        n = nn;
        sum=sum-m;
        ll p=sum*sum;
        ll q=m*m*n;
        fo(i,1,mark) p = p + a[i]*a[i]*n;
        if (p==0) printf("0\n");
        else
        {
            ll g=__gcd(p,q);
            p/=g;
            q/=g;
            if (q==1) printf("%lld\n",p);
            else printf("%lld/%lld\n",p,q);
        }
         
    }
    return 0;
}

E.ABBA

本质上就是一个带上下界的括号匹配（貌似并不是）
题解：贪心，$A$肯定先用$AB$的$A$，$B$肯定先用$BA$的$B$
在去掉匹配的$AB/BA$后，任何状态都不能出现大于$n$个$A$或者大于$m$个$B$

#include <bits/stdc++.h>
using namespace std;
const int mod=1e9+7;
int dp[4010][2005];
int main()
{
    int a,b;
    while (~scanf("%d%d",&a,&b))
    {
        for (int i=0;i<=2*(a+b);i++)
            for (int j=0;j<=a+b;j++)
                dp[i][j]=0;
        dp[0][a]=1;
        for (int i=0;i<2*(a+b);i++)
            for (int j=0;j<=a+b;j++)
            {
                if (j>0 && j+1<=a+b) dp[i+1][j]=(dp[i][j-1]+dp[i][j+1])%mod;
                if (j==0) dp[i+1][j]=dp[i][j+1];
                if (j>0 && j+1>a+b) dp[i+1][j]=dp[i][j-1];
            }
        cout<<dp[2*(a+b)][a]<<endl;
    }
    return 0;
}

F.Random Point in Triangle

直接跑蒙特卡洛搜索
当然也可以直接猜答案

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
struct Vector
{
    ll x,y;
};
ll area2(Vector a,Vector b)
{
    ll tmp=a.y*b.x-a.x*b.y;
    if (tmp<0) tmp=-tmp;
    return tmp;
}
int main()
{
    ll a,b,c,d,e,f;
    while (~scanf("%lld%lld%lld%lld%lld%lld",&a,&b,&c,&d,&e,&f))
    {
        Vector p,q;
        p.x=c-a;
        p.y=d-b;
        q.x=e-a;
        q.y=f-b;
        ll area=area2(p,q);
        ull ans=area;
        ans=ans*11;
        printf("%llu\n",ans);
    }
    return 0;
}

I.Points Division

答案一定是一条递增的曲线，不妨令这条曲线紧贴$B$。
线段树的下标是枚举的前一个点的高度，值是他的代价。
如果未来取得比$a[i].y$更低，那么这个点代价是$a$；否则代价是$b$

#include<iostream>
#include<algorithm>
#define fo(i,a,b) for(i=a;i<=b;i++)
#define N 100005
using namespace std;
struct www{long long x,y,a,b;} a[N];
long long tag[N*4],res[N*4],temp[N];
long long n,i,m,mx;
bool cmp(const www &p,const www &q)
{
    if (p.x != q.x) return p.x < q.x; else return p.y > q.y;
}
void build(int l,int r,int rt)
{
    tag[rt] = res[rt] = 0;
    if (l == r) return;
    int mid = (l + r) >> 1;
    build(l,mid,rt<<1); build(mid+1,r,rt<<1|1);
}
void pushdown(int rt)
{
    res[rt<<1] += tag[rt];
    res[rt<<1|1] += tag[rt];
    tag[rt<<1] += tag[rt];
    tag[rt<<1|1] += tag[rt];
    tag[rt] = 0;
}
long long query(int l,int r,int rt,int L,int R)
{
    if (L <= l && r <= R) return res[rt];
    pushdown(rt);
    int mid = (l + r) >> 1;
    long long temp = 0;
    if (L <= mid) temp = max(temp,query(l,mid,rt<<1,L,R));
    if (mid+1 <= R) temp = max(temp,query(mid+1,r,rt<<1|1,L,R));
    return temp;
}
void update1(int l,int r,int rt,int pos,long long val)
{
    if (l == r) {res[rt] = val; return;}
    pushdown(rt);
    int mid = (l + r) >> 1;
    if (pos <= mid) update1(l,mid,rt<<1,pos,val);
    if (mid+1 <= pos) update1(mid+1,r,rt<<1|1,pos,val);
    res[rt] = max(res[rt<<1],res[rt<<1|1]);
}
void update2(int l,int r,int rt,int L,int R,long long val)
{
    if (L <= l && r <= R)
    {
        res[rt] += val; tag[rt] += val;
        return;
    }
    pushdown(rt);
    int mid = (l + r) >> 1;
    if (L <= mid) update2(l,mid,rt<<1,L,R,val);
    if (mid+1 <= R) update2(mid+1,r,rt<<1|1,L,R,val);
    res[rt] = max(res[rt<<1],res[rt<<1|1]);
}
int main()
{
    while (~scanf("%d",&n))
    {
        fo(i,1,n)
            scanf("%d%d%d%d",&a[i].x,&a[i].y,&a[i].a,&a[i].b);
        sort(a+1,a+n+1,cmp);
        fo(i,1,n) temp[i] = a[i].y;
        sort(temp+1,temp+n+1);
        m = unique(temp+1,temp+n+1) - temp - 1;
        fo(i,1,n)
            a[i].y = lower_bound(temp+1,temp+m+1,a[i].y) - temp;
        m++;
        build(0,m,1);
        fo(i,1,n)
        {
            mx = query(0,m,1,0,a[i].y);
            update1(0,m,1,a[i].y,mx+a[i].b);
            update2(0,m,1,0,a[i].y-1,a[i].a);
            update2(0,m,1,a[i].y+1,m,a[i].b);
        }
        cout<<res[1]<<endl;
    }
    return 0;
}

J.Fraction Comparision

用int128比较即可