[考试反思]1029csp-s模拟测试92：弱智

我只能这么评价我自己。

看这个提交时间。。。我没话可说。。。

T1半个世界都A了还是切不掉。又一次挂细节。

T2不会证明的乱搞（虽然可以证明）A了没什么可说的算是水过。

T3之前水过的题（打的次正解）然后不会二进制分组跑图导致暴力都不会。

但是名次上难得能看了一次，但其实远不应该如此。

剩下的70分钟还是一点有效的正事都没干。。。

调整状态，马上就要11月份了啊。。。

T1：数列

ex_gcd。在ax+by=c有解的情况下关于x总代价是三个一次函数，先降中间未知后升。

可以三分。但是更简单的就是找到xy的最小非负整数和最大负整数，一共4个点取最优。

这样的话这4个点其实就是函数的两个拐点最近的整点。

应该没人像我一样傻逼写最小非负整数的时候没有%p+p%p而是直接%p。我是傻子。

 #include<cstdio>
 #define int long long
 void ex_gcd(int a,int b,int &x,int &y){
     if(!b)return x=,y=,(void);
     ex_gcd(b,a%b,x,y);
     int r=x;x=y;y=r-a/b*x;
 }
 int gcd(int a,int b){return b?gcd(b,a%b):a;}
 int _abs(int x){return x>?x:-x;}
 int _min(int a,int b){return a<b?a:b;}
 main(){
     int n,a,b,p,g,ans=,X,Y,x,y,tx,ty,ta;
     scanf("%lld%lld%lld",&n,&a,&b);
     g=gcd(a,b);a/=g;b/=g;
     ex_gcd(a,b,X,Y);
     while(n--){
         scanf("%lld",&p);
         if(p%g)return puts("-1"),;
         p/=g;x=X*p;y=Y*p;ta=1e18;
         tx=(x%b+b)%b;ty=(p-tx*a)/b;ta=_min(_abs(tx)+_abs(ty),ta);
         tx-= b;ty=(p-tx*a)/b;ta=_min(_abs(tx)+_abs(ty),ta);
         ty=(y%a+a)%a;tx=(p-ty*b)/a;ta=_min(_abs(tx)+_abs(ty),ta);
         ty-=a ;tx=(p-ty*b)/a;ta=_min(_abs(tx)+_abs(ty),ta);
         ans+=ta;
     }printf("%lld\n",ans);
 }

T2：数对

如果不考虑可以换顺序的话，就是《队长快跑》带权值的原题。

线段树优化dp，单点修改，区间加，区间查询max。

考虑怎么换顺序。

两个点会相互影响决策，要么是$ a_i < b_j $且$ b_i < a_j $这样的话$i$必须在$j$前面。（题解的符号写反了）

还有一种就是完全相反。

再一种就是$ a_i < b_j $且$ a_j < b_i $这样的话两个互相包容都可以选，换顺序没有影响。

所以其实只有前两种情况真正限制了它们的关系。

所以直接按照a+b排序就能保证上述顺序对的点顺序是对的。

也可以用max/min比较一下，不过这样的话需要第二关键字，随便弄一个都对。

