算法学社
記錄難忘的征途
posts - 141,comments - 220,trackbacks - 0
1. 涨了25pt...

A

给500,000个字符串,如果两个字符串的首字母和尾字母相连就说明这两个字符串相连。请问可以组成的最长环是多少...
分析:
一开始我看这不是求最长欧拉回路么...
后来发现一个条件: 只能前面的串去连后面的串,然后最后的串与第一个串相连... 但是当时没有想清楚...

B

在一个长度小于100,000的数列{S}中,挑选K(K<100,000)个数,让这K个数大于B且最左端的数Si的i值最小....
分析:
从右到左扫一遍... 每次假设选中了Si,且维护一下右边数列第k-1大的值,用优先级队列搞。
代码
 1 #include<cstdio>
 2 #include<algorithm>
 3 #include<queue>
 4 #include<iostream>
 5 using namespace std;
 6 typedef long long ll;
 7 int num[100005];
 8 priority_queue <int,vector<int>,greater<int>  > Q;
 9 int main(){
10     ll b;int n,k;
11     while(cin >> n >> k){
12         cin >>b;
13         for(int i=0;i<n;i++)
14             cin >>num[i];
15         ll sum = 0;
16         int ans = n;
17         n--;
18         while(!Q.empty())Q.pop();
19         for(int i = n-k+1; i<n; i++){
20             if(i < 0) continue;
21             sum += num[i];
22             Q.push(num[i]);
23         }
24         for(int i = n-k;i>=0; i--){
25 //            cout<<sum<<endl;
26             sum+= num[i];
27             Q.push(num[i]);
28 //            cout<<i<<" "<<sum<<endl;
29 //            cout<<Q.top()<<endl;
30             if(sum > b) ans = i+1;
31             if(sum-Q.top() > b) {
32                 ans = 1;
33                 break;
34             }
35             sum -= Q.top();
36             Q.pop();
37         }
38         cout<<ans<<endl;
39     }
40 }
41 

C

