暑假集训test-8-31(pm)

以为可以AK，结果t3没有调出来，然后被林巨踩了。

everyday被踩，很开心。

林巨真的好巨啊，这么多天已经总计虐我75分了。

1.玩具装箱

第一眼还以为是那道斜率优化dp,结果是个签到水题。

 //Achen
 #include<algorithm>
 #include<iostream>
 #include<cstring>
 #include<cstdlib>
 #include<vector>
 #include<cstdio>
 #include<queue>
 #include<cmath>
 #include<set>
 #include<map>
 #define Formylove return 0
 #define For(i,a,b) for(int i=(a);i<=(b);i++)
 #define Rep(i,a,b) for(int i=(a);i>=(b);i--)
 const int N=;
 typedef long long LL;
 typedef double db;
 using namespace std;
 int n,m;
 LL k,f[N],a[N];

 template<typename T>void read(T &x)  {
     char ch=getchar(); x=; T f=;
     while(ch!='-'&&(ch<''||ch>'')) ch=getchar();
     if(ch=='-') f=-,ch=getchar();
     for(;ch>=''&&ch<='';ch=getchar()) x=x*+ch-''; x*=f;
 }

 #define ANS
 int main() {
 #ifdef ANS
     freopen("toy.in","r",stdin);
     freopen("toy.out","w",stdout);
 #endif
     read(n); read(m); read(k);
     For(i,,n) read(a[i]);
     memset(f,/,sizeof(f));
     f[]=;
     For(i,,n) {
         LL mi=a[i],mx=a[i];
         Rep(j,i-,max(,i-m)) {
             f[i]=min(f[i],f[j]+k+(mx-mi)*(i-j));
             mi=min(mi,a[j]);
             mx=max(mx,a[j]);
         }
     }
     printf("%lld\n",f[n]);
     Formylove;
 }

2.铁路运费

林巨写的比较优美，我比较傻就用set水过去了。

边权从1到2等于删边。

bfs树上树边和指向深度比它小的点的边有用，分别用set存一下。每次删树边的时候在另一个set里找还有没有可以用的父亲，有就连成新的树边，退出，没有就把它标记了然后删除它到它儿子的边，递归操作它儿子。

 //Achen
 #include<algorithm>
 #include<iostream>
 #include<cstring>
 #include<cstdlib>
 #include<vector>
 #include<cstdio>
 #include<queue>
 #include<cmath>
 #include<set>
 #include<map>
 #define Formylove return 0
 #define For(i,a,b) for(int i=(a);i<=(b);i++)
 #define Rep(i,a,b) for(int i=(a);i>=(b);i--)
 const int N=2e5+;
 typedef long long LL;
 typedef double db;
 using namespace std;
 int n,m,q,ans,no[N];

 template<typename T>void read(T &x)  {
     char ch=getchar(); x=; T f=;
     while(ch!='-'&&(ch<''||ch>'')) ch=getchar();
     if(ch=='-') f=-,ch=getchar();
     for(;ch>=''&&ch<='';ch=getchar()) x=x*+ch-''; x*=f;
 }

 struct edge {
     int u,v;
 }e[N];

 int ecnt,fir[N],nxt[N<<],to[N<<];
 void init() { ecnt=; memset(fir,,sizeof(fir)); }
 void add(int u,int v) {
     nxt[++ecnt]=fir[u]; fir[u]=ecnt; to[ecnt]=v;
     nxt[++ecnt]=fir[v]; fir[v]=ecnt; to[ecnt]=u;
 }

 int d[N],fid[N];
 queue<int>que;
 set<int>s[N],son[N];

 void bfs(int s) {
     For(i,,n) d[i]=n;
     d[s]=;
     que.push(s);
     while(!que.empty()) {
         int x=que.front();
         que.pop();
         for(int i=fir[x];i;i=nxt[i]) if(d[to[i]]==n) {
             d[to[i]]=d[x]+;

             son[x].insert(to[i]);
             que.push(to[i]);
         }
     }
 }

 void solve(int x,int fa) {
     son[fa].erase(x);
     s[x].erase(fa);
     while(s[x].size()) {
         int y=*s[x].begin();
         if(!no[y]) { son[y].insert(x); return; }
         s[x].erase(y);
     }
     ans++;
     no[x]=;
     while(son[x].size()) {
         int y=*son[x].begin();
         solve(y,x);
     }
 }

 #define ANS
 int main() {
 #ifdef ANS
     freopen("train.in","r",stdin);
     freopen("train.out","w",stdout);
 #endif
     read(n); read(m); read(q);
     For(i,,m) {
         read(e[i].u); read(e[i].v);
         add(e[i].u,e[i].v);
     }
     bfs();
     For(x,,n) {
         for(int i=fir[x];i;i=nxt[i]) if(d[to[i]]+==d[x])
             s[x].insert(to[i]);
     }
     For(ti,,q) {
         int id;
         read(id);
         int u=e[id].u,v=e[id].v;
         if(son[u].find(v)!=son[u].end()||son[v].find(u)!=son[v].end()) {
             if(d[u]+==d[v]) solve(v,u);
             else solve(u,v);
         }
         else {
             if(s[v].find(u)!=s[v].end()) s[v].erase(u);
             else if(s[u].find(v)!=s[u].end()) s[u].erase(v);
         }
         printf("%d\n",ans);
     }
     Formylove;
 }
 /*
 5 6 5
 1 2
 1 3
 4 2
 3 2
 2 5
 5 3
 5
 2
 4
 1
 3
 */

