首页 > 代码库 > 数三角形 bzoj 1201

数三角形 bzoj 1201

数三角形(1s 128MB)triangle

【题目描述】

小苏看到一个这样的等边三角形:该等边三角形每边的长度为n且被分成n等份,于是每条边就有n-1个等分点。而整个三角形被连接两个不同边的等分点且平行于三角形的第三边的线段分成了n2个单位等边三角形(边长为1)。下图左是n=5的情形:

小苏想知道,删除其中的一些短边后,剩下的边一共组成多少个三角形(包括所有边长为m的三角形),正立的和倒立的都算,只要三角形的3m条短边都没有被删除就算是组成一个三角形)。例如,上图右就存在19个三角形。

【输入格式】

大三角形的所有短边可以看成由(n+1)*n/2个单位三角形的边界组成。如下图的灰色三角形所示。其中第1排有1个灰色三角形,第2排有2个灰色三角形,……,第n排有n个灰色三角形。所以输入格式是这样规定的:

从文件中读入数据,文件第一行为正整数n,表示大三角形每边的长度。接下来的n行,第i+1行有i组数,从左到右每组数描述一个三角形,每组数都有3个数,这3个数非0即1,表示对应的短边是否被删除,0表示已被删除,1表示未被删除,依次按照三角形的左、右、下边的顺序来描述。所以第i+1行有3i个数,每个数是0或1。

【输出格式】

输出文件中仅包含一个整数T,表示有多少个三角形的边界都没有被删除。

【输入样例】

5

1 1 1

1 1 0 1 1 0

1 1 1 1 1 1 1 0 1

1 0 1 1 1 1 0 1 1 1 1 1

0 1 1 1 1 1 0 1 1 1 1 1 0 1 1

【输出样例】

19

【数据范围】

1<=n<=1000,1<=m<=n。

 


题解:

主要算法:树状数组

先考虑正三角,枚举底边所在的直线,一个合法三角形的底边是实线,左右两边比底边长

预处理:

c[i].l : 往左上延伸的最长长度

c[i].r : 往右上延伸的最长长度

r[i] : 向右延伸的最长长度

设一条底边的两个编号为 i 和 j

那么一个合法三角形需要满足:

(1)点 i 在点 j 的左边 : i < j -> i >= j

(2)i 到 j 之间为实边: r[i] >= j

(3)左边长于底边 : j - i <= c[i].r -> c[i].r + i >= j

(4)右边长于底边 : j - i <= c[j].l -> i >= j + c[j].l

那么转化为了关于 j 的不等式组

然后将 j 按 j - c[j].l 排序

从大到小枚举 i ,在树状数组中加入满足(4)的点 j

那么答案就是满足 (2) 和 (3) 的元素个数中的较小值减去转化过的条件(1)的元素个数

其实暴力能过······

  1 #include<algorithm>
  2 #include<iostream>
  3 #include<cstring>
  4 #include<cstdlib>
  5 #include<cstdio>
  6 #include<cmath>
  7 using namespace std;
  8 inline int Get()
  9 {
 10     int x = 0;
 11     char c = getchar();
 12     while(0 > c || c > 9) c = getchar();
 13     while(0 <= c && c <= 9)
 14     {
 15         x = (x << 3) + (x << 1) + c - 0;
 16         c = getchar();
 17     }
 18     return x;
 19 }
 20 struct shape
 21 {
 22     int l, r, m;
 23 };
 24 struct mystery
 25 {
 26     int v, i;
 27 };
 28 shape a[1233][1233], c[1233][1233];
 29 mystery s[1000233];
 30 int n;
 31 int ans;
 32 int r[1233][1233];
 33 int tr[1000233];
 34 inline void Add(int x, int y)
 35 {
 36     while(x <= n + 1)
 37     {
 38         tr[x] += y;
 39         x += x & (-x);
 40     }
 41 }
 42 inline int Sum(int x)
 43 {
 44     int sum = 0;
 45     while(x)
 46     {
 47         sum += tr[x];
 48         x -= x & (-x);
 49     }
 50     return sum;
 51 }
 52 inline bool rule(mystery a, mystery b)
 53 {
 54     if(a.v != b.v) return a.v < b.v;
 55     return a.i < b.i;
 56 }
 57 inline void Up()
 58 {
 59     for(int i = 1; i <= n; ++i)
 60     {
 61         for(int j = 1; j <= i + 1; ++j)
 62         {
 63             if(!a[i][j].l) c[i][j].r = 0;
 64             else c[i][j].r = c[i - 1][j].r + 1;
 65             if(!a[i][j - 1].r) c[i][j].l = 0;
 66             else c[i][j].l = c[i - 1][j - 1].l + 1;
 67         }
 68         for(int j = i + 1; j >= 1; --j)
 69         {
 70             if(!a[i][j].m) r[i][j] = j;
 71             else r[i][j] = r[i][j + 1];
 72         }
 73     }
 74     for(int k = 1; k <= n; ++k)
 75     {
 76         for(int i = 1; i <= k + 1; ++i)
 77         {
 78             s[i].v = i - c[k][i].l;
 79             s[i].i = i;
 80         }
 81         sort(s + 1, s + 1 + k + 1, rule);
 82         int num = 1;
 83         for(int i = 1; i <= k + 1; ++i)
 84         {
 85             while(num <= k + 1 && s[num].v <= i)
 86             {
 87                 Add(s[num].i, 1);
 88                 ++num;
 89             }
 90             if(c[k][i].r && r[k][i])
 91             {
 92                 int one = Sum(c[k][i].r + i);
 93                 int two = Sum(i);
 94                 int thr = Sum(r[k][i]);
 95                 ans += min(one, thr) - two;
 96             }
 97         }
 98         for(int i = 1; i < num; ++i) Add(s[i].i, -1);
 99     }
100 }
101 inline void Down()
102 {
103     for(int i = n; i >= 1; --i)
104     {
105         for(int j = 1; j <= i; ++j)
106         {
107             if(!a[i][j].r) c[i][j].r = 0;
108             else c[i][j].r = c[i + 1][j + 1].r + 1;
109             if(!a[i][j].l) c[i][j].l = 0;
110             else c[i][j].l = c[i + 1][j].l + 1;
111         }
112         for(int j = i; j >= 1; --j)
113         {
114             if(!a[i - 1][j].m) r[i][j] = j;
115             else r[i][j] = r[i][j + 1];
116         }
117     }
118     for(int k = 1; k <= n; ++k)
119     {
120         for(int i = 1; i <= k; ++i)
121         {
122             s[i].v = i - c[k][i].l;
123             s[i].i = i;
124         }
125         sort(s + 1, s + 1 + k, rule);
126         int num = 1;
127         for(int i = 1; i <= k; ++i)
128         {
129             while(num <= k && s[num].v <= i)
130             {
131                 Add(s[num].i, 1);
132                 ++num;
133             }
134             
135             if(c[k][i].r && r[k][i])
136             {
137                 int one = Sum(c[k][i].r + i);
138                 int two = Sum(i);
139                 int thr = Sum(r[k][i]);
140                 ans += min(one, thr) - two;
141             }
142         }
143         for(int i = 1; i < num; ++i) Add(s[i].i, -1);
144     }
145 }
146 int main()
147 {
148     n = Get();
149     for(int i = 1; i <= n; ++i)
150         for(int j = 1; j <= i; ++j)
151         {
152             a[i][j].l = Get();
153             a[i][j].r = Get();
154             a[i][j].m = Get();
155         }
156     Up();
157     Down();
158     printf("%d", ans);
159 }

 

数三角形 bzoj 1201