C++博客-oyjpArt ACM/ICPC算法程序设计空间-随笔分类-ACM/ICPChttp://www.cppblog.com/sicheng/category/2855.html// I am new in programming, welcome to my blog <br> I am oyjpart(alpc12, 四城) zh-cnSun, 06 Jul 2008 07:28:58 GMTSun, 06 Jul 2008 07:28:58 GMT60Topcoder 北京时间表http://www.cppblog.com/sicheng/archive/2008/07/06/55467.htmloyjpartoyjpartSun, 06 Jul 2008 07:12:00 GMThttp://www.cppblog.com/sicheng/archive/2008/07/06/55467.htmlhttp://www.cppblog.com/sicheng/comments/55467.htmlhttp://www.cppblog.com/sicheng/archive/2008/07/06/55467.html#Feedback0http://www.cppblog.com/sicheng/comments/commentRss/55467.htmlhttp://www.cppblog.com/sicheng/services/trackbacks/55467.htmlSRM410 2008-07-19 24:00
SRM411 2008-07-24 09:00
SRM412 2008-07-31 19:00
SRM413 2008-08-06 23:00
SRM414 2008-08-16 24:00
SRM415 2008-08-27 09:00





oyjpart 2008-07-06 15:12 发表评论
]]>SRM406 PTS500 FoldThePaperhttp://www.cppblog.com/sicheng/archive/2008/06/18/53831.htmloyjpartoyjpartWed, 18 Jun 2008 03:29:00 GMThttp://www.cppblog.com/sicheng/archive/2008/06/18/53831.htmlhttp://www.cppblog.com/sicheng/comments/53831.htmlhttp://www.cppblog.com/sicheng/archive/2008/06/18/53831.html#Feedback0http://www.cppblog.com/sicheng/comments/commentRss/53831.htmlhttp://www.cppblog.com/sicheng/services/trackbacks/53831.html

Problem Statement

     You have a rectangular piece of paper that's divided into 1x1 cells, each of which has an integer value. The paper will be described by a vector <string> paper. The ith element of paper will be a space delimited list of integers, where the jth integer of the ith element of paper represents the value of the jth cell of the ith row of the paper.



You want to perform a sequence of folds on the paper, where you may fold anywhere along an axis that is in between two rows or columns of the paper. After performing a fold, we wish to model the folded paper as a new, flat piece of paper. We will do this by considering two overlapping cells as a single cell, with a value that is the sum of the individual cells.



You wish to perform a sequence of folds such that the value of some single cell in the resulting piece of paper is as large as possible. Return this value.

Definition

    
Class: FoldThePaper
Method: getValue
Parameters: vector <string>
Returns: int
Method signature: int getValue(vector <string> paper)
(be sure your method is public)
    

Constraints

- paper will contain between 1 and 12 elements, inclusive. - Each element of paper will be a single-space delimited list of integers with no leading or trailing spaces. - Each element of paper will contain between 1 and 12 integers, inclusive. - Each element of paper will contain the same number of integers. - Each element of paper will contain between 1 and 50 characters, inclusive. - Each integer in paper will be between -100 and 100, inclusive. - Each integer in paper will have no leading zeros. - An integer in paper equal to zero will not have a preceding negative sign.

Examples

0)
     
{
"1 1 1",
"1 1 1"
}
 
Returns: 6
We can collapse every cell onto the upper-left cell.
1)
     
{
"1 -1",
"1 -1"
}
 
Returns: 2
We should perform only the fold between the two rows, and take the resulting left column.
2)
     
{
"1 -1 -1 1",
"-1 -1 -1 -1",
"-1 -1 -1 -1",
"1 -1 -1 1"
}
 
Returns: 4
Folding between the middle rows then the middle columns allows us to combine the four corner cells.
3)
     
{
"20 13 -2 100",
"-12 0 4 -3",
"4 1 -36 21"
}
 
Returns: 131

4)
     
{
"0"
}
 
Returns: 0

This problem statement is the exclusive and proprietary property of TopCoder, Inc. Any unauthorized use or reproduction of this information without the prior written consent of TopCoder, Inc. is strictly prohibited. (c)2003, TopCoder, Inc. All rights reserved.


题目大意是有一个12*12的矩阵，现在可以对这个矩阵横向或纵向折叠，出在重叠位置的数相加。
求折叠过程中任意位置产生的最大数。