3.城墙

对于每个点维护它往右往下能到的最长没有x的长度记作a[i][j]，往左往上能到达的没有x的最长长度b[i][j]。沿着对角线扫，用树状数组维护一下，把a放进树状数组，用b和L限制查询范围。

 //Achen
 #include<algorithm>
 #include<iostream>
 #include<cstring>
 #include<cstdlib>
 #include<vector>
 #include<cstdio>
 #include<queue>
 #include<cmath>
 #include<set>
 #include<map>
 #define Formylove return 0
 #define For(i,a,b) for(int i=(a);i<=(b);i++)
 #define Rep(i,a,b) for(int i=(a);i>=(b);i--)
 const int N=;
 typedef long long LL;
 typedef double db;
 using namespace std;
 int n,m,L,tot;
 int mp[N][N],s[N][N],x[N][N],z[N][N],y[N][N],a[N][N],b[N][N];

 template<typename T>void read(T &x)  {
     char ch=getchar(); x=; T f=;
     while(ch!='-'&&(ch<''||ch>'')) ch=getchar();
     if(ch=='-') f=-,ch=getchar();
     for(;ch>=''&&ch<='';ch=getchar()) x=x*+ch-''; x*=f;
 }

 void pre() {
     For(i,,n) {
         For(j,,m) z[i][j]=mp[i][j]?:z[i][j-]+;
         Rep(j,m,) y[i][j]=mp[i][j]?:y[i][j+]+;
     }
     For(j,,m) {
         For(i,,n) s[i][j]=mp[i][j]?:s[i-][j]+;
         Rep(i,n,) x[i][j]=mp[i][j]?:x[i+][j]+;
     }
     For(i,,n) For(j,,m) {
         a[i][j]=min(y[i][j],x[i][j]);
         b[i][j]=min(z[i][j],s[i][j]);
     }
 }

 int sum[N],up;
 void add(int x,int v) {
     if(!x) return;
     for(int i=x;i<=up;i+=(i&(-i)))
         sum[i]+=v;
 }

 int qry(int x) {
     int rs=;
     for(int i=x;i;i-=(i&(-i)))
         rs+=sum[i];
     return rs;
 }

 vector<int>vc[N],vc2[N];
 void solve() {
     LL ans=;
     For(ty,,m) {
         int x=,y=ty;
         memset(sum,,sizeof(sum));
         For(i,,up) vc[i].clear(),vc2[i].clear();
         For(i,,up) {
             if(x+i>n||y+i>m) break;
             if(i-L+>=&&i-L+>i-b[x+i][y+i]) {
                 vc2[i-L+].push_back(i);
                 if(i-b[x+i][y+i]>=)
                     vc[i-b[x+i][y+i]].push_back(i);
             }
         }
         For(i,,up) {
             if(x+i>n||y+i>m) break;
             add(a[x+i][y+i]+i,);
             int upp=vc[i].size();
             For(j,,upp-)
                 ans-=qry(up)-qry(vc[i][j]);
             upp=vc2[i].size();
             For(j,,upp-)
                 ans+=qry(up)-qry(vc2[i][j]);
         }
     }
     For(tx,,n) {
         int x=tx,y=;
         memset(sum,,sizeof(sum));
         For(i,,up) vc[i].clear(),vc2[i].clear();
         For(i,,up) {
             if(x+i>n||y+i>m) break;
             if(i-L+>=&&i-L+>i-b[x+i][y+i]) {
                 vc2[i-L+].push_back(i);
                 if(i-b[x+i][y+i]>=)
                     vc[i-b[x+i][y+i]].push_back(i);
             }
         }
         For(i,,up) {
             if(x+i>n||y+i>m) break;
             add(a[x+i][y+i]+i,);
             int upp=vc[i].size();
             For(j,,upp-)
                 ans-=qry(up)-qry(vc[i][j]);
             upp=vc2[i].size();
             For(j,,upp-)
                 ans+=qry(up)-qry(vc2[i][j]);
         }
     }
     printf("%lld\n",ans);
 }

 #define ANS
 int main() {
 #ifdef ANS
     freopen("rampart.in","r",stdin);
     freopen("rampart.out","w",stdout);
 #endif
     read(n); read(m); read(L); read(tot);
     up=max(n,m);
     For(i,,tot) {
         int x,y;
         read(x); read(y);
         mp[x][y]=;
     }
     pre();
     solve();
     Formylove;
 }