C/C++の宿題を片付けます 77代目
■ このスレッドは過去ログ倉庫に格納されています
【質問者へ】
回答者の便宜のため、質問の際は以下を行うことを推奨します。
・質問は【質問テンプレ】を利用してください。
・問題文は、出題されたまま全文を書いてください。
・計算問題は数式をあげ、どのような計算をするのか詳しく説明してください。
・エラーは、その詳細と発生した行を書きましょう。エラーメッセージはコピペしてください。
・後から問題に付け足しするのはやめましょう。付け足しは作業を無駄にしがちです。
・なりすましを防ぐため、トリップを使ってください。名前欄に、「#」に続けて任意の文字列を入力して
投稿すると、その文字列を知らない他人に騙られることを防ぐことができます。
【質問テンプレ】
[1] 授業単元：
[2] 問題文(含コード&リンク)：
[3] 環境
[3.1] OS： (Windows/Linux/等々)
[3.2] コンパイラ名とバージョン： (gcc 3.4 VC 6.0等)
[3.3] 言語： (C/C++/どちらでも可 のいずれか)
[4] 期限： ([yyyy年mm月dd日hh:mmまで] または [無期限] のいずれか)
[5] その他の制限： (どこまで習っているか、標準ライブラリは使ってはいけない等々)
【アップローダー==ラウンジ】(質問が長い時はココ使うと便利 回答者もコードが長ければここに)
http://kansai2channeler.hp.infoseek.co.jp/cgi-bin/joyful/joyful.htm
【C 関数検索 man on WWW】
http://www.linux.or.jp/JM/index.html
【前スレ】
C/C++の宿題を片付けます 76代目
http://pc8.2ch.net/test/read.cgi/tech/1162771475/
【過去スレ一覧】
http://makimo.to/cgi-bin/search/search.cgi?q=C%2B%2B%82%CC%8Fh%91%E8&andor=AND&sf=2&H=&view=table&D=tech&shw=2000
0267デフォルトの名無しさん
2006/11/25(土) 01:20:020268デフォルトの名無しさん
2006/11/25(土) 01:42:19有難うございました！まだまだ勉強不足のようなので、
初めから勉強し直したいと思います！
0270デフォルトの名無しさん
2006/11/25(土) 02:03:46ファイル上に領域を確保してる。
MapViewOfFile (unix 系では mmap) を使うと、ファイルをメモリ上に対応付ける事ができ、
該当するメモリへの読み書きが直接ファイルに対する読み書きになる。
技術的にはスワップファイルと同じ仕組み。
free() する場合と違い、対応付けを解除しても内容はファイル上に保持されるので、アクセス
する必要のある配列を必要のある期間だけ対応付け、必要がなくなったら即対応付けを解除
することで、アドレス空間を節約できる。(なので一度に7配列アクセスしたら当然破綻する)
# ちなみにかつての 16bit Windows は同じような仕組みで少ない実メモリをやりくりしてた。
0271デフォルトの名無しさん
2006/11/25(土) 02:21:15コピペも張られてるみたいだし
【初心者歓迎】C/C++室 Ver.33【環境依存OK】
http://pc8.2ch.net/test/read.cgi/tech/1164376287/
0272デフォルトの名無しさん
2006/11/25(土) 05:12:21[2] 問題文：以下のサンプルを参考に、文字で cos カーブ、sin カーブを描くプログラムの作成
#include <stdio.h>
#include <math.h>
// cc -lm
#define PI 3.1415926
int main( void )
{int deg;
for (deg=0 ; deg <= 360; deg += 20)
printf("%3d: %6.3lf %6.3lf\n",deg ,cos( (double)deg * PI / 180.0 )
,sin( (double)deg * PI / 180.0 ));
}
0 * +
15 * +
30 * +
45 +
60 + *
75 + *
90 + *
105 + *
120 + *
135 + *
150 + *
165 + *
180 + *
195 + *
210 + *
[3.1] OS：Linux
[3.2] コンパイラ名とバージョン：gcc
[3.3] 言語：C
[4] 期限：11月29日
0273デフォルトの名無しさん
2006/11/25(土) 05:29:34/* cc -lm */
#include <stdio.h>
#include <string.h>
#include <math.h>
#define PI 3.1415926
int main(void){
int deg;
for (deg=0 ; deg <= 210; deg += 15){
const int x1=(int)(sin(deg*PI/180.0)*25);
const int x2=(int)(cos(deg*PI/180.0)*25);
int i;
printf("%3d ",deg);
for(i=-25;i<=25;++i)
putchar(i==x1?'*':i==x2?'+':' ');
puts("");
}
return 0;
}
0274デフォルトの名無しさん
2006/11/25(土) 10:48:07[2] 問題文：入力した文字を元に、それを拡大したアスキーアートを描くプログラムの作成
[3.1] OS：Linux
[3.2] コンパイラ名とバージョン：gcc
[3.3] 言語：C
[4] 期限：11月29日
0275デフォルトの名無しさん
2006/11/25(土) 10:54:01AAは自前で用意すんの？
0276デフォルトの名無しさん
2006/11/25(土) 12:25:21[2] 問題文(含コード&リンク)：
簡単なパイプを実現するシェルプログラムを作りなさい。
@ ls | more
. . .
@
のようにコマンドを'|'でつなげると、前のコマンドの
標準出力を後ろのコマンドの標準出力とします。
簡単のため、コマンドは２つだけとします。
パイプの実現にはpipeシステムコールを
使います。
int fds[2];
pipe(fds);
write(fds[1], . . . .);
read(fds[0], . . .);
[3] 環境
[3.1] OS：UNIX
[3.2] コンパイラ名とバージョン：gcc 3.4
[3.3] 言語：C
[4] 期限：2006年11月25日23:00まで
[5] その他の制限： (どこまで習っているか、標準ライブラリは使ってはいけない等々)
ポインタ、構造体、システムコールなどは習いました。
http://kansai2channeler.hp.infoseek.co.jp/cgi-bin/joyful/img/3033.txt
授業でやった↑のサンプルを改造して作るらしいです。
0277デフォルトの名無しさん
2006/11/25(土) 12:53:34まず、サンプルの実行結果だけど、
-----
input the number of vertices n and the number of edges m
4 3
input tail and head of an edge m times
1 2
1 3
2 3
adjacency list
1: 1 2
2: 3
3:
4:
path: a tail head
1: 10 0 0
2: 1 1 2
3: 3 2 3
-----
隣接リスト（adjacency list）で１に対して１、２って正しいのか？
２、３が正しいと思うんだが、
それと最期のpath: と a って何？ a[1]=10になってるしサンプルバグってね？
0278277
2006/11/25(土) 13:04:250280デフォルトの名無しさん
2006/11/25(土) 13:28:340281241
2006/11/25(土) 14:21:31レスありがとうございました。
#include <stdio.h>
#include <string.h>
#define COLSIZE 80
int main(int argc, char *argv[])
{
char buffer[COLSIZE + 1], new[COLSIZE + 1];
FILE *fp;
fgets(new, COLSIZE + 1, stdin);
a[strlen(new)] = '\0';
fp = fopen("new", "r");
fgets(buffer, COLSIZE + 1, fp);
while(fgets(buffer, COLSIZE + 1, fp) != NULL){
printf("%s", buffer);
}
fclose(fp);
return 0;
}
考えた結果こうなったのですが違ってますよね？
どこが違うのでしょうか？
0282デフォルトの名無しさん
2006/11/25(土) 14:25:45ところどころに出てくる
> COLSIZE + 1
には何か特別な意味があるのか？
0283デフォルトの名無しさん
2006/11/25(土) 14:32:131.
＞fp = fopen("new", "r");
"new" では、「newという内容の文字列」になってしまう。
指定したいのは new という「変数」。
2.
一番目の文字を指す時 new[ ? ] の ? に入る数字は？また長さ一の文字列の最後の文字はどこにある?
二番目の文字を指す時 new[ ? ] の ? に入る数字は？また長さ二の文字列の最後の文字はどこにある?
・
n 番目の文字を指す時 new[ ? ] の ? に入る数字は？また長さ n の文字列の最後の文字はどこにある?
0284282
2006/11/25(土) 14:42:560285デフォルトの名無しさん
2006/11/25(土) 14:44:25>a[strlen(new)] = '\0';
new[strlen(new)] = '\0'; もしくは
new[strlen(new)-1] = '\0';
fopen() 直後のfgetは余計
0286デフォルトの名無しさん
2006/11/25(土) 14:56:08あんた馬鹿?
0287デフォルトの名無しさん
2006/11/25(土) 15:13:23[2] 問題文：2つの最大40桁の自然数を読み込み、その和を返すプログラムを作れ。
[3] 環境
[3.1] OS： Linux
[3.2] コンパイラ名： gcc
[3.3] 言語： C
[4] 期限：来週の木曜日（11月30日）まで
[5] その他の制限：構造体、ポインタを軽くならったところです。
入力は文字列で行い、それを型変換し計算するという流れになるかと思いますが、
文字列の扱いに困っております。
0288デフォルトの名無しさん
2006/11/25(土) 15:37:05以前 fj に似たようなのが流れてたなぁと思って検索したら出てきた。
ttp://katsu.watanabe.name/ancientfj/article.php?mid=%3C784lk3%245t%241%40news1.dti.ne.jp%3E
0289デフォルトの名無しさん
2006/11/25(土) 15:42:25適当に書いてみた。
#include<stdio.h>
#include<string.h>
void std_reverse(char*begin,char*end){
char c;while(begin != end && begin != --end)c=*begin,*begin++=*end,*end=c;
}
char*my_strrev(char*p){std_reverse(p,p+strlen(p));return p;}
int main(void){
char a[41]={0},b[41],c,*p=a,*q=b;
scanf("%s%s",a,b);
my_strrev(a);my_strrev(b);
for(c=0;*p||*q||c;c=*p > '9',*p++ -= c * 10)
*p = (*p?*p:'0') + (*q?*q++-'0':0) + c;
puts(my_strrev(a));
return 0;
}
0290デフォルトの名無しさん
2006/11/25(土) 15:46:240291デフォルトの名無しさん
2006/11/25(土) 16:10:21少なくともこれは全く無意味だな。
0292デフォルトの名無しさん
2006/11/25(土) 16:12:140293デフォルトの名無しさん
2006/11/25(土) 16:18:22もしかして、真偽値を確かめてから偽なら後を実行するというアレか！？
0294デフォルトの名無しさん
2006/11/25(土) 16:18:590295デフォルトの名無しさん
2006/11/25(土) 16:39:42その文脈でもしくはと書く以上、
：私にはどちらを依頼人が望んでいるかは分からないが、
：依頼人の望む動作は A と B のどちらかであるように私には思える。
と解釈するのが常識的であって、つまり本人は A にも B にも意味があると
思っていると判断されて当然では。
逆に、私はわざとやったとは考えてないが、もし本人が、一方は全くのハズレで
意味がないと重々承知した上でわざとこのように書いたのなら、その軽重は
別として、相手を混乱させてやろうといった意思があったことになるのでは。
0296デフォルトの名無しさん
2006/11/25(土) 16:49:59回答 彼は頼もしくはない
0297デフォルトの名無しさん
2006/11/25(土) 18:10:21オナニーはひとりで汁
0298デフォルトの名無しさん
2006/11/25(土) 18:14:00まぁ、みせっこもたまにはいいよね
0299デフォルトの名無しさん
2006/11/25(土) 18:23:13[2] 問題文(含コード&リンク)： ttp://kansai2channeler.hp.infoseek.co.jp/cgi-bin/joyful/img/3036.txt
[3] 環境
[3.1] OS： Linux
[3.2] コンパイラ名とバージョン： gcc
[3.3] 言語： C言語
[4] 期限： ２００６年１１月２６日２０時００分まで
[5] その他の制限： 問題で与えられたソースを元にせよという指示がありました。
よくわからない部分は全力で調べるのでよろしくお願いします
0300デフォルトの名無しさん
2006/11/25(土) 18:27:52[2] 問題文：
文字列と数字を与え、文字列の長さが数字より大きいかどうかを判定し、
その後に文字列の長さがどれだけ数字より大きいかを求めて表示するようなプログラムを作成せよ。
文字列の長さの方が小さい場合は何もせず終了せよ。
[3] 環境
[3.1] OS： Linux
[3.2] コンパイラ名： gcc
[3.3] 言語： C
[4] 期限： 11/29まで
[5] その他の制限： まだ始めたばっかなんで全然わかりません。
0301デフォルトの名無しさん
2006/11/25(土) 18:32:14スレ住人おすすめの本かサイトを教えてくれ。
0302デフォルトの名無しさん
2006/11/25(土) 18:33:48brainf*ckがいいんじゃないか？
ぐぐってみろ
0303デフォルトの名無しさん
2006/11/25(土) 18:42:40受け取った引数xの全桁の合計値を返す関数
int sum_dgt(int x) {/* */} を作り実行しなさい。
例えばx=234のときは9となる。（正の整数とする）
[3] 環境
[3.1] OS： Linux
[3.2] コンパイラ名： gcc
[3.3] 言語： C
[4] 期限： 11/28まで
おねがいします
0304281
2006/11/25(土) 18:43:33はい、そういうことです。
>>283氏と>>285氏のレスで解りました。
ありがとうございます。
ただ、>>285氏の
＞fopen() 直後のfgetは余計
というので、なぜfopen() 直後のfgetが必要ないのかは、
いまいち解りませんが。
0306デフォルトの名無しさん
2006/11/25(土) 19:21:44死んでくれ
0307151
2006/11/25(土) 19:23:170308デフォルトの名無しさん
2006/11/25(土) 19:24:06わかっているさ
0309デフォルトの名無しさん
2006/11/25(土) 19:24:11俺より下だったら答えてあげる
0310151
2006/11/25(土) 19:29:110311デフォルトの名無しさん
2006/11/25(土) 19:30:13死ねよKO
0312151
2006/11/25(土) 19:31:58慶應ではないです
特定されると嫌なので書きたくないですが、慶應よりは下です
0313デフォルトの名無しさん
2006/11/25(土) 19:32:29ボンボンというだけで慶應に入って
おまけに卒業に必要な単位は他人任せか
腐れ外道が
0314151
2006/11/25(土) 19:36:09幼稚舎にいた頃、アルバイトで来ていた日大の学生さんを見て、日大はお猿さん養成所も同然だと母が言ってましたが。
0315デフォルトの名無しさん
2006/11/25(土) 19:37:44こっそり本人が書いたんじゃないかと思えてきたよｗ
0316デフォルトの名無しさん
2006/11/25(土) 19:40:380317デフォルトの名無しさん
2006/11/25(土) 19:42:29信じてるよ
0318デフォルトの名無しさん
2006/11/25(土) 19:43:37俺も信じてるよ
0319デフォルトの名無しさん
2006/11/25(土) 19:44:17[1] 授業単元：情報処理工学１
[2] 問題文(含コード&リンク)： 下記
[3] 環境
[3.1] OS：XP
[3.2] コンパイラ名とバージョン：VC++
[3.3] 言語： C++
[4] 期限： 明日まで
[5] その他の制限：
class Vector{
public:
void f(){cout << "f()" << endl;}
};
int main(){
vector<Vector*> vec;
Vector *v;
for(int i = 0; i < 10; i++) { v = new Vector; vec.push_back(v); }
for(vector<Vector*>::Iterator i = vec.begin(); i < vec.end(); i++)
i.***->f();
}
というプログラムを作ってるんですが、動的にクラスを確保して、それを次々にvectorに格納して、最後に表示をするっていうふうにしたんです。
ですが、***の部分をどうやったらいいのか、わかりません。クラスではなく、整数であれば、i をそのまま使えばいいんですけど、クラスだとうまくいきません。
どうしたらいいでしょうか？？
0320デフォルトの名無しさん
2006/11/25(土) 19:44:26前スレのは知らんがこのスレのを消化していこうぜ
下らん流れはさっさと断ち切るに限る
0321151 ◆nEYxVzEafI
2006/11/25(土) 19:46:590323デフォルトの名無しさん
2006/11/25(土) 19:50:57//受け取った引数xの全桁の合計値を返す関数
int sum_dgt(int x)
{
int digit = 1;
int sum = 0;
int i ;
int xx ;
int xxx = x;
for(i = 1; x > 1 ;i++){
x /=10;
digit *= 10;
}
for(;digit > 0; digit/=10){
xx = xxx;
xx /= digit;
sum += xx;
xxx-=xx*digit;
}
return sum;
}
0325デフォルトの名無しさん
2006/11/25(土) 19:53:32>>151 （今日中）
>>274 （29日まで）
>>276 （今日の23時まで）
>>299 （明日の20時まで）
>>300 （29日まで）
>>319 （明日まで）
期限切れを消しといた
0326デフォルトの名無しさん
2006/11/25(土) 19:54:32もっとゆとりを持って依頼に来てくれよｗ
0328デフォルトの名無しさん
2006/11/25(土) 20:02:12│>>90 （今日中）
│>>299 （明日の20時まで）
│>>319 （明日まで）
↓>>151 >>274 >>300 （29日まで）
まあまったり待てや
0329デフォルトの名無しさん
2006/11/25(土) 20:03:380330デフォルトの名無しさん
2006/11/25(土) 20:06:56#include <stdio.h>
#include <string.h>
//'\0'は数えないものとする
void Compar(int x, char s[])
{
unsigned int len;
len = strlen(s);
if(x <= len)
printf("文字の長さと数値の差は%d \n",len - x);
else
return;
}
0331デフォルトの名無しさん
2006/11/25(土) 20:09:00│>>299 （明日の20時まで）
│>>319 （明日まで）
↓>>151 >>274 （29日まで）
0332323,300
2006/11/25(土) 20:10:570333名無しさん
2006/11/25(土) 20:13:471を入力されるとIMEがONになって、'あ' と表示するプログラムの作成中です。
以下のコードではIMEがONになってくれなくて、'a' と表示されます。
ご指摘お願いします。
void fnConvert(MSG *pMsg)
{
static int iCnt= 1;
if( (pMsg->wParam == 0x31) && (pMsg->message == WM_KEYDOWN) )
{
pMsg->wParam= 0x41;
break;
}
}
0334デフォルトの名無しさん
2006/11/25(土) 20:14:23俺は274でもやるか
本人降臨してほしいが
0335デフォルトの名無しさん
2006/11/25(土) 20:17:010336デフォルトの名無しさん
2006/11/25(土) 20:18:32しかし宿題出された瞬間にここに出されても萎えるという
0337デフォルトの名無しさん
2006/11/25(土) 20:33:21http://kansai2channeler.hp.infoseek.co.jp/cgi-bin/joyful/img/3037.c
http://kansai2channeler.hp.infoseek.co.jp/cgi-bin/joyful/img/3038.c
0338デフォルトの名無しさん
2006/11/25(土) 21:03:07[2] 問題文(含コード&リンク)： あるテキストファイルに含まれるA～Zの個数をカウントするプログラムをつくれ。
[3] 環境
[3.1] OS： Windows
[3.2] コンパイラ名とバージョン： lcc
[3.3] 言語： C
[4] 期限： 11月28日まで
お願いします。どう検索すればいいかも分からないので教えていただきたいです。
答えがほしいですが、ヒントをいただければ自分でも努力はします。
0339デフォルトの名無しさん
2006/11/25(土) 21:05:42fgetc()を使ってEOFになるまで一文字ずつ読み込み、isalpha() && isupper()が
真になる時に配列count[26]のうち対応する文字のカウントを＋１する。
0340デフォルトの名無しさん
2006/11/25(土) 21:15:57質問されても面倒なので答えを書いておく
#include <stdio.h>
#include <stdlib.h>
#include <ctype.h>
int main(void)
{
int chr[26] = {0}, c, i;
FILE *fp;
if ((fp =fopen("count1.c", "r")) == NULL) /* ファイル名は適宜変更 */
exit(1);
while ((c = fgetc(fp)) != EOF)
if (isalpha(c) && isupper(c))
chr[c - 'A']++; /* Windows上で動かす事を前提とする */
for (i = 0; i < 26; i++)
printf("%c : %dcount(s)\n", i + 'A', chr[i]);
return 0;
}
0341デフォルトの名無しさん
2006/11/25(土) 21:16:380342デフォルトの名無しさん
2006/11/25(土) 21:23:57[2] 問題文(含コード&リンク)： テキスト(英文)を読み込んで単語を全て取り出し、重複する単語は1つにまとめ、その単語がいくつあったかabc順に並べて表示するプログラムを書きなさい。(.,"などは含まない)
[3] 環境
[3.1] OS：WindowsXP
[3.2] lcc
[3.3] 言語：C
[4] 期限：2006年11月26日24:00まで
プログラミング初心者です。申し訳ありませんがどなたかお願いします。
0343デフォルトの名無しさん
2006/11/25(土) 21:26:08#include <stdio.h>
#include <ctype.h>
#define NUMBER 26
main(int argc, char **argv){FILE *fp; int i, c, chr[NUMBER]={0};
if((fp=fopen(argv[1],"r")) == NULL) puts("need 1 parameter");
else {
while ((c=fgetc(fp))!=EOF) if (isalpha(c)) chr[toupper(c)-'A']++;
fclose(fp);}
for (i=0; i<NUMBER; i++) printf("%c: %3d time%s\n",i+'A',chr[i], chr[i]>1?"s":"");
getchar();return 0;}
0344デフォルトの名無しさん
2006/11/25(土) 21:26:55ていうか、K&Rもんばっかだな
0345287
2006/11/25(土) 21:30:50ありがとうございます。が、なんか文字化けしちゃって私の環境では見れませんね・・・
わざわざ張ってくださったのに申し訳ないです
>>289
これはこれは、わざわざ作ってくださってありがとうございます
しかし私のポインタの知識が未熟なせいか、今現在はまだ読み解けない部分がありますね
けれども、これはとても参考になります。ほんとにありがとうございました
0346デフォルトの名無しさん
2006/11/25(土) 21:32:33base64知らないのか
0347デフォルトの名無しさん
2006/11/25(土) 21:41:190348デフォルトの名無しさん
2006/11/25(土) 21:50:500349デフォルトの名無しさん
2006/11/25(土) 21:54:34ホント悪かった。
0350デフォルトの名無しさん
2006/11/25(土) 21:55:31│>>319 >>342 （明日まで）
↓>>151 >>274 （29日まで）
未解消
0351287
2006/11/25(土) 21:58:09すみません知りませんでした
一応ググって復元してみたら無事プログラムソースを見ることができました
お手数かけまして申し訳ありません
0352デフォルトの名無しさん
2006/11/25(土) 22:00:19[2] 問題文(含コード&リンク)：
y=sin2x*cosx(0≦x≦π/4)をgnuplotで図示せよ。
gnuplotの関数表示機能は使わず、C言語でプログラムを組んで、その出力
データをgnuplotで読み込んで図示させること。刻み数は1000とする。
[3] 環境
[3.1] OS： Linux
[3.2] コンパイラ名とバージョン：gcc
[3.3] 言語： C
[4] 期限： 来週の火曜日の20時までです。
[5] その他の制限：特になし。
出力データをファイルに書き込み、gnuplotでそのファイルを読み込み
グラフを表示させる、といことです。
よろしくお願いします。
0353デフォルトの名無しさん
2006/11/25(土) 22:05:28俺遺伝的アルゴリズムやったことないんだ
計算結果はdoubleで返すことを期待されてると思うんだが
double型にビット演算ってできないよな？
何か良いアイデア知ってる火といないかな
0354デフォルトの名無しさん
2006/11/25(土) 22:07:25class Vector {
public:
void f(){ std::cout << "f()" << std::endl; }
};
int main()
{
std::vector<Vector*> vec;
Vector *v;
for (int i = 0; i < 10; i++) {
v = new Vector;
vec.push_back(v);
}
for (std::vector<Vector*>::iterator i = vec.begin(); i < vec.end(); ++i)
(*i)->f();
}
0355デフォルトの名無しさん
2006/11/26(日) 00:15:05gnuplotに読み込ますデータの仕様くらい用意してから書き込めよ。
0356デフォルトの名無しさん
2006/11/26(日) 00:16:38つ union
0357デフォルトの名無しさん
2006/11/26(日) 00:40:25お こんな要素があったのか バッチリだな
ありがとう
0358352
2006/11/26(日) 00:55:37すいません、
＞gnuplotに読み込ますデータの仕様くらい用意してから書き込めよ。
というのは、出力データをどんな感じでファイルに書き込むか、という
ことですか？
xiとyiとしたら
0.0000 0.14356
0.0100 0.16344
. .
. .
. .
という感じです。
0359デフォルトの名無しさん
2006/11/26(日) 01:05:49#include<stdio.h>
#include<math.h>
#define N 1000
double f(double x){
return sin(2*x)*cos(x);
}
int main(void){
double x, dx=(M_PI/4)/N;
int i;
for(i=0;i<=N;i++){
x=i*dx;
printf("%f %f\n", x, f(x));
}
return 0;
}
0360デフォルトの名無しさん
2006/11/26(日) 02:23:49excuse me for writing in English, i haven't installed japanese IME
so i can't type in japanese
i'm working on my project in my computer science class and
found a very interesting result. that is recursive merge-sort
function is faster than non-recursive version.
i posted the same question on cprogramming.com so that
the details and my source code is provided at
http://cboard.cprogramming.com/showthread.php?p=612918#post612918
0361353
2006/11/26(日) 04:03:12さすがに人に見せられるレベルじゃないのでもう少し手を入れるが
もし締め切りが近かったらアレなので89はまだいたら締め切り示せ
0363その１
2006/11/26(日) 09:13:28#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#include <float.h>
#include <math.h>
#define N 500//世代内の個体数
#define G 1000//最大世代数
#define Kousa 0.5//交差率。0～1
#define Totuzen 0.1//突然変異率。0～1だが小さめに設定すること
#define Pi 3.14159265358979
union uGENE
{
int i[2];
double d;
};
double hyouka (double x);
void kousa(int *a,int *b);
void mutation(int *a,int *b);
int mycomp(const void *_a, const void *_b);
int choice(int i);
0364その２
2006/11/26(日) 09:14:14int g,i;
union uGENE future[N+1];
double r,present[N];
srand((unsigned) time(NULL));
for(i=0;i<N;i++) present[i]=3.0*((double)rand()/RAND_MAX)-1.0;
for(g=0;g<G;g++){
printf("%d世代目計算中・・・\n",g);
qsort(present, N, sizeof(double), mycomp);
for(i=0;i<N;i++) {
r = (double)rand()/RAND_MAX;
future[i].d=present[choice(i)];
if(r<Kousa){
//交叉
future[i+1].d=present[choice(i+1)];
kousa(future[i].i,future[i+1].i);
kousa(future[i].i+1,future[i+1].i+1);
i++;
}
else if(r<Kousa+Totuzen){
//突然変異
mutation(future[i].i,future[i].i+1);
}
//何もしなければコピーになる
}
//世代交代
for(i=0;i<N;i++) present[i]=future[i].d;
}
qsort(present, N, sizeof(double), mycomp);
printf("x=%fの時f(x)は%f",present[0],hyouka(present[0]));
}
0365その３
2006/11/26(日) 09:15:20//評価関数はx{sin(10πx)}+2cos(10πx)+1.0をそのまま用いるが
//0以下の場合はマシンエプシロンの10倍を返す。
double hyouka (double x){
if(x<-1||x>2) return DBL_EPSILON;
double y,tmp;
tmp=10.0*Pi*x;
y=x*sin(tmp)+2*cos(tmp)+1.0;
if (y<=0) return DBL_EPSILON*10;
else return y;
}
//交叉関数。一様交叉を行う
void kousa(int *a,int *b){
//交叉するビットを示すマスクを作る
unsigned int mask=0;
//マスクの各ビットにランダムに0/1を入れる
for(int i=0;i<32;i++){
mask+=rand()%2;
mask<<=1;
}
//aとbで内容が同じビットは交叉する必要がない
mask=mask&(*a^*b);
*a^=mask;
*b^=mask;
}
0366その４
2006/11/26(日) 09:16:00void mutation(int *a,int *b){
unsigned int mask=1;
int r=rand()%64;
if(r>=32){
r-=32;
mask<<=r;
*a^=mask;
}
else{
mask<<=r;
*b^=mask;
}
}
int mycomp(const void *_a, const void *_b)
{
double a = *(double *)_a;
double b = *(double *)_b;
double tmp=hyouka(b)-hyouka(a);
if (tmp<0) {
return -1;
} else if (tmp > 0) {
return 1;
} else {
return 0;
}
}
■ このスレッドは過去ログ倉庫に格納されています