很多大牛fail了，我一个DFS+剪枝也超时了，一共32人pass sys test，1000pts无人ac，此套题难度还是很大的。

基本思路是状态压缩DP，横向(1<<12)*纵向(1<<12)*加和。

但是这样会超时。关键是没有利用到折叠的信息。

预先生成某个位置的状态（由那些位置叠加而来），就可以减少检查量，就可以ac了。

如何生成这些状态呢？没错，又是一个DP. 呵呵。




oyjpart 2008-06-18 11:29 发表评论
]]>点错题，看错题，做错题！可以去死了！！！！！！！！！http://www.cppblog.com/sicheng/archive/2008/06/05/52300.htmloyjpartoyjpartThu, 05 Jun 2008 13:54:00 GMThttp://www.cppblog.com/sicheng/archive/2008/06/05/52300.htmlhttp://www.cppblog.com/sicheng/comments/52300.htmlhttp://www.cppblog.com/sicheng/archive/2008/06/05/52300.html#Feedback0http://www.cppblog.com/sicheng/comments/commentRss/52300.htmlhttp://www.cppblog.com/sicheng/services/trackbacks/52300.html2.看错题：交了发现居然还有220多分，原来是500分的题。汗。。。做250的吧。。估计交了也没多少分了。管他咯，做咯。结果怎么做怎么不对，样例死都调不过。到最后才发现题目都读错了。。一个数是上面的两个数之和，我看成上面的右边所有数的和了。。题目说的是upper right neiboring digits,就那个digit's'，汗，又读错了，亏我还说自己英语好。真tm丢脸。
3.唉 就这样就算了，结果呢，第二题又被cha掉了。汗。我的心那个叫苍凉啊。
4.到practise room 一测俺的第二题，全是 "Received: caught signal SIGKILL"，那个叫不解啊。我在自己机器上G++运行的好好的，干嘛的就SIGKILL了呢。。唉。。

可以去死了！！！！！！



oyjpart 2008-06-05 21:54 发表评论
]]>线性筛素数方法http://www.cppblog.com/sicheng/archive/2008/06/04/52172.htmloyjpartoyjpartWed, 04 Jun 2008 08:39:00 GMThttp://www.cppblog.com/sicheng/archive/2008/06/04/52172.htmlhttp://www.cppblog.com/sicheng/comments/52172.htmlhttp://www.cppblog.com/sicheng/archive/2008/06/04/52172.html#Feedback6http://www.cppblog.com/sicheng/comments/commentRss/52172.htmlhttp://www.cppblog.com/sicheng/services/trackbacks/52172.html
const int N = 25600000;
bool a[N];
int p[N];
int n;

void Prime1() {
    memset(a, 0, n * sizeof(a[0]));
    int num = 0, i, j;
    for(i = 2; i < n; ++i) if(!a[i]) {
        p[num++] = i;
        for(j = i+i; j < n; j +=i) {
            a[j] = 1;
        }
    }
}

void Prime2() {
    memset(a, 0, n*sizeof(a[0]));
    int num = 0, i, j;
    for(i = 2; i < n; ++i) {
        if(!(a[i])) p[num++] = i;
        for(j = 0; (j<num && i*p[j]<n); ++j) {
            a[i*p[j]] = 1;
            if(!(i%p[j])) break;
        }
    }
}

测试:

筛 [0, 100000) 范围内的素数
第一种素数筛法 0 毫秒
第二种素数筛法 0 毫秒

筛 [0, 200000) 范围内的素数
第一种素数筛法 15 毫秒
第二种素数筛法 0 毫秒

筛 [0, 400000) 范围内的素数
第一种素数筛法 16 毫秒
第二种素数筛法 15 毫秒

筛 [0, 800000) 范围内的素数
第一种素数筛法 47 毫秒
第二种素数筛法 16 毫秒

筛 [0, 1600000) 范围内的素数
第一种素数筛法 62 毫秒
第二种素数筛法 63 毫秒

筛 [0, 3200000) 范围内的素数
第一种素数筛法 297 毫秒
第二种素数筛法 109 毫秒

筛 [0, 6400000) 范围内的素数
第一种素数筛法 922 毫秒
第二种素数筛法 266 毫秒

筛 [0, 12800000) 范围内的素数
第一种素数筛法 2187 毫秒
第二种素数筛法 563 毫秒