给一个点数为100,000的树,每次询问可以将u和v之间的路的每条边权加1。最后请输出所有边的边权。
分析:
动态树和树链剖分都可,其实还有更简单的方法。不过没时间想了... 直接甩模板了
树链剖分版本:
  1 // template
  2 #include<iostream>
  3 #include<algorithm>
  4 #include<cassert>
  5 #include<cstdio>
  6 #include<cstdlib>
  7 #include<cstring>
  8 using namespace std;
  9 template <typename T> inline void chkmax(T &a, T b){if(a<b) a=b;}
 10 // graph
 11 const int V = 100005;
 12 const int E = 200005;
 13 int head[V],pnt[E],nxt[E],flag[E];
 14 int n,e;
 15 void addedge(int u,int v){
 16     nxt[e] = head[u];
 17     pnt[e] = v;
 18     head[u] = e;
 19     e ++;
 20 }
 21 // dsu 
 22 int parent[V];
 23 int find(int x){ return x == parent[x] ? x : parent[x] = find(parent[x]);}
 24 // seg_ment tree
 25 int seg[V<<2], M;
 26 void find(int l,int r){
 27     //cout<<l<<" "<<r<<endl;
 28     for(l += M-1, r += M+1; l^r^1; l>>=1, r>>=1){
 29         if(~l&1) seg[l^1]++;
 30         if(r&1) seg[r^1]++;
 31     }
 32 }
 33 // prepare
 34 int deep[V],size[V],heavy[V],P[V];
 35 void dfs(int u,int f){
 36     size[u] = 1;
 37     int mx = 0, s = -1;
 38     for(int i=head[u]; i!=-1;i = nxt[i]){
 39         if( pnt[i] == f) continue;
 40         int v = pnt[i];
 41         P[v] = i^1;
 42         deep[v] = deep[u] + 1;
 43         dfs(v,u);
 44         if(size[v] > mx){
 45             mx = size[v];
 46             s = i;
 47         }
 48         size[u] += size[v];
 49     }
 50     heavy[u] = s;
 51     if(s!=-1) parent[pnt[s]] = u;
 52 }
 53 void prepare(){
 54     for(int i=0;i<n;i++) parent[i] = i;
 55     deep[0] = 0;
 56     P[0] = -1;
 57     dfs(0,0);
 58     for(int i=30;i;i--) if((1<<i) > n+1) M = 1<<i;
 59     for(int i=0;i<2*M;i++) seg[i] = 0;
 60     int len = 1;
 61     for(int u = 0; u<n; u++) if(heavy[u] == -1){
 62         int v = u;
 63         while(v && pnt[heavy[pnt[P[v]]]] == v){
 64             flag[P[v]] = flag[P[v]^1] = len ++;
 65             v = pnt[P[v]];
 66         }
 67     }
 68 }
 69 // operator
 70 int lca(int u,int v){
 71     while(1){
 72         int a = find(u), b = find(v);
 73         if(a == b) return deep[u]<deep[v] ? u : v;
 74         else if(deep[a] > deep[b]) u = pnt[P[a]];
 75         else v = pnt[P[b]];
 76     }
 77 }
 78 int cnt[E];
 79 void query(int u,int v){
 80     //cout<<"query: "<<u<<" "<<v<<endl; 
 81     while(u != v){
 82         int l = P[u];
 83         if(pnt[heavy[pnt[P[u]]]] == u){
 84             int p = find(u);
 85             if(deep[p] < deep[v]) p = v;
 86             int r = heavy[p];
 87             assert(flag[l] <= flag[r]);
 88             find(flag[l],flag[r]);
 89             u = p;
 90         }
 91         else {
 92             u = pnt[l];
 93             cnt[l]++; cnt[l^1]++;
 94 //            cout<<l<<endl;
 95         }
 96     }
 97 }
 98 void ask(int a,int b){
 99     int p = lca(a,b);
100 //    cout<<a<<" "<<b<<" "<<p<<endl;
101     query(a,p); query(b,p);
102 }
103 // main
104 int main(){
105     int k;
106     while(~scanf("%d",&n)){
107         e = 0;
108         memset(head,-1,sizeof(head));
109         memset(cnt,0,sizeof(cnt));
110         for(int i=0;i<n-1;i++){
111             int u,v;
112             scanf("%d%d",&u,&v);
113             u--; v--;
114             addedge(u,v);
115             addedge(v,u);
116         }
117         prepare();
118         cin >>k;
119         while(k--){
120             int u,v;
121             scanf("%d%d",&u,&v);
122             u--;v--;
123             ask(u,v);
124         }
125         for(int i=0;i<e;i+=2){
126             int v = pnt[i], u = pnt[i^1];
127             if(heavy[v]!=(i^1) && heavy[u]!=i){
128                 cout<<cnt[i]<<" ";
129             }
130             else {
131                 int pos = flag[i]+M,ans = 0;
132                 while(pos) {ans += seg[pos]; pos >>=1;}
133                 cout<<ans<<" ";
134             }
135         }
136         cout<<endl;
137     }
138 }
139 
posted on 2012-05-28 09:07 西月弦 阅读(755) 评论(8)  编辑 收藏 引用 所属分类: 比赛感言

FeedBack:
# re: codeforces #121 div1 [未登录]
2012-05-28 23:09 | David
你好,请问有没有详解动态树或树链剖分的论文?  回复  更多评论
  
# re: codeforces #121 div1
2012-05-29 09:31 | 西月弦
http://fanhq666.blog.163.com/blog/static/819434262011518104215977/
动态树
树链剖分看这篇论文
漆子超《分治算法在树的路径问题中的应用》
@David  回复  更多评论
  
# re: codeforces #121 div1 [未登录]
2012-05-30 00:10 | David
@西月弦
非常感谢!  回复  更多评论
  
# re: codeforces #121 div1
2012-06-27 13:05 | 蒟蒻
请问C动态树怎么整...?大牛能否给个思路...link cut tree刚学,不会维护区间和之类的...  回复  更多评论
  
# re: codeforces #121 div1
2012-06-27 14:43 | 西月弦
@蒟蒻
不好意思,我不会LCT...  回复  更多评论
  
# re: codeforces #121 div1
2012-06-27 19:14 | 蒟蒻
@西月弦
谢谢...大牛的树链剖分的模版在哪找的...自己写的吗?

  回复  更多评论
  
# re: codeforces #121 div1
2012-06-27 19:36 | 西月弦
@蒟蒻
恩,自己写的,比较丑额。。。。  回复  更多评论
  
# re: codeforces #121 div1
2012-06-28 12:39 | 蒟蒻
@西月弦
谢谢...  回复  更多评论
  

只有注册用户登录后才能发表评论。
网站导航: 博客园   IT新闻   BlogJava   知识库   博问   管理