C++博客-LSM ACM/ICPC算法程序设计空间

marker切奇数题，我切偶数题

LSM — Mon, 17 Sep 2007 12:17:00 GMT

1、http://acm.pku.edu.cn/JudgeOnline/problem?id=1474
2、http://acm.pku.edu.cn/JudgeOnline/problem?id=2794
3、http://acm.pku.edu.cn/JudgeOnline/problem?id=1366
4、http://acm.pku.edu.cn/JudgeOnline/problem?id=1392
5、http://acm.pku.edu.cn/JudgeOnline/problem?id=2169

6、http://acm.pku.edu.cn/JudgeOnline/problem?id=3150
7、http://acm.pku.edu.cn/JudgeOnline/problem?id=1400
8、http://acm.pku.edu.cn/JudgeOnline/problem?id=2065
9、http://acm.pku.edu.cn/JudgeOnline/problem?id=1860
10、http://acm.pku.edu.cn/JudgeOnline/problem?id=2168

11、http://acm.pku.edu.cn/JudgeOnline/problem?id=2791
12、http://acm.pku.edu.cn/JudgeOnline/problem?id=2532
13、http://acm.pku.edu.cn/JudgeOnline/problem?id=1712
14、http://acm.pku.edu.cn/JudgeOnline/problem?id=1709
15、http://acm.pku.edu.cn/JudgeOnline/problem?id=1457

16、http://acm.pku.edu.cn/JudgeOnline/problem?id=1134
17、http://acm.pku.edu.cn/JudgeOnline/problem?id=1137
18、http://acm.pku.edu.cn/JudgeOnline/problem?id=1481

LSM 2007-09-17 20:17 发表评论

计算几何相关模板（更新中……）

LSM — Wed, 22 Aug 2007 10:39:00 GMT

计算几何相关模板（更新中……）

最近在学计算几何，边学，边整理模板，有错的话请大家指出！

1//计算几何模板 ~ alpc02
2const double PRECISION = 1e-8;
3struct Point {
4    double x, y;
5};
6int dblcmp(double d) {
7    return (fabs(d) < PRECISION) ? 0:(d>0 ? 1:-1);
8} //三叉口函数，避免精度误差
9double length(double x, double y) {
10    return sqrt(x*x + y*y);
11} //向量长度
12double dotdet(double x1, double y1, double x2, double y2) {
13    return x1*x2 + y1*y2;
14} //点积
15double det(double x1, double y1, double x2, double y2) {
16    return x1*y2 - x2*y1;
17} //叉积
18int cross(const Point &a, const Point &c, const Point &d) {
19    return dblcmp( det(a.x-c.x, a.y-c.y, d.x-c.x, d.y-c.y) );
20} //右手螺旋定则，1——a在cd右侧，-1——a在cd左侧，0——三点共线
21bool between(const Point &a, const Point &c, const Point &d) {
22    return dblcmp( dotdet(c.x-a.x, c.y-a.y, d.x-a.x, d.y-a.y) ) != 1;
23} //在cross(a,c,d)==0的基础上，可判断点a是否在cd内部
24int segIntersect(const Point &a, const Point &b, const Point &c, const Point &d) {
25    int a_cd = cross(a,c,d);
26    if(a_cd == 0 && between(a,c,d))    return 2;
27    int b_cd = cross(b,c,d);
28    if(b_cd == 0 && between(b,c,d))    return 2;
29    int c_ab = cross(c,a,b);
30    if(c_ab == 0 && between(c,a,b))    return 2;
31    int d_ab = cross(d,a,b);
32    if(d_ab == 0 && between(d,a,b))    return 2;
33    if ((a_cd ^ b_cd) == -2 && (c_ab ^ d_ab) == -2)
34        return 1;
35    return 0;
36} //两线段相交情况：0——不相交，1——规范相交，2——不规范相交（交于端点或重合）
37void intersectPoint(const Point &a, const Point &b, const Point &c, const Point &d, Point &e) {
38    double sc, sd;
39    sc = fabs( det(b.x-a.x, b.y-a.y, c.x-a.x, c.y-a.y) );
40    sd = fabs( det(b.x-a.x, b.y-a.y, d.x-a.x, d.y-a.y) );
41    e.x = (sc * d.x + sd * c.x) / (sc + sd);
42    e.y = (sc * d.y + sd * c.y) / (sc + sd);
43} //两线段规范相交时，求交点坐标
44int linesegIntersect(const Point &a, const Point &b, const Point &c, const Point &d) {
45    int c_ab = cross(c,a,b);
46    if(c_ab == 0)    return 2;
47    int d_ab = cross(d,a,b);
48    if(d_ab == 0)    return 2;
49    if(c_ab ^ d_ab == -2)
50        return 1;
51    return 0;
52} //直线ab和线段cd相交情况：0——不相交，1——规范相交，2——不规范相交（交于端点或重合）
53int lineIntersect(const Point &a, const Point &b, const Point &c, const Point &d) {
54    if(dblcmp(det(b.x-a.x, b.y-a.y, d.x-c.x, d.y-c.y)) != 0)
55        return 1;
56    if(cross(a,c,d) == 0)
57        return 2;
58    return 0;
59} //两直线相交情况：0——平行，1——规范相交，2——不规范相交（重合）
60