筛 [0, 25600000) 范围内的素数
第一种素数筛法 4828 毫秒
第二种素数筛法 1187 毫秒

证明：任何一个合数只被标记一次。
      可以试着执行下这个程序的流程，就明白了

怎么样 还行吧？
什么，觉得这个程序效率上没多大提升，没有什么用？
把a[]改成int类型，然后
void Prime2() {
    memset(a, 0, n*sizeof(a[0]));
    int num = 0, i, j;
    for(i = 2; i < n; ++i) {
        if(!(a[i])) p[num++] = i;
        for(j = 0; (j<num && i*p[j]<n && (p[j]<=a[i]||a[i]==0)); ++j) {
            a[i*p[j]] = p[j];
        }
    }
}
这样一来a[i]将记录i的最小质因子
那么[0, n)内的数的因式分解就可以... 嘿嘿
o(质因子个数)求任意数因式分解：
void factor(int x) {
    while(a[x] != 0) {
        printf("%d\n", a[x]);
        x /= a[x];
    }
    printf("%d\n", x);
} 

然后用这个做了上次杭州比赛的GCD那题，虽然其实就是个容斥原理，可是我等白菜就是不会做。唉。
第一名8题，我们4题，这个差距大的有点想吐。
题目http://acm.hdu.edu.cn/showproblem.php?pid=1695
//Solution by alpc12:

#include <string.h>
#include <stdio.h>

const int N = 100010;

typedef __int64 LL;
#define I64Format "%I64d\n"
inline int count(int x) {int ret = 0;while(x != 0) {ret ++; x &= (x-1);}return ret;}

int a[N], p[18000];

void pre() {
    memset(a, 0, sizeof(a));
    int num = 0, i, j;
    for(i = 2; i < N; ++i) {
        if(!a[i])  p[num++] = i;
        for(j = 0; j < num && i * p[j] < N && (p[j]<=a[i] || a[i]==0); ++j) {
            a[p[j] * i] = p[j];
        }
    }
}

void go(int x, int y) {
    if(x == 0) { printf("0\n"); return; }
    int i, j;
    LL ans = 0;
    for(i = 1; i <= y; ++i) {
        if(!a[i]) {
            ans += (i<=x?(i-1):x);
        } else {
            int fac[20], nfac = 0, z = i;
            while(a[z] != 0) {
                fac[nfac++] = a[z];
                z /= a[z];
            }
            fac[nfac++] = z;
            int k = 1;
            for(j = 1; j < nfac; ++j) {
                if(fac[j] != fac[j-1])
                    fac[k++] = fac[j];
            }
            nfac = k;
            int now = 0;
            int xx = x;
            if(x >= i) xx=i-1;
            int mask;
            for(mask = 1; mask < (1<<nfac); ++mask) {
                int d = count(mask), mul = 1;
                for(j = 0; j < nfac; ++j) if((mask&(1<<j)) != 0) {
                    mul *= fac[j];
                }
                if(d&1) now += xx/mul;
                else now -= xx/mul;
            }
            ans += xx-now;
        }
    }
    printf(I64Format, 1+ans);
}

int main() {

//    freopen("t.in", "r", stdin);

    int ntc, a, b, c, d, k;
    scanf("%d", &ntc);
    pre();
    int tc = 0;
    while(ntc--) {
        printf("Case %d: ", ++tc);
        scanf(" %d %d %d %d %d", &a, &b, &c, &d, &k);
        if(k == 0) printf("0\n");
        else {
            a = b/k;
            b = d/k;
            if(a > b) go(b, a);
            else go(a, b);
        }
    }
    return 0;
}




oyjpart 2008-06-04 16:39 发表评论
]]>基本参数搜索http://www.cppblog.com/sicheng/archive/2008/06/03/52034.htmloyjpartoyjpartTue, 03 Jun 2008 07:45:00 GMThttp://www.cppblog.com/sicheng/archive/2008/06/03/52034.htmlhttp://www.cppblog.com/sicheng/comments/52034.htmlhttp://www.cppblog.com/sicheng/archive/2008/06/03/52034.html#Feedback13http://www.cppblog.com/sicheng/comments/commentRss/52034.htmlhttp://www.cppblog.com/sicheng/services/trackbacks/52034.html

上次百度之星第三题竟然不会做，很是惭愧啊，脑袋生锈了。

