[考試反思]1029csp-s模擬測試92：弱智

 我只能這麼評價我本身。

看這個提交時間。。。我沒話可說。。。

T1半個世界都A了仍是切不掉。又一次掛細節。

T2不會證實的亂搞（雖然能夠證實）A了沒什麼可說的算是水過。

T3以前水過的題（打的次正解）而後不會二進制分組跑圖致使暴力都不會。

可是名次上可貴能看了一次，但其實遠不該該如此。

剩下的70分鐘仍是一點有效的正事都沒幹。。。

調整狀態，立刻就要11月份了啊。。。

 

T1：數列

ex_gcd。在ax+by=c有解的狀況下關於x總代價是三個一次函數，先降中間未知後升。

能夠三分。可是更簡單的就是找到xy的最小非負整數和最大負整數，一共4個點取最優。

這樣的話這4個點其實就是函數的兩個拐點最近的整點。

應該沒人像我同樣傻逼寫最小非負整數的時候沒有%p+p%p而是直接%p。我是傻子。

 1 #include<cstdio>
 2 #define int long long
 3 void ex_gcd(int a,int b,int &x,int &y){
 4     if(!b)return x=1,y=0,(void)0;
 5     ex_gcd(b,a%b,x,y);
 6     int r=x;x=y;y=r-a/b*x;
 7 }
 8 int gcd(int a,int b){return b?gcd(b,a%b):a;}
 9 int _abs(int x){return x>0?x:-x;}
10 int _min(int a,int b){return a<b?a:b;}
11 main(){
12     int n,a,b,p,g,ans=0,X,Y,x,y,tx,ty,ta;
13     scanf("%lld%lld%lld",&n,&a,&b);
14     g=gcd(a,b);a/=g;b/=g;
15     ex_gcd(a,b,X,Y);
16     while(n--){
17         scanf("%lld",&p);
18         if(p%g)return puts("-1"),0;
19         p/=g;x=X*p;y=Y*p;ta=1e18;
20         tx=(x%b+b)%b;ty=(p-tx*a)/b;ta=_min(_abs(tx)+_abs(ty),ta);
21         tx-= b;ty=(p-tx*a)/b;ta=_min(_abs(tx)+_abs(ty),ta);
22         ty=(y%a+a)%a;tx=(p-ty*b)/a;ta=_min(_abs(tx)+_abs(ty),ta);
23         ty-=a ;tx=(p-ty*b)/a;ta=_min(_abs(tx)+_abs(ty),ta);
24         ans+=ta;
25     }printf("%lld\n",ans);
26 }

