[考試反思]1103csp-s模擬測試99：好夢

時間 2019-11-11

標籤考試反思 1103csp csp 模擬測試简体版

原文原文鏈接

可能此次考得好的緣由就是熬夜頹廢到不算太晚？（啪）html

可是是真心困。c++

考前跟akt說：我但願今天考一點那種不用動腦子，就是一直碼的題。app

而後開門T1一道線段樹維護單調棧的板子我就。。。了ide

當時調了一上午啊，我就以爲考場上不可能調出來，可是T2和T3都不會作，因此回去調板子。測試

說是「調」其實很不嚴謹，由於我其實就是一直在寫，連編譯都沒有編譯過。優化

一個半小時以後差很少要寫完了（可是尚未寫完）。感受心累，並不以爲本身寫的是對的（模板爆零專家）ui

因而沒打完就棄了。spa

T3感受數據很難造而可行方案不少，因而亂打了一個。3d

而後以爲T2的題目有7個參數很難寫，就想都沒想打了一個很是蠢的$O(n^3log n)$的暴力。（能夠很是簡單的優化到$O(n^2 log n)$）code

最後給T1寫了一個20分的暴力，把本身寫的線段樹註釋了。離考試結束還有10分鐘。

檢查完文件閒得無聊。因而去繼續在註釋裏寫T1的線段樹。

原本就沒差多少因而寫完了，腦子一抽想試一下。

結果就過樣例了，然而我很不放心以爲它過不了，因而還測試點分治了一下。

而後我板子竟然沒爆炸2333太快樂了。

而後T3亂搞A了2333太快樂了。

而後T2的$O(n^3logn)$過了1000的點2333太快樂了。

而後就莫名其妙地rank1了2333太快樂了。

真的不是在用腦子考試，是在用RP考試。。。

T1：陶陶摘蘋果

題目已經說了顯然是快速傅麗葉變換再加一個狄利克雷卷積

能夠發現這是一個很是裸的線段樹優化單調棧，能夠去看我《椎》的題解，寫的挺詳細的。

複雜度是$O(nlog^2n)$的。不是最優的可是也可過。

 1 #include<bits/stdc++.h>
 2 using namespace std;
 3 int cl[800005],cr[800005],w[800005],mx[800005],mxr,n,m,h[100005],rv[800005];
 4 int ask(int p,int l,int r,int W){
 5     if(W>=mx[p])return 0;
 6     if(cl[p]==cr[p]){mxr=max(mxr,mx[p]);return 1;}
 7     if(l<=cl[p]&&cr[p]<=r){
 8         if(W>mx[p<<1|1])return ask(p<<1,l,r,W);
 9         int re=ask(p<<1|1,l,r,W)+rv[p];
10         mxr=max(mxr,mx[p]);return re;
11     }
12     int ans=0;
13     if(r>cr[p<<1])ans=ask(p<<1|1,l,r,W);
14     if(l<=cr[p<<1])ans+=ask(p<<1,l,r,mxr);
15     return ans;
16 }
17 void update(int p){
18     mx[p]=max(mx[p<<1],mx[p<<1|1]);
19     mxr=0;rv[p]=ask(p<<1,cl[p<<1],cr[p<<1],mx[p<<1|1]);
20 }
21 void modify(int p,int pos,int H){
22     if(cl[p]==cr[p]){mx[p]=H;return;}
23     modify(pos<=cr[p<<1]?p<<1:p<<1|1,pos,H);
24     update(p);
25 }
26 void build(int p,int l,int r){
27     cl[p]=l;cr[p]=r;
28     if(l==r){w[p]=1;mx[p]=h[l];return;}
29     build(p<<1,l,l+r>>1);
30     build(p<<1|1,(l+r>>1)+1,r);
31     update(p);
32 }
33 int main(){
34     freopen("TaoPApp.in","r",stdin);
35     freopen("TaoPApp.out","w",stdout);
36     scanf("%d%d",&n,&m);
37     for(int i=1;i<=n;++i)scanf("%d",&h[n+1-i]);
38     build(1,1,n);
39     for(int i=1;i<=m;++i){
40         int p,nw,r,mx=-1,cnt=0;scanf("%d%d",&p,&nw);p=n+1-p;r=h[p];
41         if(n<=5000&&m<=5000){
42             h[p]=nw;
43             for(int j=n;j;--j)if(h[j]>mx)mx=h[j],cnt++;
44             printf("%d\n",cnt);h[p]=r;
45         }
46         else modify(1,p,nw),mxr=0,printf("%d\n",ask(1,1,n,0)),modify(1,p,r);
47     }
48 }