后来从HUST上面找了道类似的题目，AC了。


The perfect hamilton path

Time Limit: 5 Sec  Memory Limit: 128 MB
Submissions: 72  Solved: 16

Description

There are N(2 <= N <= 13) cities and M bidirectional roads among the cities. There exist at most one road between any pair of the cities. Along every road, there are G pretty girls and B pretty boys(1 <= G,B <= 1000).
You want to visit every city exactly once, and you can start from any city you want to. The degree of satisfaction is the ratio of the number of the pretty girls to the number of the pretty boys. You want to know the highest degree of satisfation.

Input

There are multiply test cases.
First line: two integers N, M;
The following M lines: every line with four integers i, j, G, B, response that there is a road between i and j with G and B.

Output

The highest degree of the satisfation, rounded to the third place after the decimal point.

Sample Input

3 3
1 2 5 3
2 3 7 4
3 1 13 11

Sample Output

1.714

HINT

Source

dupeng


题目的意思是找到一个sigma(G)/sigma(B)最大的hamilton回路。
典型的参数搜索。二分或者迭代答案就可以了。

Solution：

#include <stdio.h>
#include <queue>
#include <cmath>
using namespace std;

const double EPS = 1e-4;
const int N = 15;
const int M = N * N;

#define Max(a, b) (a>b?a:b)

inline int dblcmp(double a, double b) {
    if(fabs(a-b) < EPS) return 0;
    return a < b ? -1 : 1;
}

struct Node 
{
    int x, mask;
    double s;
    Node() {}
    Node(int mm, int xx, double ss) {
        x = xx;
        mask = mm;
        s = ss;
    }
};

int n, m;

double adj[N][N];
int X[M], Y[M], G[M], B[M];

double dp[1<<N][N];

double go(double ans) {
    int i, j;
    for(i = 0; i < n; ++i) {
        adj[i][i] = 0;
        for(j = i+1; j < n; ++j) {
            adj[i][j] = adj[j][i] = -10e300;
        }
    }
    for(i = 0; i < m; ++i) {
        adj[X[i]-1][Y[i]-1] = G[i]-ans * B[i];
        adj[Y[i]-1][X[i]-1] = adj[X[i]-1][Y[i]-1];
    }

    for(i = 0; i < (1<<n); ++i) {
        for(j = 0; j < n; ++j)
            dp[i][j] = -10e100;
    }
    queue<Node> Q;
    for(i = 0; i < n; ++i) {
        Q.push(Node(1<<i, i, 0.0));
        dp[1<<i][i] = 0;
    }
    while(Q.size()) {
        int f = Q.front().mask, x = Q.front().x;
        double s = Q.front().s;
        double& d = dp[f][x];
        Q.pop();
        if(s < d) continue;
        for(i = 0; i < n; ++i) if((f&(1<<i)) == 0) {
            if(dp[f|1<<i][i] < s + adj[x][i]) {
                dp[f|1<<i][i] = s + adj[x][i];
                Q.push(Node(f|1<<i, i, s + adj[x][i]));
            }
        }
    }

    double max = -10e100;
    for(i = 0; i < n; ++i) {
        max = Max(max, dp[(1<<n)-1][i]);
    }
    return max;
}

int main()
{
    // freopen("t.in", "r", stdin);

    int i;
    double ans;
    while(scanf("%d %d", &n, &m) != EOF) {
        double min = 2000, max = 0;
        for(i = 0; i < m; ++i) {
            scanf("%d %d %d %d", &X[i], &Y[i], &G[i], &B[i]);
            if(B[i] < min) min = B[i];
            if(G[i] > max) max = G[i];
        }
        double lo = 0, hi = max/min;
        int ok = 0;
        for(i = 0; ; ++i) {
            double mid = lo + (hi-lo)/2;
            if(dblcmp((ans=go(mid)), 0.0) > 0) {
                lo = mid;
            } else if(dblcmp(ans, 0.0) == 0) {
                printf("%.3lf\n", mid);
                ok = 1;
                break;
            } else {
                hi = mid;
            }
        }

        if(!ok) { int a = 0; a = 1/a; }
    }

    return 0;
}

 




