C语言问题

Moyu18_06_12 2014-05-14 11:13:09
#include<stdio.h>
#include<math.h>
#define E 0.00001
typedef struct {
double x;
double y;
}point;
typedef struct{
point p1;
point p2;
}center;
double distance(point a,point b);
center cluster(point* p,int i);
point cen(point* p,int i,int t);
center part(point* p,int i,point p1,point p2);
int main()
{
point arr[4]={{1,2},{2,4},{2,7},{3,8}};
center c=cluster(arr,4);
printf("%f,%f\n%f,%f",c.p1.x ,c.p1.y,c.p2.x,c.p2.y);
return 0;
}
center cluster(point* p,int i)
{
point a, b,a0,b0;
int n=(i*(i-1))/2;
double max=0;
int x=0;
double d[n];
double e1,e2;
center c0;

for(int m=0;m<i-1;m++)
for(int n=1;n<i;n++)
{
d[x]=distance(p[m],p[n]);
if(d[x]>max)
{

max=d[x];
a=p[m];// 把距离最远的点付给a,b
b=p[n];
}
x++;
}
do
{
a0=a;
b0=b;
c0=part(p,i,a,b);
a=c0.p1 ;
b=c0.p2;
e1=distance(a0,a);
e2=distance(b0,b);
}while((e1>E)|(e2>E));
center c={a,b};
return c;
}
//求两点之间距离函数
double distance(point a,point b)
{
double dte;
dte=sqrt((b.x-a.x)*(b.x-a.x)+(b.y-a.y)*(b.y-a.y));
return dte;
}
//聚为两类
center part(point* p,int i,point p1,point p2)
{
point arr1[i];
point arr2[i];
int t1=0;
int t2=0;
for(int n=0;n<i;n++)
{
if(distance(p1,p[n]) < distance(p2,p[n]))
{
arr1[n]=p[n];
t1++;
}

else
{
arr2[n]=p[n];
t2++;
}
}
p1=cen(arr1,i,t1);
p2=cen(arr2,i,t2);
center c={p1,p2};
return c;
}
//求中心点函数
point cen(point* p,int i,int t)
{
point pt;
double X=0;


double Y=0;
for(int n=0;n<i;n++)
{
X+=p[n].x;
Y+=p[n].y;
}
pt.x =X/t;
pt.y =Y/t;
return pt;
}
这个程序是聚类的一个算法，但不知为什么，调试结果正确，运行就得到乱七八糟的结果，用的IDE是DEV-c++，求高手指点。。。
...全文
187 7 打赏 收藏 转发到动态 举报
写回复
用AI写文章
7 条回复
切换为时间正序
请发表友善的回复…
发表回复
Moyu18_06_12 2014-05-14
  • 打赏
  • 举报
 回复
回复4楼：改了逻辑运算符后还不行，不懂你在后面说的。。。。。。求指点
赵4老师 2014-05-14
  • 打赏
  • 举报
 回复
引用 3 楼 u012514171 的回复:
怎么感觉还是一样的乱，你手打的吧？
astyle自动整理c,c++,java源代码格式命令行工具源代码。 用astyle --help查看使用帮助。 http://download.csdn.net/detail/zhao4zhong1/5605415
赵4老师 2014-05-14
  • 打赏
  • 举报
 回复
}while((e1>E)|(e2>E)); 似应改为 }while((e1>E)||(e2>E)); 有时不将“调用函数名字＋各参数值，进入函数后各参数值，中间变量值，退出函数前准备返回的值，返回函数到调用处后函数名字＋各参数值＋返回值”这些信息写日志到文件中是无论如何也发现不了问题在哪里的，包括捕获各种异常、写日志到屏幕、单步或设断点或生成core文件、……这些方法都不行！ 写日志到文件参考下面： 
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#ifdef WIN32
    #include <windows.h>
    #include <io.h>