试了一下用map迭代器遍历来直接离散化，没锅，挺帅的。

 #include<cstdio>
 #include<map>
 #include<algorithm>
 using namespace std;
 map<int,int>M;
 struct P{
     int a,b;long long w;
     friend bool operator<(P x,P y){
         int mx=max(x.b,x.a),my=max(y.b,y.a);
         return mx<my||(mx==my&&x.b<y.b);
     }
 }p[];
 struct Segment_Tree{
     int cl[],cr[];long long w[],lz[];
     void build(int p,int l,int r){
         cl[p]=l;cr[p]=r;
         if(l==r)return;
         build(p<<,l,l+r>>);
         build(p<<|,(l+r>>)+,r);
     }
     void down(int p){
         lz[p<<]+=lz[p];lz[p<<|]+=lz[p];
         w[p<<]+=lz[p];w[p<<|]+=lz[p];
         lz[p]=;
     }
     long long ask(int p,int l,int r){
         if(l<=cl[p]&&cr[p]<=r)return w[p];
         if(lz[p])down(p);
         return max(l<=cr[p<<]?ask(p<<,l,r):0ll,r>cr[p<<]?ask(p<<|,l,r):0ll);
     }
     void chg(int p,int pos,long long v){//if(p==1)printf("->:%d %lld\n",pos,v);
         if(cl[p]==cr[p])return w[p]=max(w[p],v),(void);
         if(lz[p])down(p);
         if(pos<=cr[p<<])chg(p<<,pos,v);
         else chg(p<<|,pos,v);
         w[p]=max(w[p<<],w[p<<|]);
     }
     void modify(int p,int l,int r,long long v){//if(p==1)printf("+:%d %d %lld\n",l,r,v);
         if(l<=cl[p]&&cr[p]<=r)return lz[p]+=v,w[p]+=v,(void);
         if(lz[p])down(p);
         if(l<=cr[p<<])modify(p<<,l,r,v);
         if(r>cr[p<<])modify(p<<|,l,r,v);
         w[p]=max(w[p<<],w[p<<|]);
     }
 }T;
 int main(){int n,cnt=;//freopen("pair.in","r",stdin);
     scanf("%d",&n);
     for(int i=;i<=n;++i)scanf("%d%d%lld",&p[i].a,&p[i].b,&p[i].w),M[p[i].a],M[p[i].b];
     for(map<int,int>::iterator it=M.begin();it!=M.end();++it)(*it).second=++cnt;
     for(int i=;i<=n;++i)p[i].a=M[p[i].a],p[i].b=M[p[i].b];
     T.build(,,cnt); sort(p+,p++n);
     //for(int i=1;i<=n;++i)printf("%d %d %lld\n",p[i].a,p[i].b,p[i].w);
     for(int i=;i<=n;++i)
         if(p[i].a<p[i].b)T.modify(,p[i].a+,p[i].b,p[i].w),T.chg(,p[i].a,T.ask(,,p[i].a)+p[i].w);
         else T.chg(,p[i].a,T.ask(,,p[i].b)+p[i].w);
     printf("%lld\n",T.ask(,,cnt));
 }

T3：最小距离

思维严重僵化。Dijkstra原来可以跑多源最短路啊

然后这题就没了，所有特殊点为源点跑最短路，记录来源。

如果一条边两端的点来源不同，那么就是把两个源点连起来了，更新答案。

 #include<cstdio>
 #include<queue>
 using namespace std;
 #define LL long long
 LL min(LL a,LL b){return a<b?a:b;}
 struct P{
     int p,s;LL d;
     friend bool operator<(P a,P b){
         return a.d>b.d;
     }
 };priority_queue<P>q;
 int n,m,Q,fir[],l[],to[],w[],cnt=,ps[],tp,s[];
 LL dt[],ans[];
 void link(int a,int b,int v){l[++cnt]=fir[a];fir[a]=cnt;to[cnt]=b;w[cnt]=v;}
 int main(){//freopen("distance.in","r",stdin);
     scanf("%d%d%d",&n,&m,&Q);
     for(int i=;i<=n;++i)dt[i]=ans[i]=;
     for(int i=;i<=Q;++i)scanf("%d",&ps[i]),q.push((P){ps[i],s[ps[i]]=ps[i],dt[ps[i]]=});
     for(int i=,a,b,v;i<=m;++i)scanf("%d%d%d",&a,&b,&v),link(a,b,v),link(b,a,v);
     while(!q.empty()){
         int p=q.top().p,S=q.top().s;LL d=q.top().d;q.pop();
         if(d!=dt[p])continue;
         for(int j=fir[p];j;j=l[j])if(dt[to[j]]>d+w[j])
             q.push((P){to[j],s[to[j]]=S,dt[to[j]]=d+w[j]});
     }
     for(int i=;i<=m;++i)if(s[to[i<<]]!=s[to[i<<|]])
         ans[s[to[i<<|]]]=min(ans[s[to[i<<|]]],dt[to[i<<]]+w[i<<]+dt[to[i<<|]]),
         ans[s[to[i<<|]]]=min(ans[s[to[i<<|]]],dt[to[i<<]]+w[i<<]+dt[to[i<<|]]);
     for(int i=;i<=Q;++i)printf("%lld ",ans[ps[i]]);
 }