oyjpart 2008-06-03 15:45 发表评论
]]>湖大邀请赛http://www.cppblog.com/sicheng/archive/2008/05/27/51254.htmloyjpartoyjpartTue, 27 May 2008 03:32:00 GMThttp://www.cppblog.com/sicheng/archive/2008/05/27/51254.htmlhttp://www.cppblog.com/sicheng/comments/51254.htmlhttp://www.cppblog.com/sicheng/archive/2008/05/27/51254.html#Feedback3http://www.cppblog.com/sicheng/comments/commentRss/51254.htmlhttp://www.cppblog.com/sicheng/services/trackbacks/51254.html阅读全文

oyjpart 2008-05-27 11:32 发表评论
]]>1周又过去了http://www.cppblog.com/sicheng/archive/2008/05/19/50349.htmloyjpartoyjpartMon, 19 May 2008 01:03:00 GMThttp://www.cppblog.com/sicheng/archive/2008/05/19/50349.htmlhttp://www.cppblog.com/sicheng/comments/50349.htmlhttp://www.cppblog.com/sicheng/archive/2008/05/19/50349.html#Feedback0http://www.cppblog.com/sicheng/comments/commentRss/50349.htmlhttp://www.cppblog.com/sicheng/services/trackbacks/50349.html爸爸，妈妈，女友三人在同一周过生日，真够我忙乎
软件工程大作业，由于去杭州比赛拖欠的一大笔作业都要赶
拼题任务做不完，结果导致我回家看望妈妈的期间其实都在电脑前赶着做题
还有一个上级命令下来的任务，根本没给我拒绝的机会，唉，这念头，任人宰割啊...
本来就很郁闷了，结果回长沙的火车还晚点，造成很大麻烦。。

昨天做了一场组队练习赛，结果令我非常不爽。
昨天的Substitution Cipher题，就是一个很简单的TOPO序，WA了。后来准备让alpc33敲Lattice Animals ，又怕他敲不完自己抢过键盘来敲。当时本来就时间不多，我不去细查Substitution Cipher 这题，真是决策错误。结束之后跟alpc55交流，结果是题目里面的输入有空格，要gets读，真TMD变态！昨晚改了下居然还WA，今早发现没有改读入目标串为gets，一改就过了。这道题真的做的我很郁闷。题目为什么不好好说呢！非要让你郁闷一下！当然如果那时候我细看题目，也许能分析出来A掉，可是我觉得题目真的没必要这样搞。怎么不学学topcoder，题目的样例都出的很厚道，把输入输出搞得明明白白。

还有就是我觉得alpc33要快点熟悉STL了，好像常见的容器都不怎么会用。在场上因为这个东西浪费很多时间+2次罚时。这个任务就交给我吧，哈哈。




oyjpart 2008-05-19 09:03 发表评论
]]>杭州邀请赛 ACM在何方?http://www.cppblog.com/sicheng/archive/2008/05/14/49858.htmloyjpartoyjpartWed, 14 May 2008 12:40:00 GMThttp://www.cppblog.com/sicheng/archive/2008/05/14/49858.htmlhttp://www.cppblog.com/sicheng/comments/49858.htmlhttp://www.cppblog.com/sicheng/archive/2008/05/14/49858.html#Feedback0http://www.cppblog.com/sicheng/comments/commentRss/49858.htmlhttp://www.cppblog.com/sicheng/services/trackbacks/49858.html杭州是个很美的城市，可是我却无心欣赏美景。
ACM啊 ACM啊 何止一个竞赛

废话少说。直接进入比赛正题。

Rank Team Solved Penalty 1001 1002 1003 1004 1005 1006 1007 1008 1009 1010
1 TEAM096
HyperHexagon
复旦大学		 8 23:12:04 02:15:56
(-1)		 00:28:33
(-1)		 04:39:03
(-2)		 01:51:31
(-1)		 02:02:47
(-1)
04:29:43
(-1)
03:25:40 00:58:51
(-2)
2 TEAM022
GCC
武汉大学		 7 18:29:12 02:42:16
(-1)		 00:42:41
(-1)		 (-2) 04:38:07 02:28:29
01:51:49
(-2)
03:28:51 01:16:59
3 TEAM106
Gungnir
浙江大学		 6 13:44:54 01:41:38
(-1)		 00:23:12
(-1)
03:24:43
(-1)

04:20:46 (-2) 01:53:28
(-1)		 00:41:07
4 TEAM051
Orbicular
浙江大学		 6 17:46:21 04:24:09 00:18:18
02:00:16
(-1)