#else
    #include <unistd.h>
    #include <sys/time.h>
    #include <pthread.h>
    #define  CRITICAL_SECTION   pthread_mutex_t
    #define  _vsnprintf         vsnprintf
#endif
//Log{
#define MAXLOGSIZE 20000000
#define MAXLINSIZE 16000
#include <time.h>
#include <sys/timeb.h>
#include <stdarg.h>
char logfilename1[]="MyLog1.log";
char logfilename2[]="MyLog2.log";
static char logstr[MAXLINSIZE+1];
char datestr[16];
char timestr[16];
char mss[4];
CRITICAL_SECTION cs_log;
FILE *flog;
#ifdef WIN32
void Lock(CRITICAL_SECTION *l) {
    EnterCriticalSection(l);
}
void Unlock(CRITICAL_SECTION *l) {
    LeaveCriticalSection(l);
}
#else
void Lock(CRITICAL_SECTION *l) {
    pthread_mutex_lock(l);
}
void Unlock(CRITICAL_SECTION *l) {
    pthread_mutex_unlock(l);
}
#endif
void LogV(const char *pszFmt,va_list argp) {
    struct tm *now;
    struct timeb tb;

    if (NULL==pszFmt||0==pszFmt[0]) return;
    _vsnprintf(logstr,MAXLINSIZE,pszFmt,argp);
    ftime(&tb);
    now=localtime(&tb.time);
    sprintf(datestr,"%04d-%02d-%02d",now->tm_year+1900,now->tm_mon+1,now->tm_mday);
    sprintf(timestr,"%02d:%02d:%02d",now->tm_hour     ,now->tm_min  ,now->tm_sec );
    sprintf(mss,"%03d",tb.millitm);
    printf("%s %s.%s %s",datestr,timestr,mss,logstr);
    flog=fopen(logfilename1,"a");
    if (NULL!=flog) {
        fprintf(flog,"%s %s.%s %s",datestr,timestr,mss,logstr);
        if (ftell(flog)>MAXLOGSIZE) {
            fclose(flog);
            if (rename(logfilename1,logfilename2)) {
                remove(logfilename2);
                rename(logfilename1,logfilename2);
            }
        } else {
            fclose(flog);
        }
    }
}
void Log(const char *pszFmt,...) {
    va_list argp;

    Lock(&cs_log);
    va_start(argp,pszFmt);
    LogV(pszFmt,argp);
    va_end(argp);
    Unlock(&cs_log);
}
//Log}
int main(int argc,char * argv[]) {
    int i;
#ifdef WIN32
    InitializeCriticalSection(&cs_log);
#else
    pthread_mutex_init(&cs_log,NULL);
#endif
    for (i=0;i<10000;i++) {
        Log("This is a Log %04d from FILE:%s LINE:%d\n",i, __FILE__, __LINE__);
    }
#ifdef WIN32
    DeleteCriticalSection(&cs_log);
#else
    pthread_mutex_destroy(&cs_log);
#endif
    return 0;
}
//1-78行添加到你带main的.c或.cpp的那个文件的最前面
//81-85行添加到你的main函数开头
//89-93行添加到你的main函数结束前
//在要写LOG的地方仿照第87行的写法写LOG到文件MyLog1.log中
lanlvy 2014-05-14
  • 打赏
  • 举报
 回复
怎么感觉还是一样的乱，你手打的吧？
lanlvy 2014-05-14
  • 打赏
  • 举报
 回复
