七月流水账

7.1-7.3

7.4

看了一会儿李老大给的Caffe的教程

重配python + qt + vtk，大概是原来装了一个xxx（忘记名字），然后位数有些不对，搞了好一阵

果然还是卸载重装好><

然后又看了会儿Matlab神经网络工具箱怎么用

还没看太明白就被抓去搞界面了>w<（终于稍微有点点懂怎么把代码和GUI的东西结合起来了，实现还不太明白）

间隙看了线性筛

7.5

读取dcm的标签信息，需要先配置好simplevtk

import SimpleITK as sitkinputImage = sitk.ReadImage(str(file_name))image_laterality = inputImage.GetMetaData('0020|0060')print image_laterality

View Code

需要什么信息先用ImageJ打开图片看下对应标签的编号

然后感觉是VTK这一套东西都不懂，做起事情来很茫然..

---------------------------------------------------------------------

晚上和以前一起实习的聊了聊，发现好多好多都没有学阿

计算机网络，操作系统，计算机组成原理，编译原理算基础类的过一遍

语言类的C++ primer之类的过一遍

linux XXXXX编程

机器学习的常用算法，用过还有会推导，用到什么框架

leetcode

（web服务器xxxxxxxx）

7.6

准备老王中 & 继续完善软件的功能

创建一个新的界面

7.7

觉得周围的人都很厉害，自己菜的不行...

DCNN原来是由多个CNN组成的模型

区域生长

7.8

LDA 将高维样本数据采用线性组合的方式投影到一维空间，并找到使得类内间距最小，类间间距最大的组合系数

划了一天水叭

7.9

难过了好久...不知道该怎么办于是去实验室接着老王了（其实是蹭网蹭空调蹭水喝><

主体部分大概好了...但是总觉得还差很多TAT

7.10

放射组学

------- 昏割线 -------------------------------------

剪掉头发重新好好当狗了要加油喔><

-------------------------------------------------------

使用PyQt 编写第一个Python GUI程序

但是教程给出的示例代码在我的电脑会报错

原因是电脑上装的Qt建一个ui的时候，会搞出一个发送者和接受者的信号，和槽，但是在代码里面没有写到

所以再这样加上两个函数就好

 1 import sys
 2 from PyQt4 import QtCore, QtGui, uic
 3 
 4 qtCreatorFile = "C:\\Users\\lab\\Desktop\\test\\haha.ui"  # Enter file here.
 5 
 6 Ui_MainWindow, QtBaseClass = uic.loadUiType(qtCreatorFile)
 7 
 8 class MyApp(QtGui.QMainWindow, Ui_MainWindow):
 9     def __init__(self):
10         QtGui.QMainWindow.__init__(self)
11         Ui_MainWindow.__init__(self)
12         self.setupUi(self)
13 
14     def accept(self):
15         print "111"
16 
17     def reject(self):
18         print "22"
19 
20 
21 if __name__ == "__main__":
22     app = QtGui.QApplication(sys.argv)
23     window = MyApp()
24     window.show()
25     sys.exit(app.exec_())

View Code

Python调用另外的py文件中的类

7.11

编写 UI 和后端的逻辑！

python中的编码问题有点长先码在这里叭

具体现在遇到的问题也很好解决，就是输出中文字符串乱码的话，在前面加一个u转为unicode编码就能正常显示了

真的好久没有摔过狗吃屎了哈哈哈哈哈，只是感觉手指现在打代码有点抬不起qwq

RNN 今晚讲课有讲，但是感觉听的很恍惚...先码叭

Problem A. Ambiguous Cipher