02:43:26
(-1)		 (-3) 04:57:51
(-2)		 01:42:21
(-1)
5 TEAM010
Heroes
南京大学		 6 19:32:42 04:30:31
(-4)		 01:00:20
04:52:12

03:27:45
(-3)
01:39:06
(-1)		 00:42:48
(-2)
6 TEAM003
Zenith
复旦大学		 6 20:15:26 04:05:44
(-6)		 01:37:50
(-6)
02:28:09
(-1)		 04:18:06
(-1)


02:46:22 00:19:15
7 TEAM002
acm.hdu.edu.cn
杭州电子科技大学		 5 13:49:03 02:26:27
(-1)		 00:34:35
(-1)
04:34:11
(-3)

03:30:27 (-13)
01:03:23
8 TEAM101
RedField-Z+S
友情参赛队		 5 14:51:31 03:22:30
(-2)		 00:33:43

04:32:29


03:17:12
(-5)		 00:45:37
9 TEAM082
alpcT2
国防科技大学		 4 05:32:58 01:06:22 00:27:06 (-3) 02:55:47




01:03:43
10 TEAM043
Jackrabbit
复旦大学		 4 12:28:46 03:35:58 01:22:04
(-4)
04:27:29
(-2)




01:03:15
11 TEAM062
alpcT1
国防科技大学		 4 14:00:20 02:58:40
(-4)		 02:07:35
(-3)
04:09:14




02:04:51
(-1)
12 TEAM056
TEAM III
杭州电子科技大学		 4 14:40:03 (-2) 04:16:47
(-5)
03:44:16

(-4)
02:01:01
(-1)		 02:37:59
13 TEAM067
Backbone--
华东理工大学		 4 15:48:18 03:28:01
(-2)		 01:32:17
(-4)



(-2)
04:56:07
(-1)		 02:31:53
(-3)
14 TEAM100
AngryHair
浙江工业大学		 4 16:23:57 02:31:31
(-3)		 01:00:54
(-3)
04:53:01
(-6)



(-12) 02:58:31
(-3)
15 TEAM026
Bear
北京理工大学		 3 07:57:27 (-5) 01:02:45
(-2)





03:15:01
(-3)		 01:59:41
16 TEAM099
FZU1
福州大学		 3 08:03:52 03:51:55
(-4)		 00:39:42
(-1)
(-11)



(-9) 01:52:15
17 TEAM001
ECNU-01
华东师范大学		 3 08:09:53
04:54:56
(-2)
00:39:59




01:54:58
18 TEAM032
Swordsmen
杭州电子科技大学		 3 09:58:11
02:05:48
(-1)
03:44:58
(-2)




03:07:25
19 TEAM007
Eternal Arbiter
武汉大学		 3 11:38:15 01:33:16 03:25:01
(-11)
(-7)




02:59:58
20 TEAM015
TJU-Team2
天津大学		 3 11:43:31 04:16:24
(-2)		 00:59:35
(-11)



(-1) 03:47:32
(-6)
21 TEAM071
ppg
浙江师范大学		 3 14:28:04 04:31:19
(-2)		 04:26:09
(-7)
(-1)




02:10:36
(-1)
22 TEAM033
SkyFree
杭州电子科技大学		 3 14:58:18
03:13:16
(-6)
02:59:50
(-4)

(-2)

04:05:12
(-4)
23 TEAM110
warcraft
浙江师范大学		 3 15:02:35 04:16:09
(-5)		 03:50:08
(-11)
(-3)




01:36:18
24 TEAM077
hummer
华中科技大学		 2 01:05:33 (-4) 00:21:43

(-1)



00:43:50
25 TEAM116
TJU-Team1
天津大学		 2 02:37:46 (-18) 00:18:01 (-10)





01:39:45
(-2)
26 TEAM030
Splinglet
浙江大学宁波理工学院		 2 02:47:37 (-5) 00:39:21
(-1)
(-4)




01:28:16
(-1)
27 TEAM012
alpcT3
国防科技大学		 2 03:33:09 (-7) 00:54:38





02:38:31 (-5)
28 TEAM005
Ice-Cream
杭州电子科技大学		 2 03:46:20 (-5) 01:36:14
(-1)
(-1)




