2019.10.29 csp-s模拟测试92 反思总结

今天快乐的墨雨笙因为什么而几乎爆零了呢？

顾此失彼+不会对拍+无脑的复杂度

今天高兴的墨雨笙又因为什么调了一个下午呢？

不明题意+不想范围+板子低级错误

R.I.P.

T1：

//唉
//害怕TLE，所以没有用while
//害怕不正确，所以两种情况又算了两次
//细节全都考虑不清，何苦呢我……？
#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;
const long long inf=1e15;
long long sum,ans,n,a,b,x,y,xx,yy,x2,y2,flag;
long long exgcd(long long a,long long b,long long &x,long long &y){
    if(!b){
        x=,y=;
        return a;
    }
    long long q=exgcd(b,a%b,x,y);
    long long z=x;
    x=y,y=z-(a/b)*y;
    return q;
}
int main()
{
//    freopen("1.in","r",stdin);
    scanf("%lld%lld%lld",&n,&a,&b);
    long long gcd=exgcd(a,b,x2,y2);
    if(b/gcd<a/gcd)swap(a,b),swap(x2,y2);//!
    for(int i=;i<=n;i++){
        long long val;
        scanf("%lld",&val);
        if(flag)continue;
        if(val<)val=-val;//不，这个不用考虑也可以
        ans=inf;
        if(val%a==)ans=min(ans,val/a);
        if(val%b==)ans=min(ans,val/b);
        x=x2,y=y2;
        if(val%gcd==){
            x*=val/gcd,y*=val/gcd;
            long long b0=b/gcd,a0=a/gcd;
            xx=(x%b0+b0)%b0;//!
            yy=(val-xx*a)/b;
            while(abs(xx-b0)+abs(yy+a0)<abs(xx)+abs(yy))xx-=b0,yy+=a0;//!
            while(abs(xx+b0)+abs(yy-a0)<abs(xx)+abs(yy))xx+=b0,yy-=a0;
            ans=min(ans,(long long)abs(xx)+(long long)abs(yy));
            //实际上只有四个地方的点有可能产生贡献：
            //x为最小正整数解、x为最小负整数解、y为最小正整数解、y为最小负整数解
            //否则一定可以继续缩小
        }
        if(ans==inf){
            flag=;
        }
        else sum+=ans;
    }
    if(flag)printf("-1");
    else printf("%lld\n",sum);
    return ;
 }

T2：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=3e5+;
int n,c[*N],cnt,T[N*],L[N*],R[N*],tot,lst;
long long f[N],val,maxx[N],sum[N*],ans;
struct node{
    int a,b,w;
}s[N];
bool cmp(node x,node y){
    return x.a==y.a?x.b<y.b:x.a<y.a;
}
void build(int &p,int l,int r){
    p=++tot;
    if(l==r)return;
    int mid=(l+r)/;
    build(L[p],l,mid);
    build(R[p],mid+,r);
}
long long query(int p0,int p,int l,int r,int l0,int r0){
    if(l0<=l&&r<=r0){
        return sum[p]-sum[p0];
    }
    int mid=(l+r)/;
    if(r0<=mid)return query(L[p0],L[p],l,mid,l0,r0);
    else if(l0>mid)return query(R[p0],R[p],mid+,r,l0,r0);
    else{
        return query(L[p0],L[p],l,mid,l0,mid)+query(R[p0],R[p],mid+,r,mid+,r0);
    }
}
void update(int p0,int &p,int l,int r,int pos,int val){
    p=++tot;
    sum[p]=sum[p0]+val,L[p]=L[p0],R[p]=R[p0];
    if(l==r)return;
    int mid=(l+r)/;
    if(pos<=mid)update(L[p0],L[p],l,mid,pos,val);
    else update(R[p0],R[p],mid+,r,pos,val);
}
int main()
{
//    freopen("pair.in","r",stdin);
    scanf("%d",&n);
    for(int i=;i<=n;i++){
        scanf("%d%d%d",&s[i].a,&s[i].b,&s[i].w);
        c[++cnt]=s[i].a,c[++cnt]=s[i].b;
    }
    sort(c+,c+cnt+);
    cnt=unique(c+,c+cnt+)-c-;
    sort(s+,s+n+,cmp);
    build(T[],,cnt);
    for(int i=;i<=n;i++){
        int a1=lower_bound(c+,c+cnt+,s[i].a)-c;
        int b1=lower_bound(c+,c+cnt+,s[i].b)-c;
        while(lst<a1){
            maxx[lst+]=maxx[lst];
            T[lst+]=T[lst];
            lst++;
        }
        if(b1<a1)val=query(T[b1],T[a1],,cnt,a1,cnt);
        else val=;
        f[i]=maxx[min(b1,a1)]+val+s[i].w;
        maxx[a1]=max(maxx[a1],f[i]);
        update(T[a1],T[a1],,cnt,b1,s[i].w);
        ans=max(ans,f[i]);
    }
    printf("%lld\n",ans);
    return ;
}

