第五周 6.14-6.20

6.14-6.15

什么都没干。

6.16

好久没码题了。快结课了会空很多。但在时间上有些难以抉择。

最主要的是后面还有一系列问题要困扰我。放假家里还有事情。

烦阿烦阿该何去何从。

校门外的树

因为是段增的。一开始想到了挖掘机那种方法。

但是段问的时候就挂了。

学了两颗树处理段增段问的方法。

 1 # include <iostream>
 2 # include <cstdio>
 3 # include <cstring>
 4 # include <algorithm>
 5 using namespace std;
 6 # define maxn 50000+5
 7 int n,c1[maxn],c2[maxn];
 8 
 9 int lowbit(int s)
10 {
11     return s&(-s);
12 }
13 
14 void add1(int i,int x)
15 {
16     while(i<=n) {c1[i]+=x; i+=lowbit(i);}
17     return;
18 }
19 
20 void add2(int i,int x)
21 {
22     while(i<=n) {c2[i]+=x; i+=lowbit(i);}
23     return;
24 }
25 
26 int sum1(int i)
27 {
28     int ans=0;
29     while(i>0) {ans+=c1[i]; i-=lowbit(i);}
30     return ans;
31 }
32 
33 int sum2(int i)
34 {
35     int ans=0;
36     while(i>0) {ans+=c2[i]; i-=lowbit(i);}
37     return ans;
38 }
39 
40 int main(void)
41 {
42     int m;
43     while(cin>>n>>m)
44     {
45         memset(c1,0,sizeof(c1));
46         memset(c2,0,sizeof(c2));
47         for(int i=0;i<m;i++)
48         {
49             int k,l,r;
50             scanf("%d%d%d",&k,&l,&r);
51             if(k==1) {add1(l,1); add2(r,1);}
52             else if(k==2) printf("%d
",sum1(r)-sum2(l-1));
53         }
54     }
55     return 0;
56 }

Aguin

6.17

SuperBrother打鼹鼠

二维树状数组

 1 # include <iostream>
 2 # include <cstdio>
 3 # include <cstring>
 4 using namespace std;
 5 # define maxn 1030
 6 int n,c[maxn][maxn];
 7 
 8 int lowbit(int s)
 9 {
10     return s&(-s);
11 }
12 
13 void add(int x,int y,int v)
14 {
15     for(int i=x;i<=n;i+=lowbit(i))
16         for(int j=y;j<=n;j+=lowbit(j))
17             c[i][j]+=v;
18     return;
19 }
20 
21 int sum(int x,int y)
22 {
23     int ans=0;
24     for(int i=x;i>0;i-=lowbit(i))
25         for(int j=y;j>0;j-=lowbit(j))
26             ans+=c[i][j];
27     return ans;
28 }
29 
30 int main(void)
31 {
32     while(cin>>n)
33     {
34         memset(c,0,sizeof(c));
35         int m;
36         while(cin>>m)
37         {
38             if(m==1)
39             {
40                 int x,y,k;
41                 scanf("%d%d%d",&x,&y,&k);
42                 add(x+1,y+1,k);
43             }
44             else if(m==2)
45             {
46                 int x1,y1,x2,y2;
47                 scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
48                 printf("%d
",sum(x2+1,y2+1)-sum(x1,y2+1)-sum(x2+1,y1)+sum(x1,y1));    
49             }
50             else if(m==3) break;
51         }
52     }
53     return 0;
54 }

Aguin

6.18

P1512 SuperBrother打鼹鼠

依旧此题。昨晚就在写。刷牙的时候感觉不对。今天推了重写。

还是去vijos交。自家oj不知怎么一直爆M阿。

二维线段树。点更新。面询问。

  1 # include <iostream>
  2 # include <cstdio>
  3 using namespace std;
  4 # define maxn 1030
  5 int n;
  6 
  7 struct nodey
  8 {
  9     int yl,yr,sum;
 10 };
 11 
 12 struct nodex
 13 {
 14     int xl,xr;
 15     nodey treey[maxn*4];
 16     
 17     void buildtree(int i,int l,int r)
 18     {
 19         treey[i].yl=l;treey[i].yr=r;treey[i].sum=0;
 20         if(l<r)
 21         {
 22             buildtree(2*i,l,(l+r)/2);
 23             buildtree(2*i+1,(l+r)/2+1,r);
 24         }
 25     //    printf("build x:[%d,%d]y:[%d,%d]sum:%d
",xl,xr,treey[i].yl,treey[i].yr,treey[i].sum);
 26         return;
 27     }
 28     
 29     void update(int i,int y,int v)
 30     {
 31         treey[i].sum+=v;
 32     //    printf("update x:[%d,%d]y:[%d,%d]sum:%d
",xl,xr,treey[i].yl,treey[i].yr,treey[i].sum);
 33         if(treey[i].yl==treey[i].yr) return;
 34         else if(y<=(treey[i].yl+treey[i].yr)/2) update(2*i,y,v);
 35         else update(2*i+1,y,v);
 36         return;
 37     }
 38     
 39     int query(int i,int l,int r)
 40     {
 41         int ans=0;
 42         if(treey[i].yl>=l&&treey[i].yr<=r) return treey[i].sum;
 43         if(l<=(treey[i].yl+treey[i].yr)/2) ans+=query(2*i,l,r);
 44         if(r>=(treey[i].yl+treey[i].yr)/2+1) ans+=query(2*i+1,l,r);
 45         return ans;
 46     }
 47     
 48 }treex[maxn*4];
 49 
 50 void buildtreex(int i,int l,int r)
 51 {
 52     treex[i].xl=l;treex[i].xr=r;
 53     treex[i].buildtree(1,1,n);
 54     if(l<r)
 55     {
 56         buildtreex(2*i,l,(l+r)/2);
 57         buildtreex(2*i+1,(l+r)/2+1,r);
 58     }
 59     return;
 60 }
 61 
 62 void updatex(int i,int x,int y,int v)
 63 {
 64     treex[i].update(1,y,v);
 65     if(treex[i].xl==treex[i].xr) return;
 66     if(x<=(treex[i].xl+treex[i].xr)/2) updatex(2*i,x,y,v);
 67     else updatex(2*i+1,x,y,v);
 68     return;
 69 }
 70 
 71 int queryx(int i,int x1,int y1,int x2,int y2)
 72 {
 73     int ans=0;
 74     if(treex[i].xl>=x1&&treex[i].xr<=x2) return treex[i].query(1,y1,y2);
 75     if(x1<=(treex[i].xl+treex[i].xr)/2) ans+=queryx(2*i,x1,y1,x2,y2);
 76     if(x2>=(treex[i].xl+treex[i].xr)/2+1) ans+=queryx(2*i+1,x1,y1,x2,y2);
 77     return ans;
 78 }
 79 
 80 int main(void)
 81 {
 82     while((scanf("%d",&n))!=EOF)
 83     {
 84         buildtreex(1,1,n);
 85         int m;
 86         while(1)
 87         {
 88             scanf("%d",&m);
 89             if(m==3) break;
 90             if(m==1)
 91             {
 92                 int x,y,k;
 93                 scanf("%d%d%d",&x,&y,&k);
 94                 updatex(1,x+1,y+1,k);
 95             }
 96             else if(m==2)
 97             {
 98                 int x1,y1,x2,y2;
 99                 scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
100                 printf("%d
",queryx(1,x1+1,y1+1,x2+1,y2+1));
101             }
102         }
103     }
104     return 0;
105 }

Aguin

update可以写成非递归形式

 1 # include <iostream>
 2 # include <cstdio>
 3 using namespace std;
 4 # define maxn 1030
 5 int n,locx[maxn],locy[maxn];
 6 
 7 struct nodey
 8 {
 9     int yl,yr,sum;
10 };
11 
12 struct nodex
13 {
14     int xl,xr;
15     nodey treey[maxn*4];
16     
17     void buildtree(int i,int l,int r)
18     {
19         treey[i].yl=l;treey[i].yr=r;treey[i].sum=0;
20         if(l<r)
21         {
22             buildtree(2*i,l,(l+r)/2);
23             buildtree(2*i+1,(l+r)/2+1,r);
24         }
25         else locy[l]=i;
26     //    printf("build x:[%d,%d]y:[%d,%d]sum:%d
",xl,xr,treey[i].yl,treey[i].yr,treey[i].sum);
27         return;
28     }
29     
30     int query(int i,int l,int r)
31     {
32         int ans=0;
33         if(treey[i].yl>=l&&treey[i].yr<=r) return treey[i].sum;
34         if(l<=(treey[i].yl+treey[i].yr)/2) ans+=query(2*i,l,r);
35         if(r>=(treey[i].yl+treey[i].yr)/2+1) ans+=query(2*i+1,l,r);
36         return ans;
37     }
38     
39 }treex[maxn*4];
40 
41 void buildtreex(int i,int l,int r)
42 {
43     treex[i].xl=l;treex[i].xr=r;
44     treex[i].buildtree(1,1,n);
45     if(l<r)
46     {
47         buildtreex(2*i,l,(l+r)/2);
48         buildtreex(2*i+1,(l+r)/2+1,r);
49     }
50     else locx[l]=i;
51     return;
52 }
53 
54 void UPDATE(int x,int y,int v)
55 {
56     for(int i=locx[x];i;i/=2)
57         for(int j=locy[y];j;j/=2)
58             treex[i].treey[j].sum+=v;
59     return;
60 }
61 
62 int queryx(int i,int x1,int y1,int x2,int y2)
63 {
64     int ans=0;
65     if(treex[i].xl>=x1&&treex[i].xr<=x2) return treex[i].query(1,y1,y2);
66     if(x1<=(treex[i].xl+treex[i].xr)/2) ans+=queryx(2*i,x1,y1,x2,y2);
67     if(x2>=(treex[i].xl+treex[i].xr)/2+1) ans+=queryx(2*i+1,x1,y1,x2,y2);
68     return ans;
69 }
70 
71 int main(void)
72 {
73     while((scanf("%d",&n))!=EOF)
74     {
75         buildtreex(1,1,n);
76         int m;
77         while(1)
78         {
79             scanf("%d",&m);
80             if(m==3) break;
81             if(m==1)
82             {
83                 int x,y,k;
84                 scanf("%d%d%d",&x,&y,&k);
85                 UPDATE(x+1,y+1,k);
86             }
87             else if(m==2)
88             {
89                 int x1,y1,x2,y2;
90                 scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
91                 printf("%d
",queryx(1,x1+1,y1+1,x2+1,y2+1));
92             }
93         }
94     }
95     return 0;
96 }

Aguin

在自家oj。树状数组好快跑完。

然而二维线段树爆M。不明所以阿。

vijos上树状数组。二维线段树。都差不多。

6.19

听了郏老大的平衡树课。然而觉得这些离我还太遥远。

作为一只小菜比。基础方面还有很大空缺。

影响最大的就是不会dp。

截取一段不知道说什么但是看起来很牛掰的网文

修炼dp。从数字三角形开始。

数字三角形

 1 # include <iostream>
 2 # include <cstdio>
 3 # include <algorithm>
 4 # include <cstring>
 5 using namespace std;
 6 int N,num[26][26],dp[26][26];
 7 
 8 int solve(int i,int j)
 9 {
10     if(dp[i][j]>=0) return dp[i][j];
11     return dp[i][j]=num[i][j]+( (i==N) ? 0 : max(solve(i+1,j),solve(i+1,j+1)) );
12 }
13 
14 int main(void)
15 {
16     cin>>N;
17     memset(dp,-1,sizeof(dp));
18     for(int i=1;i<=N;i++)
19         for(int j=1;j<=i;j++)
20             scanf("%d",&num[i][j]);
21     cout<<solve(1,1)<<endl;
22     return 0;
23 }

Aguin

6.20

打了个bc

hdu5273 Dylans loves sequence

比赛的时候也是发现N小想用蠢蠢的办法结果没有调完就结束了。

赛后看题解。思路基本上和提到的第一种做法一样。

（hack的时候也有偷看别人代码有看到树状数组的但大部分都是用简单方法）

赛后写完过了。然而觉得中间写的很复杂。离散化直接用了map。

 1 # include <iostream>
 2 # include <cstdio>
 3 # include <cstring>
 4 # include <map>
 5 # include <algorithm>
 6 using namespace std;
 7 int N,a[1005],b[1005],c[1005];
 8 
 9 struct node
10 {
11     int no,L,R,ans;
12 }query[100005];
13 
14 bool cmp1(node a,node b)
15 {
16     if(a.L!=b.L) return a.L<b.L;
17     else return a.R<b.R;
18 }
19 
20 bool cmp2(node a,node b)
21 {
22     return a.no<b.no;
23 }
24 
25 int lowbit(int s)
26 {return s&(-s);}
27 
28 void add(int i,int x)
29 {
30     while(i<=N){c[i]+=x; i+=lowbit(i);}
31     return;
32 }
33 
34 int sum(int i)
35 {
36     int ans=0;
37     while(i>0) {ans+=c[i];i-=lowbit(i);}
38     return ans;
39 }
40 
41 int main(void)
42 {
43     int Q;cin>>N>>Q;
44     for(int i=1;i<=N;i++) scanf("%d",a+i),b[i]=a[i];
45     sort(b+1,b+1+N);
46     map<int,int> m;
47     for(int i=1;i<=N;i++) m[b[i]]=i;
48     for(int i=1;i<=Q;i++)
49     {
50         scanf("%d%d",&query[i].L,&query[i].R);
51         query[i].no=i; query[i].ans=0;
52     }
53     sort(query+1,query+1+Q,cmp1);    
54     for(int i=1;i<=Q;i++)
55     {
56         memset(c,0,sizeof(c));
57         for(int j=query[i].L;j<=query[i].R;j++) {query[i].ans+=sum(N)-sum(m[a[j]]);add(m[a[j]],1);}
58         for(int j=i+1;j<=Q;j++)
59         {
60             if(query[j].L==query[j-1].L)
61             {
62                 query[j].ans=query[j-1].ans;
63                 for(int k=query[j-1].R+1;k<=query[j].R;k++)
64                 {
65                     query[j].ans+=sum(N)-sum(m[a[k]]);
66                     add(m[a[k]],1);
67                 }
68                 i=j;
69             }
70             else {i=j-1;break;}
71         }
72     }
73     sort(query+1,query+1+Q,cmp2);
74     for(int i=1;i<=Q;i++) printf("%d
",query[i].ans);
75     return 0;
76 }

Aguin

相关阅读:
linux命令学习笔记：cut详解
 浏览器url传参中文时得到null的解决方法
 jQuery给控件赋值....
Myeclipse 错误An internal error has occurred 解决办法
 浏览器发送URL的编码特性
 DEBUG -- CLOSE BY CLIENT STACK TRACE问题的两种解决方案，整理自网络
 解决Maven中OutOfMemory错误
 java中的URLEncoder和URLDecoder类;中文在地址栏中的处理
 关于CLOSE BY CLIENT STACK TRACE
Hibernate复合主键映射
原文地址：https://www.cnblogs.com/Aguin/p/4575922.html