01:50:06
29 TEAM048
Coderseekers
华东理工大学		 2 04:10:39 (-4) 02:09:48
(-4)
(-2)




00:40:51
30 TEAM072
FFR
东南大学		 2 05:39:37
01:14:37
(-1)





04:05:00 (-6)
31 TEAM065
SunShine1
上海师范大学		 2 05:41:12 (-5) 01:54:27
(-3)
(-1)




02:46:45
32 TEAM060
TJU-Team3
天津大学		 2 05:48:44
02:25:19
(-4)		 (-4)





02:03:25
33 TEAM058
LifeCodere
浙江大学宁波理工学院		 2 05:54:40 (-10) 02:02:29
(-6)





(-3) 01:32:11
(-1)
34 TEAM094
Bug2
上海大学		 2 06:19:50 (-3) 02:45:04
(-3)		 (-3)





02:14:46
(-1)
35 TEAM075
xYtForWard
浙江大学		 2 06:21:59 (-7) 00:24:45






04:37:14
(-4)
36 TEAM038
SkyLine
湖州师范学院		 2 07:13:29


03:15:17
(-4)




02:18:12
(-1)
37 TEAM091
ECNU-03
华东师范大学		 2 07:15:58
01:37:09
(-2)
04:38:49
(-1)




(-6)
38 TEAM037
FT
东南大学		 2 07:18:49
00:54:19
(-1)






04:24:30
(-5)
39 TEAM078
KOP
华东师范大学软件学院		 2 07:25:00
03:17:42
(-3)






02:47:18
(-1)
40 TEAM089
X_BOX
华东交通大学		 2 07:32:26
02:17:38
(-4)






03:34:48
(-1)
41 TEAM081
Killers
南京大学		 2 07:45:46
01:30:10
(-2)





(-1) 04:35:36
(-3)
42 TEAM009
Dota
华东理工大学		 2 07:47:02 (-1) 02:53:34
(-6)






01:53:28
(-3)
43 TEAM017
Thinkmore
华东师范大学软件学院		 2 07:49:34 (-5) 01:50:43
(-4)






03:58:51
(-2)
44 TEAM014
GreatGrass
华中科技大学		 2 07:52:21 (-2) 03:18:30
(-4)






02:13:51
(-3)
45 TEAM042
dragon
浙江师范大学		 2 08:21:09
03:36:19
(-7)
(-1)




02:24:50
46 TEAM059
Cheers
浙江工业大学		 2 08:23:44 (-4) 01:19:23
(-4)






04:24:21
(-4)
47 TEAM105
zust-Rookie
浙江科技学院		 2 09:49:31
03:21:13
(-5)





(-1) 03:28:18
(-4)
48 TEAM064
ZstuYs
浙江理工大学		 2 10:27:59
02:48:47
(-5)






04:39:12
(-4)
49 TEAM057
Slash
武汉大学		 2 11:06:07 (-2) 02:46:46
(-8)






03:59:21
(-5)
50 TEAM084
WZU
温州大学		 2 11:29:04
04:49:00
(-15)






01:40:04
51 TEAM024
Miracle2008
浙江工业大学		 2 12:49:13
04:42:11
(-13)






03:07:02
(-2)
52 TEAM085
Daydayup
杭州电子科技大学		 2 13:37:42
04:27:46
(-10)






04:09:56
(-5)


会议记录：昨晚我和吴家旭商量，让锻炼新人，所以我们都让新人做了第一题。我 们选择了B题，一个错误的算法上去WA了，我的责任，因为是我和llz合计的算法， WA之后看了他的代码找到了一处错误，交上去还是WA。后来我就上去敲J题，因为J 题过了1个人，B题就交给他俩，后来发现是算法错误。我的J题敲完之后WA，调了一 下，随后我就下来看代码了，后来我发现自己有错误，改正之后AC了，此时时间已经 过去两个小时，一直在郁闷的我们看到了一点希望。很快B题在改正之后过了。然后 我们寻找下一道可以做的题，我们看中了A题（过的人最多），经过讨论llz认同了我 的做法，我就上去敲了，这题敲的非常不好，我要做深刻的检讨，虽然当时我不慌张 ，可是由于几个莫名其妙的RE，让我浪费了大量罚时。后来发现只是一个数组下标越 界。但是过了之后，我们的形势还是有所好转，然后我开始看别D题、觉得可以做没 想的很清楚，我就跟林杨飞讨论一下，我觉得时间不多我就上去敲了，封榜后1A了， 我们很惊喜。这个时候的时间是2：22，我对我们组做了一个决策，去做场上过的最多 的I和G题，虽然llz对E有想法，但是他没想清楚我没让他上去敲，接着我推了I题发现属 于我能力之外，于是我看看G题，我和llz觉得是一个很弱的欧拉函数，llz就上去写 欧拉函数产生器、写着写着发现这样做不对，然后我又回到了I，llz继续想G，这样到 比赛结束。
个人而言，我觉得我属于慢热型选手，比赛开始之后很久才进入状态，平常也经常暴 露这种问题；或者说是自己状态调整的不好，如果在组队的时候考虑让比较稳的人敲 第一题，同时我要加强自己写题时的大脑清醒度。
对于队伍，我们今天的决策大部分都是我做的，我们敲题的顺序完全没有错误，应该说我们敲出四题是正常水平，可是在罚时上有很大的空间，如果能早些出题，就可以考虑E题了，就有可能出5题。

