herrDeng網內搜尋
自訂搜尋
熱門文章
-
計算你的BMI
-
1. 利用遞迴input n算2^n 2. 用C算GCD(3333,456)
-
請用 C/C++/java寫一簡易程式
-
33, 45, 87, 99, 27 算平均
-
一、 試利用switch case由程式中直接輸入一個1~7 之間的整數day，代表星期一到星期日。若day 的值是1，則印出 "星期一"，若day 的值是2，則印出 "星期二"，若day 的值是7，則印出 "星期日"，...
-
BMI (Body Mass Index)= 體重 (kg) / 身高 ^2(m 2 ) 在台灣，行政院衛生署乃根據其相關研究，於2002年4月公佈台灣成人肥胖標準： BMI＜18.5 為過輕， 18.5≦BMI＜24 為正常體重， 24≦BMI＜27 為過重， BM...
-
需要練習的Java程式：
-
用 C/C++算上學期加權平均
36 則留言:
#include <iostream>
#include <cstdlib>
#include <random>
#include <ctime>
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution<double> X(70, 10);
double a[1000];
double sum = 0;
double sum2 = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9)cout << endl;
}
sum = sum / 1000.0;
cout << "平均數:" << sum << endl;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i]) - (sum * sum);
}
sum2 = sum2 / 1000;
double s = sqrt(sum2);
cout << "標準差:" << s << endl;
cout << (clock()-seed)/1000.0 <<"sec." << endl;
system("Pause");
return 0;
}
#include <iostream>
#include <cstdlib>
#include <random>
#include <ctime>
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution<double> X(70, 10);
double a[1000];
double sum = 0;
double sum2 = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9)cout << endl;
}
sum = sum / 1000.0;
cout << "平均數:" << sum << endl;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i]) - (sum * sum);
}
sum2 = sum2 / 1000;
double s = sqrt(sum2);
cout << "標準差:" << s << endl;
cout << (clock()-seed)/1000.0 <<"sec." << endl;
system("Pause");
return 0;
}
#include <iostream>
#include <cstdlib>
#include <random>
#include <ctime>
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution<double> X(70, 10);
double a[1000];
double sum = 0;
double sum2 = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9)cout << endl;
}
sum = sum / 1000.0;
cout << "平均數:" << sum << endl;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i]) - (sum * sum);
}
sum2 = sum2 / 1000;
double s = sqrt(sum2);
cout << "標準差:" << s << endl;
cout << (clock()-seed)/1000.0 <<"sec." << endl;
system("Pause");
return 0;
}
#include <iostream>
#include<cstdlib>
#include<random>
#include<ctime>
#include<cmath>
using namespace std;
int main(){
int seed = clock();
mt19937_64 gen(seed);
normal_distribution<double> x(70, 10);
int a[1000];
double sum = 0;
double s = 0;
double sum2;
for (int i = 0; i < 1000; i++){
a[i] = (int)(round(x(gen)));
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9){
cout << endl;
}
}
sum2 = sum / 1000.0;
for (int i = 0; i < 1000; i++){
// s += (a[i] - sum2)*(a[i] - sum2);
s += (a[i] * a[i]) - (sum2*sum2);
}
s /= 1000;
s = sqrt(s);
cout << sum2;
cout << endl;
cout << s;
cout << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return (0);
}
#include
#include
#include
#include
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution X(70, 10);
double a[1000];
double sum = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9) cout << endl;
}
double av = sum / 1000.0;
cout << av << "\n";
double sum2 = 0;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i] - av*av);
}
double var = sum2 / 1000.0;
double s = sqrt(var);
cout << s << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return 0;
}
#include<iostream>
#include<cstdlib>
#include<random>
#include<ctime>
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution<double> x(70, 10);
int a[1000];
double s = 0;
for (int i = 0; i < 1000; i++){
a[i] = x(gen);
cout << a[i] << "\t";
s = a[i]+s;
if (i % 10 == 9)cout << endl;
}
cout << (clock() - seed) / 10000.0 << "sec.\n";
double av = s / 1000.0;
double sum2 = 0;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] - av)*(a[i] - av);
}
double var = sum2 / 1000.0;
double c = sqrt(var);
cout << c << "標準差";
cout << endl;
system("Pause");
return 0;
}
#include <iostream>
#include <cstdlib>
#include <random>
#include <ctime>
#include <cmath>
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution<double> dis(70, 10);
double a[1000];
double sum = 0;
double sum2 = 0;
for (int i = 0; i < 1000; i++){
a[i] = dis(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9) cout << endl;
}
double av = sum/1000.0;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i]- av*av);
}
double var = sum2 / 1000.0;
cout << sum / 1000.0 << endl;
double s = sqrt(var);
cout << s << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
cout << endl;
system("Pause");
return 0;
}
#include
#include
#include
#include
#include
using namespace std;
int main(){
int seed = clock();
mt19937_64 gen(seed);
normal_distribution x(70, 10);
int a[1000];
double sum = 0;
double s = 0;
double sum2;
for (int i = 0; i < 1000; i++){
a[i] = (int)(round(x(gen)));
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9){
cout << endl;
}
}
sum2 = sum / 1000.0;
for (int i = 0; i < 1000; i++){
s += (a[i] * a[i]) - (sum2*sum2);
}
s /= 1000;
s = sqrt(s);
cout << sum2;
cout << endl;
cout << s;
cout << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return (0);
}
#include
#include
#include
#include
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution X(70, 10);
double a[1000];
double sum = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9) cout << endl;
}
double av = sum / 1000.0;
cout << av << "\n";
double sum2 = 0;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i] - av*av);
}
double var = sum2 / 1000.0;
double s = sqrt(var);
cout << s << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return 0;
}
#include
#include
#include
#include
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution X(70, 10);
double a[1000];
double sum = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9) cout << endl;
}
double av = sum / 1000.0;
cout << av << "\n";
double sum2 = 0;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i] - av*av);
}
double var = sum2 / 1000.0;
double s = sqrt(var);
cout << s << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return 0;
}
#include
#include
#include
#include
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution X(70, 10);
double a[1000];
double sum = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9) cout << endl;
}
double av = sum / 1000.0;
cout << av << "\n";
double sum2 = 0;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i] - av*av);
}
double var = sum2 / 1000.0;
double s = sqrt(var);
cout << s << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return 0;
}
#include <iostream>
#include <cstdlib>
#include <random>
#include <ctime>
#include <cmath>
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution<double>X(70, 10);
double a[1000];
double sum = 0;
for (int i = 0; i < 100; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9)cout << endl;
}
double av = sum / 1000.0;
cout << av << "\n";
double sum2 = 0;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i] - av*av);
}
double var = sum2 / 1000.0;
double s = sqrt(var);
cout << s << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return 0;
}
#include <iostream>
#include<cstdlib>
#include<random>
#include<ctime>
#include<cmath>
using namespace std;
int main(){
int seed = clock();
mt19937_64 gen(seed);
normal_distribution<double> x(70, 10);
int a[1000];
double sum = 0;
double s = 0;
double sum2;
for (int i = 0; i < 1000; i++){
a[i] = (int)(round(x(gen)));
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9){
cout << endl;
}
}
sum2 = sum / 1000.0;
for (int i = 0; i < 1000; i++){
s += (a[i] * a[i]) - (sum2*sum2);
}
s /= 1000;
s = sqrt(s);
cout << sum2;
cout << endl;
cout << s;
cout << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return (0);
}
#include
#include
#include
#include
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution X(70, 10);
double a[1000];
double sum = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9) cout << endl;
}
double av = sum / 1000.0;
cout << av << "\n";
double sum2 = 0;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i] - av*av);
}
double var = sum2 / 1000.0;
double s = sqrt(var);
cout << s << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return 0;
}
#include <iostream> //程式如果用到cout或cin就必須引用iostream
#include <cstdlib>
#include <random>
#include <ctime>
#include <cmath>
using namespace std;
int main(void) {
int seed = clock();
mt19937_64 gen(seed);
normal_distribution<double> dis(70, 10);
int a[1000];
double sum = 0;
for (int i = 0; i < 1000; i++){
a[i] = (int)(round(dis(gen)));
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9)cout << endl;
}
double av = sum / 1000.0;
cout << av << "\n";
double sum2 = 0;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i] - av*av);
}
double var = sum2 / 1000.0;
double s = sqrt(var);
cout << s << "\n";
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return 0;
}
#include <iostream>
#include <cstdlib>
#include <random>
#include <ctime>
#include <cmath>
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution<double> X(70, 10);
double a[1000];
double sum = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9)cout << endl;
}
double av = sum / 1000.0;
cout << av << "\n";
double sum2 = 0;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i] - av*av);
}
double var = sum2 / 1000.0;
double s = sqrt(var);
cout << s << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return 0;
}
#include <iostream>
#include <cstdlib>
#include <random>
#include <ctime>
#include <cmath>
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution<double>X(70, 10);
double a[1000];
double sum = 0;
for (int i = 0; i < 100; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9)cout << endl;
}
double av = sum / 1000.0;
cout << av << "\n";
double sum2 = 0;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i] - av*av);
}
double var = sum2 / 1000.0;
double s = sqrt(var);
cout << s << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return 0;
}
#include <iostream>
#include <cstdlib>
#include <random>
#include <ctime>
#include <cmath>
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution<double> X(70, 10);
double a[1000];
double sum = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9)cout << endl;
}
double av = sum / 1000.0;
cout << av << "\n";
double sum2 = 0;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i] - av*av);
}
double var = sum2 / 1000.0;
double s = sqrt(var);
cout << s << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return 0;
}
#include
#include
#include
#include
#include
using namespace std;
int main(){
int seed = clock();
mt19937_64 gen(seed);
normal_distributionX(70, 10);
double a[1000];
double sum = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9)cout << endl;
}
double av = sum / 1000.0;
cout << av << "\n";
double sum2 = 0;
for (int i = 0; i < 1000; i++){
sum2+=(a[i] * a[i] - av*av);
}
double var = sum2 / 1000.0;
double s = sqrt(var);
cout << s << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return 0;
}
#include <iostream>
#include <cstdlib>
#include <random>
#include <ctime>
#include <cmath>
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution<double> dis(70, 10);
double a[100];
double sum = 0;
for (int i = 0; i < 100; i++){
a[i] = (int)(round(dis(gen)));
cout << a[i] << "\t";
sum+= a[i];
if (i % 10 == 9)cout << endl;
}
double av = sum / 100.0;
cout <<"average= "<< av << "\n";
double sum2 = 0;
for (int i = 0; i < 100; i++){
sum2 += (a[i] * a[i] - av*av);
}
double var = sum2 / 100.0;
double s = sqrt(var);
cout <<"standard deviation= "<< s << endl;
cout << (clock()-seed/100.0)<<"sec. \n";
system("pause");
return 0;
}
#include
#include
#include
#include
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution X(70, 10);
double a[1000];
double sum = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9) cout << endl;
}
double av = sum / 1000.0;
cout << av << "\n";
double sum2 = 0;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i] - av*av);
}
double var = sum2 / 1000.0;
double s = sqrt(var);
cout << s << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return 0;
}
B10233139 陳鈺璇
#include
#include
#include
#include
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution X(70, 10);
double a[1000];
double sum = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9) cout << endl;
}
double av = sum / 1000.0;
cout << av << "\n";
double sum2 = 0;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i] - av*av);
}
double var = sum2 / 1000.0;
double s = sqrt(var);
cout << s << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return 0;
}
#include
#include
#include
#include
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution X(70, 10);
double a[1000];
double sum = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9) cout << endl;
}
double av = sum / 1000.0;
cout << av << "\n";
double sum2 = 0;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i] - av*av);
}
double var = sum2 / 1000.0;
double s = sqrt(var);
cout << s << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return 0;
}
#include
#include
#include
#include
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution X(70, 10);
double a[1000];
double sum = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9) cout << endl;
}
double av = sum / 1000.0;
cout << av << "\n";
double sum2 = 0;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i] - av*av);
}
double var = sum2 / 1000.0;
double s = sqrt(var);
cout << s << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return 0;
}
b10233137
#include
#include
#include
#include
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution X(70, 10);
double a[1000];
double sum = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9) cout << endl;
}
double av = sum / 1000.0;
cout << av << "\n";
double sum2 = 0;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i] - av*av);
}
double var = sum2 / 1000.0;
double s = sqrt(var);
cout << s << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return 0;
}
B10233169
#include
#include
#include
#include
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution X(70, 10);
double a[1000];
double sum = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9) cout << endl;
}
double av = sum / 1000.0;
cout << av << "\n";
double sum2 = 0;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i] - av*av);
}
double var = sum2 / 1000.0;
double s = sqrt(var);
cout << s << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return 0;
}
B10233139
#include
#include
#include
#include
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution X(70, 10);
double a[1000];
double sum = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9) cout << endl;
}
double av = sum / 1000.0;
cout << av << "\n";
double sum2 = 0;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i] - av*av);
}
double var = sum2 / 1000.0;
double s = sqrt(var);
cout << s << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return 0;
}
B10233169
#include
#include
#include
#include
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution X(70, 10);
double a[1000];
double sum = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9) cout << endl;
}
double av = sum / 1000.0;
cout << av << "\n";
double sum2 = 0;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i] - av*av);
}
double var = sum2 / 1000.0;
double s = sqrt(var);
cout << s << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return 0;
}
#include <iostream>
#include <cstdlib>
#include <random>
#include <ctime>
#include <cmath>
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution<double> X(70, 10);
double a[1000];
double sum = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9)cout << endl;
}
double av = sum / 1000.0;
cout << av << "\n";
double sum2 = 0;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i] - av*av);
}
double var = sum2 / 1000.0;
double s = sqrt(var);
cout << s << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return 0;
}
#include
#include
#include
#include
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution X(70, 10);
double a[1000];
double sum = 0;
double sum2 = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9)cout << endl;
}
sum = sum / 1000.0;
cout << "平均數:" << sum << endl;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i]) - (sum * sum);
}
sum2 = sum2 / 1000;
double s = sqrt(sum2);
cout << "標準差:" << s << endl;
cout << (clock()-seed)/1000.0 <<"sec." << endl;
system("Pause");
return 0;
}
#include
#include
#include
#include
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution X(70, 10);
double a[1000];
double sum = 0;
double sum2 = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9)cout << endl;
}
sum = sum / 1000.0;
cout << "平均數:" << sum << endl;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i]) - (sum * sum);
}
sum2 = sum2 / 1000;
double s = sqrt(sum2);
cout << "標準差:" << s << endl;
cout << (clock()-seed)/1000.0 <<"sec." << endl;
system("Pause");
return 0;
}
#include
#include
#include
#include
#include
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution X(70, 10);
double a[1000];
double sum = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9)cout << endl;
}
double av = sum / 1000.0;
cout << av << "\n";
double sum2 = 0;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i] - av*av);
}
double var = sum2 / 1000.0;
double s = sqrt(var);
cout << s << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return 0;
}
ImageMatrix convert2green(ImageMatrix m){
ImageMatrix mm = ImageMatrix(m.height, m.width);
for (int y = 0; y < m.height; y++)
for (int x = 0; x < m.rowsize; x += 3){
BYTE B = m.term[y*m.rowsize + x + 0];
BYTE G = m.term[y*m.rowsize + x + 1];
BYTE R = m.term[y*m.rowsize + x + 2];
mm.term[y*m.rowsize + x + 0] = mm.term[y*m.rowsize + x + 2] = 0;
mm.term[y*m.rowsize + x + 1] = G;
}
return mm;
}
ImageMatrix convert2red(ImageMatrix m){
ImageMatrix mm = ImageMatrix(m.height, m.width);
for (int y = 0; y < m.height; y++)
for (int x = 0; x < m.rowsize; x += 3){
BYTE B = m.term[y*m.rowsize + x + 0];
BYTE G = m.term[y*m.rowsize + x + 1];
BYTE R = m.term[y*m.rowsize + x + 2];
mm.term[y*m.rowsize + x + 0] = mm.term[y*m.rowsize + x + 1] = 0;
mm.term[y*m.rowsize + x +2] = G;
}
return mm;
}
ImageMatrix convert2blue(ImageMatrix m){
ImageMatrix mm = ImageMatrix(m.height, m.width);
for (int y = 0; y < m.height; y++)
for (int x = 0; x < m.rowsize; x += 3){
BYTE B = m.term[y*m.rowsize + x + 0];
BYTE G = m.term[y*m.rowsize + x + 1];
BYTE R = m.term[y*m.rowsize + x + 2];
mm.term[y*m.rowsize + x +1] = mm.term[y*m.rowsize + x + 2] = 0;
mm.term[y*m.rowsize + x + 0] =G;
}
return mm;
}
ImageMatrix convert2half(ImageMatrix m){
ImageMatrix mm = ImageMatrix(m.height, m.width);
for (int y = 0; y<m.height; y++)
for (int x = 0; x<m.rowsize; x += 3){
BYTE B = m.term[y*m.rowsize + x];
BYTE G = m.term[y*m.rowsize + x + 1];
BYTE R = m.term[y*m.rowsize + x + 2];
BYTE gray = (BYTE)(0.2989*R + 0.5870*G + 0.1140*B + 0.5);
if (y<m.height/2)
mm.term[y*m.rowsize + x] = mm.term[y*m.rowsize + x + 1] =
mm.term[y*m.rowsize + x + 2] = gray;
else{
mm.term[y*m.rowsize + x + 0] = B;
mm.term[y*m.rowsize + x + 1] = G;
mm.term[y*m.rowsize + x + 2] = R;
}
}
return mm;
}
ImageMatrix convert2green(ImageMatrix m){
ImageMatrix mm = ImageMatrix(m.height, m.width);
for (int y = 0; y < m.height; y++)
for (int x = 0; x < m.rowsize; x += 3){
BYTE B = m.term[y*m.rowsize + x + 0];
BYTE G = m.term[y*m.rowsize + x + 1];
BYTE R = m.term[y*m.rowsize + x + 2];
mm.term[y*m.rowsize + x + 0] = mm.term[y*m.rowsize + x + 2] = 0;
mm.term[y*m.rowsize + x + 1] = G;
}
return mm;
}
ImageMatrix convert2red(ImageMatrix m){
ImageMatrix mm = ImageMatrix(m.height, m.width);
for (int y = 0; y < m.height; y++)
for (int x = 0; x < m.rowsize; x += 3){
BYTE B = m.term[y*m.rowsize + x + 0];
BYTE G = m.term[y*m.rowsize + x + 1];
BYTE R = m.term[y*m.rowsize + x + 2];
mm.term[y*m.rowsize + x + 0] = mm.term[y*m.rowsize + x + 1] = 0;
mm.term[y*m.rowsize + x + 2] = G;
}
return mm;
}
ImageMatrix convert2blue(ImageMatrix m){
ImageMatrix mm = ImageMatrix(m.height, m.width);
for (int y = 0; y < m.height; y++)
for (int x = 0; x < m.rowsize; x += 3){
BYTE B = m.term[y*m.rowsize + x + 0];
BYTE G = m.term[y*m.rowsize + x + 1];
BYTE R = m.term[y*m.rowsize + x + 2];
mm.term[y*m.rowsize + x + 1] = mm.term[y*m.rowsize + x + 2] = 0;
mm.term[y*m.rowsize + x + 0] = G;
}
return mm;
}
ImageMatrix convert2half(ImageMatrix m){
ImageMatrix mm = ImageMatrix(m.height, m.width);
for (int y = 0; y<m.height; y++)
for (int x = 0; x<m.rowsize; x += 3){
BYTE B = m.term[y*m.rowsize + x];
BYTE G = m.term[y*m.rowsize + x + 1];
BYTE R = m.term[y*m.rowsize + x + 2];
BYTE gray = (BYTE)(0.2989*R + 0.5870*G + 0.1140*B + 0.5);
if (y<m.height / 2)
mm.term[y*m.rowsize + x] = mm.term[y*m.rowsize + x + 1] =
mm.term[y*m.rowsize + x + 2] = gray;
else{
mm.term[y*m.rowsize + x + 0] = B;
mm.term[y*m.rowsize + x + 1] = G;
mm.term[y*m.rowsize + x + 2] = R;
}
}
return mm;
}
#include
#include
#include
#include
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution X(70, 10);
double a[1000];
double sum = 0;
double sum2 = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9)cout << endl;
}
sum = sum / 1000.0;
cout << "平均數:" << sum << endl;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i]) - (sum * sum);
}
sum2 = sum2 / 1000;
double s = sqrt(sum2);
cout << "標準差:" << s << endl;
cout << (clock()-seed)/1000.0 <<"sec." << endl;
system("Pause");
return 0;
}
ImageMatrix convert2green(ImageMatrix m){
ImageMatrix mm = ImageMatrix(m.height, m.width);
for (int y = 0; y < m.height; y++)
for (int x = 0; x < m.rowsize; x += 3){
BYTE B = m.term[y*m.rowsize + x + 0];
BYTE G = m.term[y*m.rowsize + x + 1];
BYTE R = m.term[y*m.rowsize + x + 2];
mm.term[y*m.rowsize + x + 0] = mm.term[y*m.rowsize + x + 2] = 0;
mm.term[y*m.rowsize + x + 1] = G;
}
return mm;
}
ImageMatrix convert2red(ImageMatrix m){
ImageMatrix mm = ImageMatrix(m.height, m.width);
for (int y = 0; y < m.height; y++)
for (int x = 0; x < m.rowsize; x += 3){
BYTE B = m.term[y*m.rowsize + x + 0];
BYTE G = m.term[y*m.rowsize + x + 1];
BYTE R = m.term[y*m.rowsize + x + 2];
mm.term[y*m.rowsize + x + 0] = mm.term[y*m.rowsize + x + 1] = 0;
mm.term[y*m.rowsize + x +2] = G;
}
return mm;
}
ImageMatrix convert2blue(ImageMatrix m){
ImageMatrix mm = ImageMatrix(m.height, m.width);
for (int y = 0; y < m.height; y++)
for (int x = 0; x < m.rowsize; x += 3){
BYTE B = m.term[y*m.rowsize + x + 0];
BYTE G = m.term[y*m.rowsize + x + 1];
BYTE R = m.term[y*m.rowsize + x + 2];
mm.term[y*m.rowsize + x +1] = mm.term[y*m.rowsize + x + 2] = 0;
mm.term[y*m.rowsize + x + 0] =G;
}
return mm;
}
ImageMatrix convert2half(ImageMatrix m){
ImageMatrix mm = ImageMatrix(m.height, m.width);
for (int y = 0; y<m.height; y++)
for (int x = 0; x<m.rowsize; x += 3){
BYTE B = m.term[y*m.rowsize + x];
BYTE G = m.term[y*m.rowsize + x + 1];
BYTE R = m.term[y*m.rowsize + x + 2];
BYTE gray = (BYTE)(0.2989*R + 0.5870*G + 0.1140*B + 0.5);
if (y<m.height/2)
mm.term[y*m.rowsize + x] = mm.term[y*m.rowsize + x + 1] =
mm.term[y*m.rowsize + x + 2] = gray;
else{
mm.term[y*m.rowsize + x + 0] = B;
mm.term[y*m.rowsize + x + 1] = G;
mm.term[y*m.rowsize + x + 2] = R;
}
}
return mm;
}
張貼留言