这样看太累了，帮你引用下代码块，给下面的人看 
#include<stdio.h>
#include<math.h>
#define  E  0.00001
typedef struct {
double x;
double y;
}point;
typedef struct{
point p1;
point p2;
}center;
double distance(point a,point b);
center cluster(point* p,int i);
point cen(point* p,int i,int t);
center part(point* p,int i,point p1,point p2);
int main()
{
point arr[4]={{1,2},{2,4},{2,7},{3,8}};
center c=cluster(arr,4);
printf("%f,%f\n%f,%f",c.p1.x ,c.p1.y,c.p2.x,c.p2.y);
return 0;
}
center cluster(point* p,int i)
{
point a, b,a0,b0;
int n=(i*(i-1))/2;
double max=0;
int x=0;
double d[n]; 
double e1,e2;
center c0;

for(int m=0;m<i-1;m++)
for(int n=1;n<i;n++)
{
d[x]=distance(p[m],p[n]);
if(d[x]>max)
{

max=d[x];
a=p[m];// 把距离最远的点付给a,b 
b=p[n];
}
x++;
}
do
{
a0=a;
b0=b;
c0=part(p,i,a,b);
a=c0.p1 ;
b=c0.p2;
e1=distance(a0,a);
e2=distance(b0,b);
}while((e1>E)|(e2>E));
center c={a,b};
return c;
}
//求两点之间距离函数 
double distance(point a,point b)
{
 double dte;
 dte=sqrt((b.x-a.x)*(b.x-a.x)+(b.y-a.y)*(b.y-a.y));
 return dte;
}
//聚为两类 
center part(point* p,int i,point p1,point p2)
{
point arr1[i];
point arr2[i];
int t1=0;
int t2=0;
for(int n=0;n<i;n++)
{
if(distance(p1,p[n]) < distance(p2,p[n]))
{
arr1[n]=p[n];
t1++;
}

else
{
arr2[n]=p[n];
t2++;
}
}	
p1=cen(arr1,i,t1);
p2=cen(arr2,i,t2);
center c={p1,p2};
return c;
}
//求中心点函数 
point cen(point* p,int i,int t)
{
point pt;
double X=0;


double Y=0;
for(int n=0;n<i;n++)
{
X+=p[n].x;
Y+=p[n].y;
}
pt.x =X/t;
pt.y =Y/t;
return pt;
}
九桔猫 2014-05-14
  • 打赏
  • 举报
 回复
初始化？ 猜测，没见到具体数据。。。。。
lanlvy 2014-05-14
  • 打赏
  • 举报
 回复
引用 5 楼 zhao4zhong1 的回复:
[quote=引用 3 楼 u012514171 的回复:] 怎么感觉还是一样的乱，你手打的吧？
astyle自动整理c,c++,java源代码格式命令行工具源代码。 用astyle --help查看使用帮助。 http://download.csdn.net/detail/zhao4zhong1/5605415[/quote]十分感谢！
C&C++语言学习大全套【精品】
2012-C&C++语言大全.zip C&C++语言教程大礼包.rar C_C+++指针的13份资料.zip C_C++程序员实用大全.zip C++程序设计语言特别版(中文版) .rar C语言大全第4版+源码.rar ...你必须知道的495个C语言问题(完整版).rar
C语言精彩编程100例
C语言精彩编程100例是一本面向C语言初学者的入门教程，旨在通过一系列精心挑选的编程实例，帮助读者...通过系统的练习和深入的理解，你将能够运用C语言解决各种实际问题，为后续的软件开发或系统级编程打下坚实的基础。
C语言常见问题
C语言常见问题
汉诺塔问题详解（C语言）
汉诺塔问题是一个古典的数学问题，本文主要和大家一起用c语言解决汉诺塔问题。
兔子繁衍问题--C语言
兔子繁衍问题——C语言 问题提出 一对兔子，从出生后第三个月起每个月都生一对兔子。小兔子长到第三个月后每个月又生一对兔子。假如兔子都不死，请问第一个月出生的一对兔子，至少需要繁衍到第几个月时兔子总数才能...
C语言

70,036

社区成员

243,247

社区内容

发帖
与我相关
我的任务
社区描述
C语言相关问题讨论
社区管理员
  • C语言
  • 花神庙码农
  • 架构师李肯
加入社区
  • 近7日
  • 近30日
  • 至今

加载中

查看更多榜单
社区公告
暂无公告

试试用AI创作助手写篇文章吧

+ 用AI写文章