感言：
第一名做出8题。我佩服。但同时也给了我巨大的打击。
原来差距如此之大。我也不知道自己怎么了，那么简单的A题，runtime error几次，还以为是queue没有清空，结果就是个白痴的数组越界。
今天师傅DaringQQ也来比赛了，果然不敌师傅。

给自己的训练目标
1.脱离SCL(Source Code Library)做题。强大代码能力。
2.多比赛 训练比赛时的大脑清醒度。
3.训练稳定性。向alpc02大牛学习，一个让人放心的选手，队伍的核心。

For ALPCs：
1.首先在这里表扬下alpcT2::alpc52，不慌不忙却又思维活跃地1Y掉J题和D题，非常牛B。
2.alpc50,alpc42还需要比赛的磨练，在赛场上好好发挥。
3.所有ALPCs都没有推出I题公式，这怎么能行呢？



oyjpart 2008-05-14 20:40 发表评论
]]>SRM401 550 PTS ParticleCollisionhttp://www.cppblog.com/sicheng/archive/2008/05/07/49083.htmloyjpartoyjpartTue, 06 May 2008 18:32:00 GMThttp://www.cppblog.com/sicheng/archive/2008/05/07/49083.htmlhttp://www.cppblog.com/sicheng/comments/49083.htmlhttp://www.cppblog.com/sicheng/archive/2008/05/07/49083.html#Feedback0http://www.cppblog.com/sicheng/comments/commentRss/49083.htmlhttp://www.cppblog.com/sicheng/services/trackbacks/49083.html

Problem Statement

    

Particles (which can be considered points in 3D-space for the purposes of the problem) can move in an electro-magnetic field. If a particle is charged, its trajectory can be described as spiral, and if it is uncharged, it is just a straight line. Given two particles (one charged and one uncharged) it should be determined whether they can possibly collide or not. Two particles can possibly collide if and only if their trajectories intersect.

Some steps have already been made by the physicist to simplify the problem, so the coordinates of the charged particle are represented as follows:

x1 = cos(PI * t)

y1 = sin(PI * t)

z1 = t

and for the uncharged particle:

x2 = vx * t + x0

y2 = vy * t + y0

z2 = vz * t + z0

Here t is a parameter which can be chosen arbitrarily and independently for both trajectories.

Your method will be given 6 integers - vx, vy, vz, x0, y0 and z0, describing the trajectory of the uncharged particle. It should determine whether the two given trajectories intersect or not. If they do, it should return a vector <double> containing exactly 3 elements x, y and z - the coordinates of the point where a collision can happen. If there is more than one such point, it should return a vector <double> containing exactly three zeroes. If collision of the two particles is impossible it should return an empty vector <double>.

Definition

    
Class: ParticleCollision
Method: collision
Parameters: int, int, int, int, int, int
Returns: vector <double>
Method signature: vector <double> collision(int vx, int vy, int vz, int x0, int y0, int z0)
(be sure your method is public)
    

Notes

- PI can be considered equal to 3.14159265358979323846. - All return values with either an absolute or relative error of less than 1.0E-9 are considered correct.

Constraints

- vx, vy and vz will each be between -10 and 10, inclusive. - x0, y0 and z0 will each be between -10 and 10, inclusive.

Examples

0)
     
0
 
0
 
0