LSM 2007-08-22 18:39 发表评论

计算几何入门题（更新中……）

LSM — Wed, 22 Aug 2007 10:28:00 GMT

摘要: 计算几何入门题（更新中……） http://acm.pku.edu.cn/JudgeOnline/problem?id=1696Space Ant用到叉积，点积。这个题的代码稍微改改，貌似可以认为是凸包的o(n^2)算法，即卷包心菜。 Code highlighting produced by Actipro CodeHighlighter (freeware)http... 阅读全文

LSM 2007-08-22 18:28 发表评论

POJ 3351 —— Gerrymandering

LSM — Wed, 22 Aug 2007 03:30:00 GMT

http://acm.pku.edu.cn/JudgeOnline/problem?id=3351
Gerrymandering
难！O(n^3)的dp，降为接近O(n^2)才能过。

不难想到O(n^3)的dp，dp[k][t]代表前k个riding合并成t块（也就是说合并k-t次），所能取得的最大值（赢的riding个数减去输的riding个数）。定义w[i][j]代表i至j合并一起后的输赢状态，1或-1。
1、dp[0][0] = 0，状态转移方程dp[k][t] = max(dp[k'][t-1] + win[k'+1][k]), 0<=k'2、目标则是max(t) 使得dp[n][t] > 0，输出答案n-t。
转移的复杂度是o(n)，这样我们的复杂度就是O(n^3)，必须降低复杂度。

想法：我们求dp[k][t]时，遍历了一次dp[k'][t-1]，当我们再求dp[k+1][t]时，还需要再遍历一次dp[k'][t-1]，只不过这个时候k'范围增加了1。那我们能否在求dp[k][t]遍历dp[k'][t-1]时把这些有用信息记录下来供求dp[k+1][t]时用呢？