View Code

也有其它作法，如分治（像《走路》）。

題解的作法是考慮每一個點左右兩端，預處理出來它的相似與於單調棧的東西，而後二分位置。也挺麻煩的。

T2：開心的金明

暴力50%：

枚舉買進原料的時間，枚舉製造電腦的時間，枚舉賣出的時間，線段樹維護電腦還能存多少（區間減，區間取min），複雜度O(n^3logn)。

 1 #include<cstdio>
 2 #define inf 1234567890
 3 inline int min(int a,int b){return a<b?a:b;}
 4 int k,bf[50005],s[50005],kpn[50005],kpf[50005],ksf[50005],mf[50005],mn[50005];long long ans;
 5 struct Segment_Tree{
 6     int mn[400005],lz[400005],cl[400005],cr[400005];
 7     void build(int p,int l,int r){
 8         cl[p]=l;cr[p]=r;
 9         if(l==r)return (void)(mn[p]=kpn[l]);
10         build(p<<1,l,l+r>>1);
11         build(p<<1|1,(l+r>>1)+1,r);
12         mn[p]=min(mn[p<<1],mn[p<<1|1]);
13     }
14     void down(int p){
15         mn[p<<1]-=lz[p];mn[p<<1|1]-=lz[p];
16         lz[p<<1]+=lz[p];lz[p<<1|1]+=lz[p];
17         lz[p]=0;
18     }
19     void minus(int p,int l,int r,int w){
20         if(r<l)return;
21         if(l<=cl[p]&&cr[p]<=r)return mn[p]-=w,lz[p]+=w,(void)0;
22         if(lz[p])down(p);
23         if(l<=cr[p<<1])minus(p<<1,l,r,w);
24         if(r>cr[p<<1])minus(p<<1|1,l,r,w);
25         mn[p]=min(mn[p<<1],mn[p<<1|1]);
26     }
27     int ask(int p,int l,int r){
28         if(r<l)return inf;
29         if(l<=cl[p]&&cr[p]<=r)return mn[p];
30         if(lz[p])down(p);
31         return min(l<=cr[p<<1]?ask(p<<1,l,r):inf,r>cr[p<<1]?ask(p<<1|1,l,r):inf);
32     }
33 }T;
34 int main(){
35     freopen("happy.in","r",stdin);
36     freopen("happy.out","w",stdout);
37     scanf("%d",&k);
38     for(int i=1;i<=k;++i)scanf("%d%d%d%d",&bf[i],&s[i],&mf[i],&mn[i]);
39     for(int i=1;i<k;++i)scanf("%d%d%d",&kpn[i],&ksf[i],&kpf[i]),ksf[i]+=ksf[i-1],kpf[i]+=kpf[i-1];
40     T.build(1,1,k-1);
41     for(int sd=1;sd<=k;++sd){
42         int lft=s[sd];
43         while(lft){
44             int bbd,bpd,num,cost=inf;
45             for(int pd=sd;pd;--pd)for(int bd=pd;bd;--bd){
46                 int C=bf[bd]+ksf[pd-1]-ksf[bd-1]+kpf[sd-1]-kpf[pd-1]+mf[pd],Num=min(mn[pd],T.ask(1,pd,sd-1));
47                 if(C<cost&&Num)num=Num,cost=C,bbd=bd,bpd=pd;
48             }
49             if(cost==inf)return puts("-1"),0;
50             num=min(num,lft);
51             lft-=num;T.minus(1,bpd,sd-1,num);mn[bpd]-=num;
52             ans+=num*cost;
53         }
54     }printf("%lld\n",ans);
55 }

View Code

由於對材料沒有任何的限制，因此能夠預處理材料在每一天的價格，就能夠下降到$O(n^2log n)$。

直接上正解。

首先正解也用到了這種處理原材料的價格的方法：對於新的一天，要麼把昨天單位原材料的價格+昨天到今天存儲原材料的價格做爲今天的價格，要麼直接把今天買入的價格做爲真正價格。

而後反正你已經有原料了。用set維護結構體存每個電腦的費用以及這個費用的電腦的數量。

模擬每一天的過程。

計算今天製造電腦的價格（最小原材料價格+製造費用），放進set裏。

賣電腦。從set裏取出價格最小的電腦一直到賣夠數量，累加答案，賣不夠就無解。

屯材料。計算材料留到明天的價格，是min(留到今天的價格+存儲費用，明天再買的價格)