View Code

 

 

 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迭代器遍從來直接離散化，沒鍋，挺帥的。

 1 #include<cstdio>
 2 #include<map>
 3 #include<algorithm>
 4 using namespace std;
 5 map<int,int>M;
 6 struct P{
 7     int a,b;long long w;
 8     friend bool operator<(P x,P y){
 9         int mx=max(x.b,x.a),my=max(y.b,y.a);
10         return mx<my||(mx==my&&x.b<y.b);
11     }
12 }p[100005];
13 struct Segment_Tree{
14     int cl[800005],cr[800005];long long w[800008],lz[800008];
15     void build(int p,int l,int r){
16         cl[p]=l;cr[p]=r;
17         if(l==r)return;
18         build(p<<1,l,l+r>>1);
19         build(p<<1|1,(l+r>>1)+1,r);
20     }
21     void down(int p){
22         lz[p<<1]+=lz[p];lz[p<<1|1]+=lz[p];
23         w[p<<1]+=lz[p];w[p<<1|1]+=lz[p];
24         lz[p]=0;
25     }
26     long long ask(int p,int l,int r){
27         if(l<=cl[p]&&cr[p]<=r)return w[p];
28         if(lz[p])down(p);
29         return max(l<=cr[p<<1]?ask(p<<1,l,r):0ll,r>cr[p<<1]?ask(p<<1|1,l,r):0ll);
30     }
31     void chg(int p,int pos,long long v){//if(p==1)printf("->:%d %lld\n",pos,v);
32         if(cl[p]==cr[p])return w[p]=max(w[p],v),(void)0;
33         if(lz[p])down(p);
34         if(pos<=cr[p<<1])chg(p<<1,pos,v);
35         else chg(p<<1|1,pos,v);
36         w[p]=max(w[p<<1],w[p<<1|1]);
37     }
38     void modify(int p,int l,int r,long long v){//if(p==1)printf("+:%d %d %lld\n",l,r,v);
39         if(l<=cl[p]&&cr[p]<=r)return lz[p]+=v,w[p]+=v,(void)0;
40         if(lz[p])down(p);
41         if(l<=cr[p<<1])modify(p<<1,l,r,v);
42         if(r>cr[p<<1])modify(p<<1|1,l,r,v);
43         w[p]=max(w[p<<1],w[p<<1|1]);
44     }
45 }T;
46 int main(){int n,cnt=0;//freopen("pair.in","r",stdin);
47     scanf("%d",&n);
48     for(int i=1;i<=n;++i)scanf("%d%d%lld",&p[i].a,&p[i].b,&p[i].w),M[p[i].a],M[p[i].b];
49     for(map<int,int>::iterator it=M.begin();it!=M.end();++it)(*it).second=++cnt;
50     for(int i=1;i<=n;++i)p[i].a=M[p[i].a],p[i].b=M[p[i].b];
51     T.build(1,1,cnt); sort(p+1,p+1+n);
52     //for(int i=1;i<=n;++i)printf("%d %d %lld\n",p[i].a,p[i].b,p[i].w);
53     for(int i=1;i<=n;++i)
54         if(p[i].a<p[i].b)T.modify(1,p[i].a+1,p[i].b,p[i].w),T.chg(1,p[i].a,T.ask(1,1,p[i].a)+p[i].w);
55         else T.chg(1,p[i].a,T.ask(1,1,p[i].b)+p[i].w);
56     printf("%lld\n",T.ask(1,1,cnt));
57 }

View Code

 

 

T3：最小距離

思惟嚴重僵化。Dijkstra原來能夠跑多源最短路啊

而後這題就沒了，全部特殊點爲源點跑最短路，記錄來源。

若是一條邊兩端的點來源不一樣，那麼就是把兩個源點連起來了，更新答案。

 1 #include<cstdio>
 2 #include<queue>
 3 using namespace std;
 4 #define LL long long
 5 LL min(LL a,LL b){return a<b?a:b;}
 6 struct P{
 7     int p,s;LL d;
 8     friend bool operator<(P a,P b){
 9         return a.d>b.d;
10     }
11 };priority_queue<P>q;
12 int n,m,Q,fir[200005],l[400005],to[400005],w[400005],cnt=1,ps[200005],tp,s[200005];
13 LL dt[200005],ans[200005];
14 void link(int a,int b,int v){l[++cnt]=fir[a];fir[a]=cnt;to[cnt]=b;w[cnt]=v;}
15 int main(){//freopen("distance.in","r",stdin);
16     scanf("%d%d%d",&n,&m,&Q);
17     for(int i=1;i<=n;++i)dt[i]=ans[i]=10000000000000000;
18     for(int i=1;i<=Q;++i)scanf("%d",&ps[i]),q.push((P){ps[i],s[ps[i]]=ps[i],dt[ps[i]]=0});
19     for(int i=1,a,b,v;i<=m;++i)scanf("%d%d%d",&a,&b,&v),link(a,b,v),link(b,a,v);
20     while(!q.empty()){
21         int p=q.top().p,S=q.top().s;LL d=q.top().d;q.pop();
22         if(d!=dt[p])continue;
23         for(int j=fir[p];j;j=l[j])if(dt[to[j]]>d+w[j])
24             q.push((P){to[j],s[to[j]]=S,dt[to[j]]=d+w[j]});
25     }
26     for(int i=1;i<=m;++i)if(s[to[i<<1]]!=s[to[i<<1|1]])
27         ans[s[to[i<<1|0]]]=min(ans[s[to[i<<1|0]]],dt[to[i<<1]]+w[i<<1]+dt[to[i<<1|1]]),
28         ans[s[to[i<<1|1]]]=min(ans[s[to[i<<1|1]]],dt[to[i<<1]]+w[i<<1]+dt[to[i<<1|1]]);
29     for(int i=1;i<=Q;++i)printf("%lld ",ans[ps[i]]);
30 }

View Code