看状态转移方程：dp[k][t] = max(dp[k'][t-1] + win[k'+1][k])
问题：受win[k'+1][k]影响，dp[k'][t-1] 的最大值，加上win[k'+1][k]后并不一定是dp[k][t]的最大值
但是：仔细看会发现，win[k'+1][k]仅有两种取值，1和-1,而且当k'变化时dp[k'][t-1]的值相差2、4、6
这样的话，方程可以写成：dp[k][t] = maxdp + win[k''+1][k]，这里maxdp = max(dp[k'][t-1])，k''有一个限制，就是，dp[k''][t-1]必须等于maxdp。
一旦转换成这类方程后，dp[k'][t-1]的信息就可以重复利用，状态转移的复杂度大大降低。

实现：记录maxdp，然后建一个list，记录满足dp[k''][t-1] = maxdp的k'’，实时更新list，每次只需从list里转移状态就行。
优化：一旦从list里找到一个k''，满足win[k''+1][k]=1，那么dp[k][t] = maxdp + 1，dp[k][t]的计算就可马上结束。

1 //POJ_3351 ～ alpc02
2 #include <algorithm>
3 using namespace std;
4
5 const int INF = 12345678;
6 const int N = 1010;
7 const int P = 12;
8
9 int n, p;
10 int v[N][P];
11 int win[N][N];
12 int dp[N][N];
13
14 int Max(int a,int b){return a>b ? a:b;}
15 void input() {
16     scanf("%d", &n);
17     scanf("%d", &p);
18     int i, j, k;
19     memset(v[0], 0, sizeof(v[0]));
20     for(i=1; i<=n; i++)
21         for(j=1; j<=p; j++) {
22             scanf("%d", &k);
23             v[i][j] = v[i-1][j] + k;
24         }
25 }
26 void init() {
27     int i, j, k, vote;
28     for(i=1; i<=n; i++)
29         for(j=i; j<=n; j++) {
30             vote = v[j][1] - v[i-1][1];
31             for(k=2; k<=p; k++)
32                 if(v[j][k] - v[i-1][k] >= vote)
33                     break;
34             if(k > p)    win[i][j] = 1;
35             else        win[i][j] = -1;
36         }
37 }
38 int solve() {
39     init();
40
41     int i, j, k, t;
42     int list[N], top, mx;
43     dp[0][0] = 0;
44     for(j=1; j<=n; j++) {
45         top = 0;
46         mx = -INF;
47         for(i=j; i<=n; i++) {
48             if(dp[i-1][j-1] > mx) {
49                 mx = dp[i-1][j-1];
50                 list[0] = i-1;
51                 top = 1;
52             }
53             else if(dp[i-1][j-1] == mx)
54                 list[top++] = i-1;
55             for(t=0; t<top-1; t++)
56                 if(win[list[t]+1][i] == 1)
57                     break;
58             dp[i][j] = mx + win[list[t]+1][i];
59         }
60     }
61
62     int ans;
63     for(ans=n; ans>=1; ans--)
64         if(dp[n][ans] > 0)
65             return n-ans;
66     return -1;
67 }
68 int main() {
69     input();
70     printf("%d\n", solve());
71     return 0;
72 }

LSM 2007-08-22 11:30 发表评论

计算几何初步认识(一）

LSM — Tue, 21 Aug 2007 00:52:00 GMT

计算几何初步认识

——Wjx昨晚给我速成的计算几何

一、点。

l 点的坐标A(x1, y1)，B(x2, y2)

二、向量。

l 向量AB = (x2-x1, y2-y1) = (x3,y3) ，CD = (x4, y4)。

l 向量的模|AB| = sqrt(x3*x3 + y3*y3) 即向量的长度。

三、点积。

l 点积的结果为一个数值。

l 数值计算方法AB * CD ＝ x3*x4 + y3*y4。

l 几何意义AB * CD = |AB| * |CD| * cos(a)，a为向量AB逆时针转向CD的角度，0<=a<360，也可以认为是两向量的夹角，0<=a<=180。一般用于求夹角，a = acos( (AB * CD) / (|AB| * |CD|) )。也可：|CD| * cos(a) = AB * CD / |AB|，即向量CD在AB上的投影。

四、叉积。

l 叉积的结果为一个向量。

l AB×CD数值= x3*y4 – x4*y3。

l 方向由右手螺旋定则判定。

l 几何意义：AB×CD = |AB| * |CD| * sin(a)，取绝对值即是以AB和CD为边的平行四边形面积。

五、线段相交的判定（判定线段AB和线段CD是否相交，属于哪种相交）。

l 规范相交：交点只有一个，且不是线段的端点。

1、充要条件：点A和点B在CD的两侧并且点C和点D在AB的两侧。

2、如何判断点A在向量CD的左侧还是右侧：CA×CD的数值大于0则左侧，小于0则右侧。于是CA×CD的数值和CB×CD的数值异号，点A和点B在CD的两侧，同理可判断点C和点D是否在AB的两侧（注意，必须严格异号）。

l 不规范相交：交点为某个线段的端点，甚至两线段有一段重合。

1、充要条件，一条线段的一个端点在另一条线段上。即A在CD上或B在CD上或C在AB上或D在AB上。

2、A在线段CD上的充要条件：

a) AC×AD = 0，几何意义，AC和AD组成的平行四边形的面积为0，即A、C、D三点共线。

b) A点在CD之间，A.x处于C.x和D.x之间并且A.y处于C.y和D.y之间

LSM 2007-08-21 08:52 发表评论