个人玄学版本

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=3e5+;
int n,cnt,T[N*],L[N*],R[N*],tot;
long long f[N],c[*N],val,maxx[N],sum[N*],ans,aa[N];
struct node{
    long long a,b,w,a1,b1;
}s[N];
bool cmp(node x,node y){
    return x.a<y.a;
}
void build(int &p,int l,int r){
    p=++tot;
    sum[p]=;
    if(l==r)return;
    int mid=(l+r)/;
    build(L[p],l,mid);
    build(R[p],mid+,r);
}
long long query(int p0,int p,int l,int r,int l0,int r0){
    if(l0<=l&&r<=r0){
        return sum[p]-sum[p0];
    }
    int mid=(l+r)/;
    if(r0<=mid)return query(L[p0],L[p],l,mid,l0,r0);
    else if(l0>mid)return query(R[p0],R[p],mid+,r,l0,r0);
    else{
        return query(L[p0],L[p],l,mid,l0,mid)+query(R[p0],R[p],mid+,r,mid+,r0);
    }
}
void update(int p0,int &p,int l,int r,long long pos,long long val){
    p=++tot;
    sum[p]=sum[p0]+val,L[p]=L[p0],R[p]=R[p0];
    if(l==r)return;
    int mid=(l+r)/;
    if(pos<=mid)update(L[p0],L[p],l,mid,pos,val);
    else update(R[p0],R[p],mid+,r,pos,val);
}
int main()
{
//    freopen("pair.in","r",stdin);
    scanf("%d",&n);
    for(int i=;i<=n;i++){
        scanf("%lld%lld%lld",&s[i].a,&s[i].b,&s[i].w);
        c[++cnt]=s[i].a,c[++cnt]=s[i].b;
    }
    sort(c+,c+cnt+);
    cnt=unique(c+,c+cnt+)-c-;
    sort(s+,s+n+,cmp);
    build(T[],,cnt);
    for(int i=;i<=n;i++){
        long long b1=lower_bound(c+,c+cnt+,s[i].b)-c;
        update(T[i-],T[i],,cnt,b1,s[i].w);
        aa[i]=s[i].a;
    }
    for(int i=;i<=n;i++){
        long long a1=lower_bound(c+,c+cnt+,s[i].a)-c;
        long long b1=lower_bound(c+,c+cnt+,s[i].b)-c;
        int pos=upper_bound(aa+,aa+n+,s[i].b)-aa-;
        if(pos>=i)pos=i-;
        val=query(T[pos],T[i-],,cnt,a1,cnt);
        f[i]=maxx[pos]+val+s[i].w;
        maxx[i]=max(maxx[i-],f[i]);
        ans=max(ans,f[i]);
    }
    printf("%lld\n",ans);
    return ;
}

gls指点以后的清晰版本

T3：

#include<iostream>
#include<cstdio>
#include<queue>
using namespace std;
const int N=2e5+;
const long long inf=1e18;
int n,m,p;
int a[N],vis[N],tag[N];
int ver[*N],Next[*N],head[N],tot,xx[*N];
long long dis[N],ans[N],edge[*N];
priority_queue<pair<long long,int> >q;
void add(int x,int y,long long z){
    ver[++tot]=y;
    Next[tot]=head[x];
    head[x]=tot;
    edge[tot]=z;
    xx[tot]=x;
}
void dij(){
    for(int i=;i<=n;i++){
        if(!dis[i])q.push(make_pair(,i));
    }
    while(!q.empty()){
        while(vis[q.top().second]&&q.size())q.pop();
        if(q.empty())break;
        int x=q.top().second;
        q.pop();
        vis[x]=;
        for(int i=head[x];i;i=Next[i]){
            int y=ver[i];
            long long z=edge[i];
            if(!vis[y]&&dis[x]+z<dis[y]){
                dis[y]=dis[x]+z;
                tag[y]=tag[x];
                q.push(make_pair(-dis[y],y));
            }
        }
    }
}
int main()
{
//    freopen("1.in","r",stdin);
    scanf("%d%d%d",&n,&m,&p);
    for(int i=;i<=n;i++)dis[i]=inf;
    for(int i=;i<=p;i++){
        scanf("%d",&a[i]);
        dis[a[i]]=;
        tag[a[i]]=a[i];
    }
    for(int i=;i<=n;i++)ans[i]=inf;
    for(int i=,x,y;i<=m;i++){
        long long z;
        scanf("%d%d%lld",&x,&y,&z);
        add(x,y,z),add(y,x,z);
    }
    dij();
    for(int i=;i<=tot;i+=){
        int x=xx[i],y=ver[i];
        if(tag[x]!=tag[y]){
            long long d=dis[x]+dis[y]+edge[i];
            ans[tag[x]]=min(ans[tag[x]],d);
            ans[tag[y]]=min(ans[tag[y]],d);
        }
    }
    for(int i=;i<=p;i++)printf("%lld ",ans[a[i]]);
    return ;
}

别问，问就是死不足惜【大笑】