最近看的日剧，如果真有这种事，这姑娘濒死之际后悔吗？为了正义，坚持原则而死，为了更“好”得活而破坏自己的原则，哪一种人生更好呢？
鄙人忍不住假想，如果是我，我会怎么做……大概会跟这姑娘一样吧，苟活一生又有什么意思呢，倒不如为自己的追求去死，也算死得其所。不过，如果真有这种事发生在我身边，说不定我会去要挟，敲诈对方，而不是去告发啦……呵呵。
2010年4月24日星期六
正直地死和迷糊地活，哪个更好？
最近看的日剧，如果真有这种事，这姑娘濒死之际后悔吗？为了正义，坚持原则而死，为了更“好”得活而破坏自己的原则，哪一种人生更好呢？
鄙人忍不住假想，如果是我，我会怎么做……大概会跟这姑娘一样吧，苟活一生又有什么意思呢，倒不如为自己的追求去死，也算死得其所。不过，如果真有这种事发生在我身边，说不定我会去要挟，敲诈对方，而不是去告发啦……呵呵。
2010年4月20日星期二
喜感的学校，喜感的食堂，喜感的浴室
我校以收费之贵在大陆高校里夺得探花，而且这把交椅越坐越稳——
http://www.jnrain.com/wForum/disparticle.php?boardName=JNUInfo&ID=1145&pos=10
http://www.jnrain.com/wForum/disparticle.php?boardName=JNUInfo&ID=1162&pos=9
http://www.jnrain.com/wForum/disparticle.php?boardName=JNUInfo&ID=1148&pos=7
关于我校学生浴室热水水质问题的情况通报（二）
在4月15日、16日无锡电视台新闻综合频道《第一看点》连续报道有关我校学生浴室热水水质问题的新闻后，学校高度重视，于4月17日上午9点紧急召开突发公共卫生事件协调会，并于第一时间邀请无锡市卫生局、市疾病预防控制中心、市卫生监督所对我校学生浴室热水供水商——旺盛公司所供热水质量进行全面调查检测。
17日上午10点，无锡市卫生局立即组织市卫生监督所卫生监督员和市疾病预防控制中心技术人员对事件进行现场核实调查，并实地采集样品送疾控中心检测。在现场，市卫生部门工作人员对我校南北区学生浴室淋浴水、储热水箱水和供应商水源出厂水、热水槽罐车出罐水及自来水水表前的自来水等6件水样进行监督采样，并按照国家《生活饮用水卫生标准》（GB5749-2006）送无锡市疾病控制中心进行检测。根据无锡市卫生部门的意见，在尚未出具新的水质检测报告之前，学校决定暂停开放学生浴室，并在食堂、学生浴室等场所和校园网络上予以公告。
为保障师生身体健康，防止负面报道和媒体炒作引发上海世博会期间的不稳定事端，维护校园稳定，17日下午，学校向市委市政府紧急报送《XX大学关于浴室用水质量问题的报告》，提请市委市政府协助我校做好相关工作。
18日上午10：30，无锡市政府召集无锡市委宣传部、无锡市卫生局、市疾病预防控制中心、市卫生监督所、无锡市广电局、无锡市公安局、无锡市信访局、xx大学、无锡电视台等单位召开紧急会议。会上，市政府通报了事件的情况和最新的质检结果：市疾病预防控制中心对6件水样的9项重金属指标进行了全面检测，先期检测的铅等三项指标检测结果合格。同时，市政府要求市委宣传部密切关注广电、报业和网络媒体对此次事件的宣传报道，加强舆论监管，坚决防止出现媒体炒作事件；要求无锡电视台迅速澄清事实、公布权威结果，以正视听；要求公安部门迅速查清《第一看点》15、16日报道的事件真相；要求卫生防疫部门加强检测并尽快出具后续报告，同时加强对相关热水供应点的检查，以保障市民和我校学生身体健康。根据检测结果和市政府会议要求，学校决定恢复浴室开放，并在第一时间发布了相关信息。
18日下午14点，学校再次召开突发公共卫生事件通气会，向各学院学生工作副书记、分团委书记和学生会组织的负责人通报此事件的调查情况，听取师生意见建议。会上，无锡市卫生监督所的负责同志通报了初期调查结果：1.XX大学学生饮用热水是市政自来水经三台燃气热水炉加热后供应给学生，与浴室用热水完全是不同的供水系统。浴室用热水由无锡市旺盛公司提供；2.无锡市旺盛公司有两个热水加工点，一个采自无锡强盛钢铁厂的蒸汽，主要供应市内小型浴室；另一个采自中外合作无锡能达热电有限公司的蒸汽，供应XX大学和市职业技术学院。3. 市疾控中心对17日现场采集的6件水样的铅等三项指标检测结果合格。
下一步，学校将严格按照市委市政府要求，密切关注事态发展，并配合市卫生部门、公安部门、宣传部门做好相关工作，第一时间向师生公布事件调查结果和工作进展，请师生们关注并予以监督。同时，校领导希望学校后勤等部门进一步加强和改善管理，提高服务质量；各学院及相关部门也要做好进一步的解释说明和宣传引导，消除师生疑虑和担忧，全力维护校园稳定。
党委宣传部
2010年4月18日
斜体部分亮点，呵呵
目前，学校bbs刚刚重新对外开放，所有关于洗澡水的帖子全部被“河蟹”，只能集中在一个“集中讨论帖里”……除了类似“无耻媒体，纯属谣言的制造者”的五毛贴，现在五毛们已经被同学们形成围观之势……喜感啊。
另外，食堂的米饭终于被曝光出用了陈化粮，山人自己觉得吧，那食堂到不真一定用了陈化粮，但是地沟油之类的肯定没少用，但是没人信学校，这个公信力就好像婚前不当行为一样，没了就是没了……餐盘恶心的跟狗盘似地，地板恶心的跟盥洗室似地，价钱恶心的跟打劫的似地，还说是为学生服务，每年赔本经营多少RMB，那校门外的一溜小饭馆岂不是学雷锋倒贴我们学生。
不过我们学校脸皮不是一般的厚，该不承认肯定不会承认，好像承认了就丧失权威马上颜面扫地好像没穿衣服的皇帝承认了自己确实是裸奔一样……我看他什么时候嫌冷。
http://www.jnrain.com/wForum/disparticle.php?boardName=JNUInfo&ID=1145&pos=10
http://www.jnrain.com/wForum/disparticle.php?boardName=JNUInfo&ID=1162&pos=9
http://www.jnrain.com/wForum/disparticle.php?boardName=JNUInfo&ID=1148&pos=7
关于我校学生浴室热水水质问题的情况通报（二）
在4月15日、16日无锡电视台新闻综合频道《第一看点》连续报道有关我校学生浴室热水水质问题的新闻后，学校高度重视，于4月17日上午9点紧急召开突发公共卫生事件协调会，并于第一时间邀请无锡市卫生局、市疾病预防控制中心、市卫生监督所对我校学生浴室热水供水商——旺盛公司所供热水质量进行全面调查检测。
17日上午10点，无锡市卫生局立即组织市卫生监督所卫生监督员和市疾病预防控制中心技术人员对事件进行现场核实调查，并实地采集样品送疾控中心检测。在现场，市卫生部门工作人员对我校南北区学生浴室淋浴水、储热水箱水和供应商水源出厂水、热水槽罐车出罐水及自来水水表前的自来水等6件水样进行监督采样，并按照国家《生活饮用水卫生标准》（GB5749-2006）送无锡市疾病控制中心进行检测。根据无锡市卫生部门的意见，在尚未出具新的水质检测报告之前，学校决定暂停开放学生浴室，并在食堂、学生浴室等场所和校园网络上予以公告。
为保障师生身体健康，防止负面报道和媒体炒作引发上海世博会期间的不稳定事端，维护校园稳定，17日下午，学校向市委市政府紧急报送《XX大学关于浴室用水质量问题的报告》，提请市委市政府协助我校做好相关工作。
18日上午10：30，无锡市政府召集无锡市委宣传部、无锡市卫生局、市疾病预防控制中心、市卫生监督所、无锡市广电局、无锡市公安局、无锡市信访局、xx大学、无锡电视台等单位召开紧急会议。会上，市政府通报了事件的情况和最新的质检结果：市疾病预防控制中心对6件水样的9项重金属指标进行了全面检测，先期检测的铅等三项指标检测结果合格。同时，市政府要求市委宣传部密切关注广电、报业和网络媒体对此次事件的宣传报道，加强舆论监管，坚决防止出现媒体炒作事件；要求无锡电视台迅速澄清事实、公布权威结果，以正视听；要求公安部门迅速查清《第一看点》15、16日报道的事件真相；要求卫生防疫部门加强检测并尽快出具后续报告，同时加强对相关热水供应点的检查，以保障市民和我校学生身体健康。根据检测结果和市政府会议要求，学校决定恢复浴室开放，并在第一时间发布了相关信息。
18日下午14点，学校再次召开突发公共卫生事件通气会，向各学院学生工作副书记、分团委书记和学生会组织的负责人通报此事件的调查情况，听取师生意见建议。会上，无锡市卫生监督所的负责同志通报了初期调查结果：1.XX大学学生饮用热水是市政自来水经三台燃气热水炉加热后供应给学生，与浴室用热水完全是不同的供水系统。浴室用热水由无锡市旺盛公司提供；2.无锡市旺盛公司有两个热水加工点，一个采自无锡强盛钢铁厂的蒸汽，主要供应市内小型浴室；另一个采自中外合作无锡能达热电有限公司的蒸汽，供应XX大学和市职业技术学院。3. 市疾控中心对17日现场采集的6件水样的铅等三项指标检测结果合格。
下一步，学校将严格按照市委市政府要求，密切关注事态发展，并配合市卫生部门、公安部门、宣传部门做好相关工作，第一时间向师生公布事件调查结果和工作进展，请师生们关注并予以监督。同时，校领导希望学校后勤等部门进一步加强和改善管理，提高服务质量；各学院及相关部门也要做好进一步的解释说明和宣传引导，消除师生疑虑和担忧，全力维护校园稳定。
党委宣传部
2010年4月18日
斜体部分亮点，呵呵
目前，学校bbs刚刚重新对外开放，所有关于洗澡水的帖子全部被“河蟹”，只能集中在一个“集中讨论帖里”……除了类似“无耻媒体，纯属谣言的制造者”的五毛贴，现在五毛们已经被同学们形成围观之势……喜感啊。
另外，食堂的米饭终于被曝光出用了陈化粮，山人自己觉得吧，那食堂到不真一定用了陈化粮，但是地沟油之类的肯定没少用，但是没人信学校，这个公信力就好像婚前不当行为一样，没了就是没了……餐盘恶心的跟狗盘似地，地板恶心的跟盥洗室似地，价钱恶心的跟打劫的似地，还说是为学生服务，每年赔本经营多少RMB，那校门外的一溜小饭馆岂不是学雷锋倒贴我们学生。
不过我们学校脸皮不是一般的厚，该不承认肯定不会承认，好像承认了就丧失权威马上颜面扫地好像没穿衣服的皇帝承认了自己确实是裸奔一样……我看他什么时候嫌冷。
算法设计与分析（郑宗汉）几个代码错误
上学期用了这本书。而且还有悲剧的上机，无奈这本书的代码错误有点多，小到括号丢失，大到原则错误……毕竟是本教材，小女子就把自己上机的时候发现的几个错误贴上来，也为了方便后来的师弟师妹（抄代码？呵呵）学习。
多项式乘法分治算法和动态规划有单独写出来，最后是我们的结课设计，堆排序，基数排序，合并有序数列和冒泡都在里面，如果有人需要单独的函数，我相信从这里面摘出来比自己改要容易，呵呵，当然是运行没有问题的。
多项式乘法分治算法，貌似只有这本书用这个名称
#include stdio.h
#include stdlib.h
#include math.h
void product(float p[],float q[],float c[])
{
float a=p[0];
float b=p[1];
float e=q[0];
float d=q[1];
c[0]=a*e;
c[2]=b*d;
c[1]=(a+b)*(e+d)-c[0]-c[2];
}
void plus(float p[],float q[],float c[],int n)
{
int i;
for (i=0;i
#define MAX_TYPE 1000
#define ZERO_TYPE 0
template
struct NODE
{
int v_num;
Type len;
NODE *next;
};
template
Type fgraph(NODE node[],int route[],int n)
{
int i;
NODE *pnode;
int *path = new int[n];
Type min_cost,*cost = new Type[n];
for(i=0;i=0;i--)
{
pnode=&node[i];
while(pnode!=NULL)
{
if ((pnode->len+cost[pnode->v_num])len+cost[pnode->v_num];
path[i]=pnode->v_num;
}
pnode=pnode->next;
}
}
i=0;
while((route[i]!=n-1)&&(path[i]!=-1))
{
i++;
route[i]=path[route[i-1]];
}
min_cost=cost[0];
delete path;
delete cost;
return min_cost;
}
int main()
{
int i,rep,p,n=10;
int route[100];
NODE node[30];
int A[100]={1,30,2,20,3,25,100,4,32,5,35,6,40,100,4,30,5,28,6,10,100,4,35,5,25,6,20,100,7,21,8,23,100,7,20,8,38,100,7,29,8,40,100,9,38,100,9,28,100};
p=0;
for(i=0;i *qnode,*snode;
snode=qnode=&node[i];
while(1)
{
rep=A[p++];
if(rep>n)
{
snode->next=NULL;
delete qnode;
break;
}
else
{
qnode->v_num=rep;
qnode->len=A[p++];
qnode->next=new NODE;
snode=qnode;
qnode=qnode->next;
}
}
}
printf("%d\n",fgraph(node,route,n));
for(i=0;i<5;i++)>
#include math.h
#include stdlib.h
#include time.h
#include windows.h
void swap(int &x,int &y)
{
int temp;
temp=x;
x=y;
y=temp;
}
void Bubble(int a[],int n)
{
int i,k;
for(k=n-1;k>0;k--)
{
for(i=0;ia[i+1])
{
swap(a[i],a[i+1]);
}
}
}
}
void Merge(int a[],int p,int q,int r,int m)
{
int *bp=new int[m];
int i,j,k;
i=p;j=q+1;k=0;
while(i<=q&&j<=r) { if(a[i]<=a[j]) bp[k++]=a[i++]; else bp[k++]=a[j++]; } if(i==q+1) for(;j<=r;j++) bp[k++]=a[j]; else for(;i<=q;i++) bp[k++]=a[i]; k=0; for(i=p;i<=r;i++) a[i]=bp[k++]; delete bp; } void theBubble() { int i,q,s[10001],n=10000; FILE *fp1,*fp2; fp1=fopen("cindy.txt","r"); fp2=fopen("out.txt","w"); for(i=0;i<100;i++) q="0;q<100;q++)" t =" clock();" sec =" t;" sec2 =" clock();" p="L-">next;
if (p!=L)
{
p->prior->next=p->next;
p->next->prior=p->prior;
}
else p=NULL;
return p;
}
void add_entry(LinkedList *L,LinkedList *p){
p->prior=L->prior;
p->next=L;
L->prior->next=p;
L->prior=p;
}
int get_digital(LinkedList *p,int i){
int key ;
key=p->key;
if (i!=0)
key =key /(int)pow(10,i);
return key %10;
}
void append (LinkedList *L,LinkedList *L1){
if(L1->next!=L1){
L->prior->next=L1->next;
L1->next->prior=L->prior;
L1->prior->next=L;
L->prior=L1->prior;
}
}
void radix_sort(LinkedList *L,int k){
LinkedList *Lhead[10],*p;
int i,j;
for (i=0;i<10;i++) i="0;inext=Lhead[j];
Lhead[j]->prior=Lhead[j];
}
while (L->next!=L){
p=del_entry(L);
j=get_digital(p,i);
add_entry(Lhead[j],p);
}
for (j=0;j<10;j++) i="0;i<10;i++)" s="-1649273;" fp1="fopen(" fp2="fopen(" n="10000;" tp2 =" (LinkedList">key=-1;
tp2->prior=NULL;
tp2->next=NULL;
tp1=tp2;
l=tp1;
for(i=0;ikey=j;
tp1->next=tp2;
tp2->prior=tp1;
tp2->next=NULL;
tp1=tp2;
}
k=1;
for(i=0;k==1;i++){
if((s/(int)pow(10,i))==0){k=0;}
}
k=i-1;
tp2->next=l;
clock_t t = clock();
long sec = t;
radix_sort(l,k);
tp1=l->next;
long sec2 = clock();
printf("%ld\n",sec2-sec);
fclose(fp1);
fclose(fp2);
}
/*****the radix******/
template
void swap(Type &x,Type &y)
{
Type temp;
temp=x;
x=y;
y=temp;
}
template
void sift_down(Type H[],int n,int i)
{
BOOL done=FALSE;
if ((2*i)<=n) { while(!done&&(i=2*i)<=n) { if (i+1<=n && H[i+1]>H[i])
i=i+1;
if (H[i/2]
void make_heap(Type A[],int n)
{
int i;
A[n]=A[0];
for(i=n/2;i>=1;i--)
sift_down(A,n,i);
}
template
void heap_sort(Type A[],int n)
{
int i;
make_heap(A,n);
for (i=n;i>1;i--)
{
swap(A[1],A[i]);
sift_down(A,i-1,1);
}
}
void theheap()
{
int A[10001];
int n=10000,i,q;
FILE *fp1,*fp2;
fp1=fopen("cindy.txt","r");
fp2=fopen("out.txt","w");
for(i=0;i<100;i++) q="0;q<100;q++)" t =" clock();" sec =" t;" sec2 =" clock();">
void merge(Type a[],int p,int q,int r,int m)
{
int *bp=new int[m];
int i,j,k;
i=p;j=q+1;k=0;
while(i<=q&&j<=r) { if(a[i]<=a[j]) bp[k++]=a[i++]; else bp[k++]=a[j++]; } if(i==q+1) for(;j<=r;j++) bp[k++]=a[j]; else for(;i<=q;i++) bp[k++]=a[i]; k=0; for(i=p;i<=r;i++) a[i]=bp[k++]; delete bp; } template
void merge_sort(Type A[],int n)
{
int x,s,t=1;
while(t
s=t,t=2*s;x=0;
while(x+t
merge(A,x,x+s-1,x+t-1,t);
x=x+t;
}
if(x+s
merge(A,x,x+s-1,n-1,n-x);
}
}
void themerge()
{
long int a,b;
int i,q,n=10000;
int A[10001];
FILE *fp1,*fp2;
fp1=fopen("cindy.txt","r");
fp2=fopen("out.txt","w");
for(i=0;i<100;i++)
{
for(q=0;q<100;q++)
{
fscanf(fp1,"%d",&A[i]);
}
}
a=clock();
merge_sort(A,n);
b=clock();
printf("%ld\n",b-a);
fclose(fp1);
fclose(fp2);
}
/************************/
void main()
{
FILE *fp;
fp=fopen("cindy.txt","wr");
int m,p,ch1;
double ch;
char rep=getchar();
if(rep=='f')
{
for(m=0;m<100;m++)
{
for(p=0;p<100;p++)
{
ch=rand()%1000/10.0;
fprintf(fp,"%f",ch);
}
}
printf("DATE TYPE IS FLOAT\n");
}
if(rep=='i')
{
for(m=0;m<100;m++)
{
for(p=0;p<100;p++)
{
ch1=rand()%1000;
fprintf(fp,"%d",ch1);
}
}
printf("DATE TYPE IS INT\n");
}
theBubble();
themerge();
theheap();
theradix();
fclose(fp);
}
多项式乘法分治算法和动态规划有单独写出来，最后是我们的结课设计，堆排序，基数排序，合并有序数列和冒泡都在里面，如果有人需要单独的函数，我相信从这里面摘出来比自己改要容易，呵呵，当然是运行没有问题的。
多项式乘法分治算法，貌似只有这本书用这个名称
#include stdio.h
#include stdlib.h
#include math.h
void product(float p[],float q[],float c[])
{
float a=p[0];
float b=p[1];
float e=q[0];
float d=q[1];
c[0]=a*e;
c[2]=b*d;
c[1]=(a+b)*(e+d)-c[0]-c[2];
}
void plus(float p[],float q[],float c[],int n)
{
int i;
for (i=0;i
#define MAX_TYPE 1000
#define ZERO_TYPE 0
template
struct NODE
{
int v_num;
Type len;
NODE
};
template
Type fgraph(NODE
{
int i;
NODE
int *path = new int[n];
Type min_cost,*cost = new Type[n];
for(i=0;i
{
pnode=&node[i];
while(pnode!=NULL)
{
if ((pnode->len+cost[pnode->v_num])
path[i]=pnode->v_num;
}
pnode=pnode->next;
}
}
i=0;
while((route[i]!=n-1)&&(path[i]!=-1))
{
i++;
route[i]=path[route[i-1]];
}
min_cost=cost[0];
delete path;
delete cost;
return min_cost;
}
int main()
{
int i,rep,p,n=10;
int route[100];
NODE
int A[100]={1,30,2,20,3,25,100,4,32,5,35,6,40,100,4,30,5,28,6,10,100,4,35,5,25,6,20,100,7,21,8,23,100,7,20,8,38,100,7,29,8,40,100,9,38,100,9,28,100};
p=0;
for(i=0;i
snode=qnode=&node[i];
while(1)
{
rep=A[p++];
if(rep>n)
{
snode->next=NULL;
delete qnode;
break;
}
else
{
qnode->v_num=rep;
qnode->len=A[p++];
qnode->next=new NODE
snode=qnode;
qnode=qnode->next;
}
}
}
printf("%d\n",fgraph(node,route,n));
for(i=0;i<5;i++)>
#include math.h
#include stdlib.h
#include time.h
#include windows.h
void swap(int &x,int &y)
{
int temp;
temp=x;
x=y;
y=temp;
}
void Bubble(int a[],int n)
{
int i,k;
for(k=n-1;k>0;k--)
{
for(i=0;i
{
swap(a[i],a[i+1]);
}
}
}
}
void Merge(int a[],int p,int q,int r,int m)
{
int *bp=new int[m];
int i,j,k;
i=p;j=q+1;k=0;
while(i<=q&&j<=r) { if(a[i]<=a[j]) bp[k++]=a[i++]; else bp[k++]=a[j++]; } if(i==q+1) for(;j<=r;j++) bp[k++]=a[j]; else for(;i<=q;i++) bp[k++]=a[i]; k=0; for(i=p;i<=r;i++) a[i]=bp[k++]; delete bp; } void theBubble() { int i,q,s[10001],n=10000; FILE *fp1,*fp2; fp1=fopen("cindy.txt","r"); fp2=fopen("out.txt","w"); for(i=0;i<100;i++) q="0;q<100;q++)" t =" clock();" sec =" t;" sec2 =" clock();" p="L-">next;
if (p!=L)
{
p->prior->next=p->next;
p->next->prior=p->prior;
}
else p=NULL;
return p;
}
void add_entry(LinkedList *L,LinkedList *p){
p->prior=L->prior;
p->next=L;
L->prior->next=p;
L->prior=p;
}
int get_digital(LinkedList *p,int i){
int key ;
key=p->key;
if (i!=0)
key =key /(int)pow(10,i);
return key %10;
}
void append (LinkedList *L,LinkedList *L1){
if(L1->next!=L1){
L->prior->next=L1->next;
L1->next->prior=L->prior;
L1->prior->next=L;
L->prior=L1->prior;
}
}
void radix_sort(LinkedList *L,int k){
LinkedList *Lhead[10],*p;
int i,j;
for (i=0;i<10;i++) i="0;i
Lhead[j]->prior=Lhead[j];
}
while (L->next!=L){
p=del_entry(L);
j=get_digital(p,i);
add_entry(Lhead[j],p);
}
for (j=0;j<10;j++) i="0;i<10;i++)" s="-1649273;" fp1="fopen(" fp2="fopen(" n="10000;" tp2 =" (LinkedList">key=-1;
tp2->prior=NULL;
tp2->next=NULL;
tp1=tp2;
l=tp1;
for(i=0;i
tp1->next=tp2;
tp2->prior=tp1;
tp2->next=NULL;
tp1=tp2;
}
k=1;
for(i=0;k==1;i++){
if((s/(int)pow(10,i))==0){k=0;}
}
k=i-1;
tp2->next=l;
clock_t t = clock();
long sec = t;
radix_sort(l,k);
tp1=l->next;
long sec2 = clock();
printf("%ld\n",sec2-sec);
fclose(fp1);
fclose(fp2);
}
/*****the radix******/
template
void swap(Type &x,Type &y)
{
Type temp;
temp=x;
x=y;
y=temp;
}
template
void sift_down(Type H[],int n,int i)
{
BOOL done=FALSE;
if ((2*i)<=n) { while(!done&&(i=2*i)<=n) { if (i+1<=n && H[i+1]>H[i])
i=i+1;
if (H[i/2]
void make_heap(Type A[],int n)
{
int i;
A[n]=A[0];
for(i=n/2;i>=1;i--)
sift_down(A,n,i);
}
template
void heap_sort(Type A[],int n)
{
int i;
make_heap(A,n);
for (i=n;i>1;i--)
{
swap(A[1],A[i]);
sift_down(A,i-1,1);
}
}
void theheap()
{
int A[10001];
int n=10000,i,q;
FILE *fp1,*fp2;
fp1=fopen("cindy.txt","r");
fp2=fopen("out.txt","w");
for(i=0;i<100;i++) q="0;q<100;q++)" t =" clock();" sec =" t;" sec2 =" clock();">
void merge(Type a[],int p,int q,int r,int m)
{
int *bp=new int[m];
int i,j,k;
i=p;j=q+1;k=0;
while(i<=q&&j<=r) { if(a[i]<=a[j]) bp[k++]=a[i++]; else bp[k++]=a[j++]; } if(i==q+1) for(;j<=r;j++) bp[k++]=a[j]; else for(;i<=q;i++) bp[k++]=a[i]; k=0; for(i=p;i<=r;i++) a[i]=bp[k++]; delete bp; } template
void merge_sort(Type A[],int n)
{
int x,s,t=1;
while(t
s=t,t=2*s;x=0;
while(x+t
merge(A,x,x+s-1,x+t-1,t);
x=x+t;
}
if(x+s
merge(A,x,x+s-1,n-1,n-x);
}
}
void themerge()
{
long int a,b;
int i,q,n=10000;
int A[10001];
FILE *fp1,*fp2;
fp1=fopen("cindy.txt","r");
fp2=fopen("out.txt","w");
for(i=0;i<100;i++)
{
for(q=0;q<100;q++)
{
fscanf(fp1,"%d",&A[i]);
}
}
a=clock();
merge_sort(A,n);
b=clock();
printf("%ld\n",b-a);
fclose(fp1);
fclose(fp2);
}
/************************/
void main()
{
FILE *fp;
fp=fopen("cindy.txt","wr");
int m,p,ch1;
double ch;
char rep=getchar();
if(rep=='f')
{
for(m=0;m<100;m++)
{
for(p=0;p<100;p++)
{
ch=rand()%1000/10.0;
fprintf(fp,"%f",ch);
}
}
printf("DATE TYPE IS FLOAT\n");
}
if(rep=='i')
{
for(m=0;m<100;m++)
{
for(p=0;p<100;p++)
{
ch1=rand()%1000;
fprintf(fp,"%d",ch1);
}
}
printf("DATE TYPE IS INT\n");
}
theBubble();
themerge();
theheap();
theradix();
fclose(fp);
}
可爱的青衣
上周五的信息隐藏作业，要求要一副bmp隐写到另一幅bmp里，于是找到了这两幅图——很可爱的红色小青衣。
隐写之后，变成了——锵锵锵，更可爱的粉红色小青衣 lol 讲代码偷偷附上，可matlab小自己随便玩玩。（此代码来自信息隐藏实验教程源代码，虽然大陆连不到这里来，还是多嘴说一句，请尊重原作者，勿用于商业或违法行为^_^）
%文件名：imagehide.m
%程序员：王*霞*仙
%编写时间：2004.2.5
%函数功能：直接将秘密图像的高4bit隐藏在RGB载体图像的R、G、B层中所选的那一层的低4bit,并将秘密图像提取出来，最后显示。要求载体图像的大小大于等于秘密图像的大小。且秘密图像是二值或灰度图像。
%输入格式：
%data=imagehide('c:\Lenna.bmp','c:\woman.bmp','c:\mix.bmp','bmp',3)
%参数说明：
% cover是载体图像的地址
% massage是秘密图像的地址
% goleimage是隐藏后图像的地址
% pemission是图像的类型
%level是作为的载体的具体层。R为1，G为2，B为3。
%data是隐藏后图像的矩阵
function data=imagehide(cover,massage,goleimage,permission,level)
%提取图像信息并分层
cover = imread(cover,permission);
data=cover;
msg=imread(massage,permission);
[row,col]=size(cover);
cover1=cover(:,:,level);
%置载体图像R层的低4bit为0
for i=1:row
for j=1:col/3
cover1(i,j)=bitand(cover1(i,j),240);
end
end
%置秘密图像的低4bit为0
takemsg4=bitand(msg,240);
%将秘密图像的高4bit右移4位
shiftmsg4=bitshift(takemsg4,-4);
%图像隐藏
for i=1:row
for j=1:col/3
cover1(i,j)=bitor(cover1(i,j),shiftmsg4(i,j));
end
end
%写回并保存
data(:,:,level)=cover1;
imwrite(data,goleimage,permission);
%提取秘密图像信息,检测隐藏效果
data=imread(goleimage,permission);
[row,col]=size(data);
A=data(:,:,level);
for i=1:row
for j=1:col/3
A(i,j)=bitand(A(i,j),15);
end
end
A=bitshift(A,4);
%显示结果
subplot(221),imshow(cover);title('载体图像');
subplot(222),imshow(massage);title('秘密图像');
subplot(223),imshow(data);title('隐藏后的图像');
subplot(224),imshow(A);title('提取的秘密图像');
隐写之后，变成了——锵锵锵，更可爱的粉红色小青衣 lol 讲代码偷偷附上，可matlab小自己随便玩玩。（此代码来自信息隐藏实验教程源代码，虽然大陆连不到这里来，还是多嘴说一句，请尊重原作者，勿用于商业或违法行为^_^）
%文件名：imagehide.m
%程序员：王*霞*仙
%编写时间：2004.2.5
%函数功能：直接将秘密图像的高4bit隐藏在RGB载体图像的R、G、B层中所选的那一层的低4bit,并将秘密图像提取出来，最后显示。要求载体图像的大小大于等于秘密图像的大小。且秘密图像是二值或灰度图像。
%输入格式：
%data=imagehide('c:\Lenna.bmp','c:\woman.bmp','c:\mix.bmp','bmp',3)
%参数说明：
% cover是载体图像的地址
% massage是秘密图像的地址
% goleimage是隐藏后图像的地址
% pemission是图像的类型
%level是作为的载体的具体层。R为1，G为2，B为3。
%data是隐藏后图像的矩阵
function data=imagehide(cover,massage,goleimage,permission,level)
%提取图像信息并分层
cover = imread(cover,permission);
data=cover;
msg=imread(massage,permission);
[row,col]=size(cover);
cover1=cover(:,:,level);
%置载体图像R层的低4bit为0
for i=1:row
for j=1:col/3
cover1(i,j)=bitand(cover1(i,j),240);
end
end
%置秘密图像的低4bit为0
takemsg4=bitand(msg,240);
%将秘密图像的高4bit右移4位
shiftmsg4=bitshift(takemsg4,-4);
%图像隐藏
for i=1:row
for j=1:col/3
cover1(i,j)=bitor(cover1(i,j),shiftmsg4(i,j));
end
end
%写回并保存
data(:,:,level)=cover1;
imwrite(data,goleimage,permission);
%提取秘密图像信息,检测隐藏效果
data=imread(goleimage,permission);
[row,col]=size(data);
A=data(:,:,level);
for i=1:row
for j=1:col/3
A(i,j)=bitand(A(i,j),15);
end
end
A=bitshift(A,4);
%显示结果
subplot(221),imshow(cover);title('载体图像');
subplot(222),imshow(massage);title('秘密图像');
subplot(223),imshow(data);title('隐藏后的图像');
subplot(224),imshow(A);title('提取的秘密图像');
订阅：
评论 (Atom)