按照奇偶来想( 按照样例来想...奇数无解，偶数可以知道顺二和倒二位置的字母，然后从左到右推一遍再从右到左推一遍

wa的不明所以的原因是，字符串每次输出答案是直接printf 完整个字符串，但是可能有上一个case的在后面(所以就打印到长度就好>o<好傻喔

 1 #include <cstdio>
 2 #include <cstring>
 3 #include <iostream>
 4 #include <algorithm>
 5 using namespace std;
 6 typedef long long LL;
 7 const int maxn = 1e5+5;
 8 
 9 char s[maxn];
10 
11 int main()
12 {
13     int T;
14     freopen("A-large-practice (1).in","r",stdin);
15     freopen("A.txt","w",stdout);
16     scanf("%d",&T);
17     for (int kase = 1;kase <= T;kase++)
18     {
19         //scanf("%s",s+1);
20         cin >> (s+1);
21         int len = strlen(s+1);
22         printf("Case #%d: ",kase);
23         if (len%2) {puts("AMBIGUOUS");}
24         else
25         {
26             char t[105];
27             t[2] = s[1];
28             t[len-1] = s[len];
29 
30             for (int i = 4;i <= len;i += 2)
31             {
32                 t[i] = (s[i-1] - 'A'  - (t[i-2] - 'A' ) + 26) % 26 + 'A';
33             }
34             for (int i = len-3;i >= 1;i -= 2)
35             {
36                 t[i] = ( (s[i+1] - 'A') - (t[i+2] - 'A') + 26) % 26 + 'A';
37             }
38             
39             for (int i = 1;i <= len;i++) printf("%c",t[i]);
40             puts(" ");
41         }
42     }
43     return 0;
44 }

View Code

7.12 >_<

论文中...发现有好几天没看操作系统了....要割的感觉

7.13 >_<

论文中 & 看了一些相关论文关于实验结果的评估标准（感觉自己看东西都很不仔细...

7.14 >_<

Python roc

还是不懂qt的一些控件该怎么用

7.15 >_<

大早上兴高采烈出门结果悲剧咯><

中午看了51nod一个教学视频！

然后还是不懂Qt的那些button要怎么搞成一个数组一样- -| (

Project Euler

Multiples of 3 and 5 第一次交PE的题目！纪念！暴力或者容斥

7.16 >_<

原来py的大树算 2^2009这样只要 2 ** 2009 就好了呀！

py字符串拼接成命令执行

已经第很多次这个错误了：IndentationError: unindent does not match any outer indentation level

可以把代码粘到notepad++看看（显示空格和tab键）

button数组可以加到一个group里面就好，然后界面类写的差不多之后应该是一些细节的问题叭><

Even Fibonacci numbers 小于4000000的偶数斐波那契数的和

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
typedef long long LL;
const int maxn = 1e5+5;LL f[100];int main()
{f[1] = 1;f[2] = 2;for (int i = 3;i < 100;i++){f[i] = f[i-1] + f[i-2];if (f[i] > 4000000) break;}LL sum = 0LL;for (int i = 1;i < 100;i++){if (f[i] % 2 == 0) sum += f[i];}printf("%lld\n",sum);return 0;
}

View Code

Largest prime factor 600851475143的最大素因子

 1 #include<cstdio>
 2 #include<cstring>
 3 #include<cmath>
 4 #include<iostream>
 5 #include<algorithm>
 6 #include<vector>
 7 #include<set>
 8 using namespace std;
 9 
10 typedef long long LL;
11 const int maxn = 100005;
12 int dp[maxn];
13 int Max[maxn];
14 vector<int> p;
15 
16 int vis[maxn];
17 int T,a[105];
18 int cnt[maxn];
19 set<long long > s;
20 LL n;
21 
22 void solve( long long  x){
23     LL y = x;
24     for(long long  j = 2;1LL*j*j <= y;j++){
25         if(y%j == 0){
26             while(y%j == 0){
27                 s.insert(j);
28                 y = y/j;
29             }
30         }        
31     }
32     if(y > 1) s.insert(y);
33     
34     for(set<long long >::iterator it = s.begin();it != s.end();++it){
35         printf("*it = %d\n",*it);
36     }
37     
38 }
39 
40 
41 int main(){
42     while(scanf("%I64d",&n) != EOF){
43         s.clear();
44         solve(n);
45     }
46     return 0;
47 }

View Code

周杰伦日！撒花！

7.17 >_<

好像论文字数凑太多（要开始删一删和提炼一下，白翻那么久qwq

字符串最小表示法之前有看，但今天才板子敲了一敲

题目

 1 #include <cstdio>
 2 #include <cstring>
 3 #include <iostream>
 4 #include <algorithm>
 5 using namespace std;
 6 typedef long long LL;
 7 const int maxn = 5e6+5;
 8 
 9 char s[maxn];
10 
11 int getmin(char s[])
12 {
13     int i = 0,j = 1,k = 0;
14     int len = strlen(s);
15     while(i < len && j < len && k < len)
16     {
17         int t = s[(i+k) % len] - s[(j+k) % len];
18         if(t == 0) k++;
19         else
20         {
21             if(t > 0) i += k+1;
22             else j += k+1;
23             if (i == j) j++;
24             k = 0;
25         }
26     }
27     return i < j ? i : j;
28 }
29 
30 int main()
31 {
32     int T;
33     scanf("%d",&T);
34     while(T--)
35     {
36         scanf("%s",s);
37         printf("%d\n",getmin(s));
38     }
39     return 0;
40 }

View Code

这周的讨论课讲了深度学习...可能是我一点没学习了解过，所以感觉大部分是在讲一个框架...而没有通过一个小的demo讲工作原理- -

7.18 >_<

早上看到一个月色真美的结局的回答，戳中...瞬间泪目qwq

数据集不平衡

继续界面，看完影像报告之后发现还要加好多东西- -

隐藏一个控件（百度谷歌中英文......终于搜到了...哭）

Weakness and Poorness 三分加前缀和，讲可以三分那一段想起白薯上的三分例题><

还没写..

Problem B. X Squared

看好几天..终于明白做法了，主要是那个矩形的定义那句话没读明白所以搞的一头雾水..

矩形是指，只有四角有'X'的矩形

前提条件是只有一行一列是只有一个X的，然后其他行其他列都得有两个X

可以从最终状态来看，是有(n-1)个矩形的，无论怎么换位置，都不会破坏掉矩形（所以这个算充分性？）

必要性如果一个矩形是破损的话，但是无论怎么换X都是不能破坏矩形的所以矛盾了？........不会证明....

 1 #include <cstdio>
 2 #include <cstring>
 3 #include <iostream>
 4 #include <algorithm>
 5 #include <map>
 6 #include <vector>
 7 using namespace std;
 8 typedef long long LL;
 9 const int maxn = 1e5+5;
10 
11 char g[105][105];
12 int n,row[105],col[105],vis[105][105];
13 
14 int main()
15 {
16     int T;
17     freopen("B-large-practice.in","r",stdin);
18     freopen("B2.txt","w",stdout);
19     scanf("%d",&T);
20     for (int kase = 1;kase <= T;kase++)
21     {
22         scanf("%d",&n);
23         memset(row,0,sizeof(row));
24         memset(col,0,sizeof(col));
25         memset(vis,0,sizeof(vis));
26         for (int i = 1;i <= n;i++) 
27         {
28             scanf("%s",g[i]+1);
29             for (int j = 1;j <= n;j++)
30             {
31                 if (g[i][j] == 'X') row[i]++,col[j]++;
32             }
33         }
34         printf("Case #%d: ",kase);
35         int c1 = 0,c2 = 0,d1 = 0,d2 = 0;
36         for (int i = 1;i <= n;i++)
37         {
38             if (row[i] == 1) c1++;
39             if (row[i] == 2) c2++;
40             if (col[i] == 1) d1++;
41             if (col[i] == 2) d2++;
42         }
43         //printf("c1 = %d c2 = %d d1 = %d d2 = %d\n",c1,c2,d1,d2);
44         if (c1 != 1 || c2 != n-1 || d1 != 1 || d2 != n-1) {puts("IMPOSSIBLE");continue;}
45 
46         vector<pair<int,int> > h[105];
47         for (int j = 1;j <= n;j++)
48         {
49             int p1 = 0,p2 = 0;
50             for (int i = 1;i <= n;i++)
51             {
52                 if (g[i][j] == 'X')
53                 {
54                     if (p1 == 0) p1 = i;
55                     else p2 = i;
56                 }
57             }
58             h[j].push_back(make_pair(p1,p2));
59         }
60 
61         int cnt = 0;
62         for (int i = 1;i <= n;i++)
63         {
64             if (row[i] == 1) continue;
65             int p1 = 0,p2 = 0;
66             for (int j = 1;j <= n;j++)
67             {
68                 if (g[i][j] == 'X')
69                 {
70                     if (p1 == 0) p1 = j;
71                     else p2 = j;
72                 }
73             }
74             //printf("i = %d p1 = %d p2 = %d\n",i,p1,p2);
75             pair<int,int> l = h[p1][0];
76             pair<int,int> r = h[p2][0];
77             if (l.first == r .first && l.second == r.second && !vis[p1][p2])
78             {
79                 cnt++;
80                 vis[p1][p2] = 1;
81             }
82         }
83         //printf("cnt = %d\n",cnt);
84         if (cnt == (n-1)/2) puts("POSSIBLE");
85         else puts("IMPOSSIBLE");
86 
87     }
88     return 0;
89 }

View Code

7.19 >_<

论文删删改改...

然后界面现在遇到的问题是，把一些Button在初始化的时候连接起来，没有被连接的Button们就只能选中一个了，尝试写另一个函数可是没有卵用...想哭

D. "Or" Game

还有点没想明白...

Problem C. Magical Thinking

小数据暴力2^10或者贪心

 1 #include <cstdio>
 2 #include <cstring>
 3 #include <iostream>
 4 #include <algorithm>
 5 using namespace std;
 6 typedef long long LL;
 7 const int maxn = 1e5+5;
 8 
 9 char s[105],t[105];
10 int n,q,x;
11 
12 int main()
13 {
14     int T;
15     freopen("C-small-practice.in","r",stdin);
16     freopen("Csmall.txt","w",stdout);
17     scanf("%d",&T);
18     for (int kase = 1;kase <= T;kase++)
19     {
20         scanf("%d %d",&n,&q);
21         scanf("%s",s+1);
22         scanf("%s",t+1);
23         scanf("%d",&x);
24         printf("Case #%d: ",kase);
25         int a = 0,d = 0;
26         for (int i = 1;i <= q;i++)
27         {
28             if (s[i] == t[i]) a++;
29             else d++;
30         }
31         int ans = 0;
32         if (x <= a) ans = x + d;
33         else ans = a + d - (x-a);
34         printf("%d\n",ans);
35     }
36     return 0;
37 }

View Code

7.20 >_<

手还痛.....qwq

论文被批评写得不好了....qwq

7.21 >_<

接下来要做的是尽快把论文搞成英文以及导管的定位

前几天问wkc YITU面试都是什么....几乎都是算法的题目( 而且不像是leetcode 那样的...而且听他说了说题目还都不会的样子TAT

所以感觉自己还是差很多...很多...（所以要好好加油呀！XDDDDD

C - Multi-judge Solving 这是一道没读对题意wa了5发的水题

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
typedef long long LL;
const int maxn = 1e5+5;int n,k,a[maxn];int main()
{scanf("%d %d",&n,&k);for (int i = 1;i <= n;i++) scanf("%d",&a[i]);sort(a+1,a+n+1);int ans = 0,pos;for (int i = 1;i <= n;i++){if (2 * k >= a[i]) {k = max(k,a[i]);continue;}else{while(k*2 < a[i]) k = k << 1,ans++;k = max(k,a[i]);}//printf("k = %d ans = %d\n",k,ans);
    }printf("%d\n",ans);return 0;
}

View Code

7.22 >_<

Problem C. Magical Thinking

大数据的做法，因为只有两个人，所以一共四种情况

枚举这四种情况分别是错的个数为i j k l 在满足第一个人和第二个人得分的情况下更新答案

dp的做法还没有看

 1 #include <cstdio>
 2 #include <cstring>
 3 #include <iostream>
 4 #include <algorithm>
 5 using namespace std;
 6 typedef long long LL;
 7 const int maxn = 1e5+5;
 8 
 9 int n,q,b[105];
10 char s[105][105],t[105];
11 
12 int main()
13 {
14     int T;
15     freopen("C-large-practice.in","r",stdin);
16     freopen("Cbig.txt","w",stdout);
17     scanf("%d",&T);
18     for (int kase = 1;kase <= T;kase++)
19     {
20         scanf("%d %d",&n,&q);
21         for (int i = 1;i <= n;i++) scanf("%s",s[i] + 1);
22         scanf("%s",t + 1);
23         for (int i = 1;i <=n;i++) scanf("%d",&b[i]);
24         printf("Case #%d: ",kase);
25         int ans = 0;
26         if (n == 1)
27         {
28             int a = 0,d = 0;
29             for (int i = 1;i <= q;i++)
30             {
31                 if (s[1][i] == t[i]) a++;
32                 else d++;
33             }
34             if (b[1] <= a) ans = b[1] + d;
35             else ans = a + d - (b[1] - a);
36         }
37         else
38         {
39             int A = 0,B = 0,C = 0,D = 0;
40             for (int i = 1;i <= q;i++)
41             {
42                 if (s[1][i] == s[2][i] && s[1][i] == t[i]) A++;
43                 else if (s[1][i] != s[2][i] && s[1][i] == t[i]) B++;
44                 else if (s[1][i] != s[2][i] && s[2][i] == t[i]) C++;
45                 else D++;
46             }
47 
48             //printf("A = %d B = %d C = %d D = %d\n",A,B,C,D);
49 
50             for (int i = 0;i <= A;i++)
51             {
52                 for (int j = 0;j <= B;j++)
53                 {
54                     for (int k = 0;k <= C;k++)
55                     {
56                         for (int l = 0;l <= D;l++)
57                         {
58                             if (i + j + (C - k) + (D - l) == q - b[1])
59                                 if (i + (B - j) + k + (D - l) == q - b[2])
60                                 {
61                                     ans = max(ans,q - (i + j + k + l));
62                                 //    printf("i = %d j = %d k = %d l = %d\n",i,j,k,l);
63 
64                                 }
65                         }
66                     }
67                 }
68             }
69         }
70         printf("%d\n",ans);
71     }
72     return 0;
73 }

View Code

找不到基础(适合我这种纸张)的教程 qaq

7.23 >_<

初稿完成！

LSTM的一个栗子码着要看

7.24 >_<

又是新的一周，感觉日子过得又慢又快qaq

配置theano 和 keras 和这一篇

配置 theano 感觉有点问题

import theano 不会报错

print theano.config.blas.ldflags 也不会报错

但是 theano.test() 会报错找不到 nose.config

import keras 没有问题.... 害怕qaq

7.25 >_<

卸载py2.7 和昨天搞那一坨了，然后用 Anaconda(py3.5 ) + tensorflow

环境配的有点怕了...现在看起来好像还没有什么问题的样子

Tensorflow入门教程

LSTM预测股票每日最高票价的例子码着要看

7.26 >_<

栗子

改论文，不管是语法还是内容.... 都好多问题

粒子滤波

前两天一直在搜现成的代码例子看，然后都本地跑不通，就一直在调bug....感觉这样不好

于是还是从头看起....看不懂阿QAAAQ

主要看的这个教程

mnist 实例，卷积神经网络

7.27 >_<

跑通一个例子，菜菜！开心！

7.28 >_<

py 的插值

尝试用 GBDT

7.29 >_<

迷茫

7.30 >_<

py创建数组

自己造了一个周期性的数据验证了下GBDT，可能数据太规律，还蛮准的..（总觉得不先验证下这个东西就不敢用的感觉 - ^ -

# -*- coding:UTF-8 -*-
import pandas as pd
import numpy as np
import matplotlib.pyplot as plt
import tensorflow as tf
import scipy.interpolate as itp
from scipy import interpolate
import numpy as np
import pylab as pl
from sklearn import ensemble, cross_validationx = []
y = []for i in range(0,16 * 7):x.append(i)tmp = 0;if (i % 7 == 6):tmp = 12if (i % 7 == 5):tmp = 11if (i % 7 == 4):tmp = 5if (i % 7 == 3):tmp = 4if (i % 7 == 2):tmp = 3if (i % 7 == 1):tmp = 2if (i % 7 == 0):tmp = 1y.append(tmp)print (x)
print (y)#pl.plot(x,y,"r-")
#pl.show()# 组成 训练集 和 测试集
trainx = [[0 for i in range(8)]for j in range(92 * 90)]
trainy = [[0 for i in range(1)]for j in range(92 * 90)]
testx = [[0 for i in range(8)]for j in range(93 * 90)]for i in range(7,15 * 7):trainx[i][0] = y[i-1]#print (i,i-2)trainx[i][1] = y[i-2]trainx[i][2] = y[i-3]trainx[i][3] = y[i-4]trainx[i][4] = y[i-5]trainx[i][5] = y[i-6]trainx[i][6] = y[i-7]trainx[i][7] = i % 7trainy[i][0] = y[i]testdraw = [[0 for i in range(1)]for j in range(93 * 90)]
for i in range(0,16 * 7):testdraw[i][0] = ifor i in range(7,16 * 7):testx[i][0] = y[i-1]testx[i][1] = y[i-2]testx[i][2] = y[i-3]testx[i][3] = y[i-4]testx[i][4] = y[i-5]testx[i][5] = y[i-6]testx[i][6] = y[i-7]testx[i][7] = i % 7params = {'n_estimators': 500, 'max_depth': 4, 'min_samples_split': 1,'learning_rate': 0.01, 'loss': 'ls'}
clf = ensemble.GradientBoostingRegressor(**params)clf.fit(trainx,trainy)
pre = clf.predict(testx)#pl.plot(x,y,'k--')
#pl.plot(testdraw,pre,'r--')
#pl.show()print (pre)x_ = []
y_ = []for i in range(0,16 * 7):x_.append(i)y_.append(pre[i])print (i,pre[i])pl.plot(x,y,'k--')
pl.plot(x_,y_,'r--')
pl.show()

View Code

现在能跑的一个代码

  1 # -*- coding:UTF-8 -*-
  2 import pandas as pd
  3 import numpy as np
  4 import matplotlib.pyplot as plt
  5 import tensorflow as tf
  6 import scipy.interpolate as itp
  7 from scipy import interpolate
  8 import numpy as np
  9 import pylab as pl
 10 from sklearn import ensemble, cross_validation
 11 
 12 f = open('buquan.txt')
 13 line = f.readline()
 14 
 15 X = []
 16 Y = []
 17 
 18 while line:
 19     x = ''
 20     pos = 0
 21     for i in range(0,line.__len__()):
 22         if (line[i] == ' '):
 23             pos = i
 24             break;
 25         x = x + line[i]
 26 
 27     x = int(x)
 28     X.append(x)
 29 
 30     y = ''
 31     for i in range(pos+1,line.__len__()):
 32         y = y + line[i]
 33     y = float(y)
 34     Y.append(y)
 35 
 36     line = f.readline()
 37 
 38 f.close()
 39 
 40 
 41 x_use = []
 42 y_use = []
 43 x_debug = []
 44 y_debug = []
 45 
 46 
 47 # 用五周的数据来训练 组成 训练集 测试集
 48 for i in range(0,24 * 30 * 7 * 6):
 49     x_use.append(i)
 50     #print (i)
 51     y_use.append(Y[i])
 52 
 53 #pl.plot(x_use,y_use,'r--')
 54 #pl.show()
 55 
 56 trainx = [[0 for i in range(30)]for j in range(24 * 30 * 7 * 7)]
 57 trainy = [[0 for i in range(1)]for j in range(24 * 30 * 7 * 7)]
 58 testx = [[0 for i in range(30)]for j in range(24 * 30 * 7 * 7)]
 59 
 60 for i in range(30,24 * 30 * 7 * 5):
 61     id = x_use[i] / 720 % 7 + 2
 62     trainx[i][0] = id;
 63     x_debug.append(i)
 64     y_debug.append(y_use[i])
 65     for j in range(1,30):
 66         #print (i,j,i-(30-j))
 67         trainx[i][j] = y_use[i - (30 - j)]
 68     trainy[i][0] = y_use[i]
 69 
 70 #print (trainx)
 71 #print (trainy)
 72 
 73 for i in range(30,24 * 30 * 7 * 6):
 74     id = x_use[i] / 720 % 7 + 2
 75     testx[i][0] = id;
 76     for j in range(1,30):
 77         #print (i,j,i-(30-j))
 78         testx[i][j] = y_use[i - (30 - j)]
 79 
 80 #print (testx)
 81 
 82 
 83 
 84 params = {'n_estimators': 500, 'max_depth': 4, 'min_samples_split': 1,
 85            'learning_rate': 0.01, 'loss': 'ls'}
 86 clf = ensemble.GradientBoostingRegressor(**params)
 87 
 88 clf.fit(trainx,trainy)
 89 pre = clf.predict(testx)
 90 
 91 x_ = []
 92 y_ = []
 93 
 94 for i in range(0,24 * 30 * 7 * 6):
 95     x_.append(i)
 96     y_.append(pre[i])
 97     #print (i,pre[i])
 98 
 99 pl.plot(x_,y_,'r--')
100 pl.plot(x_use,y_use,'k--')
101 pl.plot(x_debug,y_debug,'c--')
102 pl.show()