屯電腦。若是set裏的電腦多於存儲上限，那就幹掉最貴的，留下的電腦漲價（存儲價格）。

沒了。就是貪心。

 1 #include<cstdio>
 2 #include<set>
 3 using namespace std;
 4 struct P{int n,f;friend bool operator<(P a,P b){return a.f<b.f;}};
 5 multiset<P>S,res;
 6 int k,bf[50005],s[50005],kpn[50005],kpf[50005],ksf[50005],mf[50005],mn[50005],sf,totn;
 7 long long ans;
 8 main(){
 9     freopen("happy.in","r",stdin);
10     freopen("happy.out","w",stdout);
11     scanf("%d",&k);
12     for(int i=1;i<=k;++i)scanf("%d%d%d%d",&bf[i],&s[i],&mf[i],&mn[i]);
13     for(int i=1;i<k;++i)scanf("%d%d%d",&kpn[i],&ksf[i],&kpf[i]);
14     sf=bf[1];
15     for(int i=1;;++i){
16         if(mn[i])S.insert((P){mn[i],mf[i]+sf});totn+=mn[i];
17         while(s[i]){
18             if(S.empty())return puts("-1"),0;
19             P x=*S.begin();int num=min(s[i],x.n);S.erase(S.begin());
20             ans+=1ll*x.f*num;x.n-=num;totn-=num;s[i]-=num;
21             if(x.n)S.insert((P){x.n,x.f});
22         }
23         if(i==k)break;
24         sf=min(sf+ksf[i],bf[i+1]);
25         while(totn>kpn[i]){
26             P x=*S.rbegin();int num=min(totn-kpn[i],x.n);S.erase(--S.end());
27             totn-=num;x.n-=num;
28             if(x.n)S.insert((P){x.n,x.f});
29         }
30         for(auto it=S.begin();it!=S.end();++it)res.insert((P){(*it).n,(*it).f+kpf[i]});
31         swap(S,res);res.clear();
32     }printf("%lld\n",ans);
33 }

View Code

T3：笨小猴

 1 #include<cstdio>
 2 #include<set>
 3 using namespace std;
 4 struct A{int a,b,o;friend bool operator<(A x,A y){return x.a<y.a;}};
 5 struct B{int a,b,o;friend bool operator<(B x,B y){return x.b<y.b;}};
 6 set<A>as[2];set<B>bs[2];
 7 int st[200005],n;long long ala,alb,tota,totb;
 8 int main(){
 9     freopen("grandmaster.in","r",stdin);
10     freopen("grandmaster.out","w",stdout);
11     scanf("%d",&n);
12     for(int i=1,a,b;i<=n+1;++i){
13         scanf("%d%d",&a,&b);
14         ala+=a;alb+=b;tota+=a;totb+=b;
15         as[0].insert((A){a,b,i});
16         bs[0].insert((B){a,b,i});
17     }
18     for(int i=1,a,b;i<=n;++i){
19         scanf("%d%d",&a,&b);
20         tota+=a;totb+=b;st[n+1+i]=1;
21         as[1].insert((A){a,b,n+1+i});
22         bs[1].insert((B){a,b,n+1+i});
23     }
24     if(ala>tota-ala&&alb>totb-alb){for(int i=1;i<=n+1;++i)printf("%d\n",i);return 0;}
25     if(ala<=tota-ala&&alb<=totb-alb){for(int i=n;i<=n+n+1;++i)printf("%d\n",i);return 0;}
26     int nst=0;
27     while(1){
28         ala=tota-ala;alb=totb-alb;//printf("%lld %lld\n",ala,alb);
29         if(ala<=tota-ala){
30             A c=*(--as[nst].end());B d=(B){c.a,c.b,c.o};//printf("--%d %d %d\n",c.a,c.b,c.o);
31             as[nst].erase(c);bs[nst].erase(d);
32             nst^=1;as[nst].insert(c);bs[nst].insert(d);
33             ala+=c.a;alb+=c.b;st[c.o]=nst;
34         }else{
35             B c=*(--bs[nst].end());A d=(A){c.a,c.b,c.o};//printf("--%d %d %d\n",c.a,c.b,c.o);
36             bs[nst].erase(c);as[nst].erase(d);
37             nst^=1;bs[nst].insert(c);as[nst].insert(d);
38             ala+=c.a;alb+=c.b;st[c.o]=nst;
39         }
40         if(ala>tota-ala&&alb>totb-alb){for(int i=1;i<=n+n+1;++i)if(st[i]==nst)printf("%d\n",i);return 0;}
41     }
42 }