Aug 16
不长不短,150左右的一个类。。
下面是源码:
Aug 16

C++矩阵类(class Matrix) 不指定

felix021 @ 2008-8-16 00:41 [IT » 程序设计] 评论(3) , 引用(0) , 阅读(14116) | Via 本站原创
@ 2009-05-25 其实早就发现了这个代码有严重的内存泄露问题,不过懒得改了。。。

------

实现重载了等于=  加法+  减法-  矩阵乘法*  矩阵幂^(使用快速幂算法)  输出<<  输入>>,用起来会比较方便了,不过就是有点长。

查了好久快速幂算法,网上居然没有现成的,于是只好自己再推一遍,先贴出是整数的快速幂算法:
int Pow(int x, int y){
    int ans, b, i;
    if(y == 0)return 1;
    if(y < 0) return -1;
    for(i = y, b = 0; i > 0; i >>= 1) b++;
    for(ans = 1, i = b - 1; i >= 0; i--){
        if(y & (1 << i)){
            ans = ans * ans * x;
        }else{
            ans = ans * ans;
        }
    }
    return ans;
}

下面是矩阵类的代码:
#include<iostream>
#include<math.h>
using namespace std;

template <typename T>
class Matrix{
//class Matrix by felix021
//usage: Matrix<typename> varname(row, col[, mod]);
//如果typename不是short/int/long/long long
//那么取模运算需要包含math.h头文件(使用fmod函数)
//取模运算的代码仅在 + - * 三个函数内.

    public:
        T **data;
        Matrix *temp;
        int row, col, modnum;
        //无参构造函数,没有给定矩阵的大小时,输出错误提示
        Matrix(){
            printf("Size parameter invalid!\n");
        }
        //构造函数, r是行数,c是列数,若给定m>0, 则在求和/差/积时对m取模
        Matrix(int r, int c, int m = 0){
            int i;
            temp = NULL;
            row = r, col = c;
            data = new T *[r]; //分配一个 T *[r]类型的数组给data
            for(i = 0; i < r; i++){
                data[i] = new T [c];  //给data下的每个指针分配内存, 类型为 T [c]
            }
            if(m > 0) modnum = m;
            else modnum = 0;
        }
        //析构函数,用于释放分配的内存
        ~Matrix(){
            int i;
            if(temp != NULL){delete temp; }
            for(i = 0; i < row; i++) delete data[i];
            delete data;
        }
        //赋值,给第i行第j列的元素赋值为 a
        void assign(int i, int j, const T a){
            if(i < 0 || i >= row || j < 0 || j >= col) {
                printf("Invalid row/column number.\n");
                return;
            }else{
                data[i][j] = a;
            }
        }
        //取出元素,给第i行第j列的元素的值
        T at(int i, int j){
            if(i < 0 || i >= row || j < 0 || j >= col) {
                printf("Invalid row/column number.\n");
                return;
            }else{
                return data[i][j];
            }
        }
        //重载<<运算符,可用cout输出
        friend inline ostream & operator<<(ostream &os, const Matrix &a){
            int i, j;
            for(i = 0; i < a.row; i++){
                for(j = 0; j < a.col - 1; j++){
                    os << a.data[i][j] << " ";
                }
                os << a.data[i][j];
                os << endl;
            }
            return os;
        }
        //重载>>运算符,可用cin输入
        friend inline istream & operator>>(istream &is, const Matrix &a){
            int i, j;
            for(i = 0; i < a.row; i++){
                for(j = 0; j < a.col; j++){
                    is >> a.data[i][j];
                }
            }
            return is;
        }
        //重载=运算符,要求两个矩阵的大小相同
        Matrix & operator = (const Matrix &a){
            int i, j;
            if(row != a.row || col != a.col){
                printf("Unmatch Matrix!\n");
                return *this;
            }
            for(i = 0; i < row; i++){
                for(j = 0; j < col; j++){
                    data[i][j] = a.data[i][j];
                    if(modnum > 0){
                        data[i][j] = mod(data[i][j]);
                    }
                }
            }
            return *this;
        }
        //重载+运算符,要求两个矩阵的大小相同
        Matrix & operator + (const Matrix &a){
            int i, j;
            if(row != a.row || col != a.col){
                printf("Unmatch Matrix!\n");
                return *this;
            }
            if(temp != NULL) delete temp;
            temp = new Matrix(row, col, modnum);
            for(i = 0; i < row; i++){
                for(j = 0; j < col; j++){
                    temp->data[i][j] = data[i][j] + a.data[i][j];
                    if(modnum > 0){
                        temp->data[i][j] = mod(temp->data[i][j]);
                    }
                }
            }
            return *temp;
        }
        //重载-运算符,要求两个矩阵的大小相同
        Matrix & operator - (const Matrix &a){
            int i, j;
            if(row != a.row || col != a.col){
                printf("Unmatch Matrix!\n");
                return *this;
            }
            if(temp != NULL) delete temp;
            temp = new Matrix(row, col, modnum);
            for(i = 0; i < row; i++){
                for(j = 0; j < col; j++){
                    temp->data[i][j] = data[i][j] - a.data[i][j];
                    if(modnum > 0){
                        temp->data[i][j] = mod(temp->data[i][j] + modnum);
                    }
                }
            }
            return *temp;
        }
        //重载*运算符,要求矩阵a的列数等于b的行数
        Matrix & operator * (const Matrix &a){
            int i, j, k;
            T tmp;
            if(col != a.row){
                printf("Unmatch Matrix!\n");
                return *this;
            }
            if(temp != NULL) delete temp;
            temp  = new Matrix(row, a.col, modnum);
            for(i = 0; i < row; i++){
                for(j = 0; j < a.col; j++){
                    tmp = 0;
                    for(k = 0; k < a.row; k++){
                        tmp += data[i][k] * a.data[k][j];
                        if(modnum > 0){
                            tmp = mod(tmp);
                        }
                    }
                    temp->data[i][j] = tmp;
                }
            }
            return *temp;
        }
        //重载^运算符,要求矩阵的列数等于行数
        Matrix & operator ^ (const int a){
            int i, j;
            if(row != col){
                printf("No n*n matrix!\n");
                return *this;
            }
            if(a < 0){
                printf("Invalid a(%d)!\n", a);
                return *this;
            }
            if(temp != NULL) delete temp;
            temp  = new Matrix(row, col, modnum);
            for(i = 0; i < row; i++){ //单位方阵
                for(j = 0; j < col; j++){
                    if(i == j )temp->data[i][j] = 1;
                    else temp->data[i][j] = 0;
                }
            }
            for(i = a, j = 0; i > 0; i >>= 1) j++;
            for(i = j - 1; i >= 0; i--){
                if(a & (1 << i)){
                    *temp = *temp * *temp * (*this);
                }else{
                    *temp = *temp * *temp;
                }
            }
            return *temp;
        }
        //重载求模函数,对int, long, long long, float, double都有效
        int mod(int i){
            return i % modnum;
        }
        long mod(long i){
            return i % modnum;
        }
        long long mod(long long i){
            return i % modnum;
        }
        float mod(float i){
            return fmod(i, modnum);
        }
        double mod(double i){
            return fmod(i, modnum);
        }
};

int main(){
    int i, j;
    Matrix <int> t(2, 2), t1(2, 2), t2(2, 2), a(5, 10, 7), b(10, 5), c(5, 5);
    cout << "(cin/cout test)\nPlease input 2 (2*2) Matrix: \n";
    cin >> t;
    cout << t << endl;
    t1 = t * t * t * t * t;
    t2 = t ^ 5;
    cout << "t^5 = \n" << t1 << endl << t2 << endl;
    getchar();
    
    for(i = 0; i < 5; i++){
        for(j = 0; j < 10; j++){
            a.assign(i, j, 1);
            b.assign(j, i, 2);
        }
    }
    c = a * b;
    cout << "Matrix a:" << endl << a << endl;
    cout << "Matrix b:" << endl << b << endl;
    cout << "Matrix a * b:" << endl << c << endl;
    getchar();
    return 0;
}
Aug 15

Felix's Web Proxy For FTP V1.0 不指定

felix021 @ 2008-8-15 06:34 [IT » 网络] 评论(1) , 引用(0) , 阅读(6633) | Via 本站原创
    花了点时间写了个web proxy for ftp, 也就是用php的扩展库,实现通过http的80端口来访问ftp服务器。这样做的动机是为了越过武大教育网对ftp端口的限制(也就是说以后我可以在外面通过http来访问武大的内网FTP了)。不过为了防止搜索引擎通过这个脚本抓取,所以即使是匿名可以登陆的服务器也需要先登陆后才能使用。当然,如果某人需要,简单修改一下这个脚本就能够突破这个限制。另外一个限制是可以下载的文件大小上限,因为具体的实现是通过服务器把ftp上的文件下载到本地,然后再发送给浏览器,如果可以任意下载大文件,那么服务器压力会很大,而且很容易失败。不过在必要的时候还是得下载一些大文件,所以另外设置了一个Privilege Password,这样就可以不受文件大小的限制了。因为时间关系只做了目录浏览和下载,没有做其他功能(重命名、上传、删除等)。等以后有时间了重写吧,得把那些东西写到一个Class里面去才行。不过其实有现成的ftp_class, 可以在phpWebFtp的包里面找到那个文件,我的一个关键函数就是从那个class里面提取出来的,而phpWebFtp是使用GPL2作为其分发协议的,所以我这个也用GPL2吧~

废话说多了,提供下载吧。
下载文件 (已下载 1193 次)
Aug 14
linux下的shell用惯了,在windows的cmd下面总是习惯打上ls和clear这样的命令,然后就会很郁闷。于是又把cygwin装上了。本来打算等过了15天包换期再说的,但是实在受不了windows垃圾的cmd了。。不过郁闷的是,在cygwin下面居然没有clear这个命令,那个汗啊。不过可以自己为它写一个名为clear的脚本,保存在/bin下面。其实也很简单,就是调用cmd, 运行一个cls命令。

#vim /bin/clear

---以下为clear的内容---
#!/bin/bash
cmd /c cls
---clear结束---

保存退出,OK,以后就可以尽管用clear啦!

然后把c:\cygwin\bin加到系统的path变量里面去,以后在cmd下面也可以用ls这些命令了。不过在cmd下面还是不能用clear这个命令——其实这个更容易解决了,在c:\cygwin\bin下面加入一个clear.bat, 内容为@cls, 保存,OK!
Tags: ,
Aug 14

ZZ: GDB使用教程 不指定

felix021 @ 2008-8-14 02:52 [IT » 程序设计] 评论(0) , 引用(0) , 阅读(3764) | Via 本站原创
from http://hi.baidu.com/buddhist_byr/blog/item/214bf995f4a30a48d0135e38.html
GDB使用教程


忘了从哪里找到这篇文章的,看了之后对于像我这样的初学者有很大帮助.为了今后参考,在这里把这篇文章贴出来.也希望更多初学者

能够从中受益.这篇文章所有版权归原作者所有.


GDB概述
————

GDB 是GNU开源组织发布的一个强大的UNIX下的程序调试工具。或许,各位比较喜欢那种图形界面方式的,像VC、BCB等IDE的调试,

但如果你是在 UNIX平台下做软件,你会发现GDB这个调试工具有比VC、BCB的图形化调试器更强大的功能。所谓“寸有所长,尺有所

短”就是这个道理。

一般来说,GDB主要帮忙你完成下面四个方面的功能:
Tags:
Aug 14

Priority Queue - Binary Heap 不指定

felix021 @ 2008-8-14 02:51 [IT » 程序设计] 评论(0) , 引用(0) , 阅读(4058) | Via 本站原创
from http://www.programfan.com/blog/article.asp?id=20319

没什么好写的,都在《Introduction to Algorithms》里面,老外的教材就是厉害,不像我们国内,计算机类书籍不少,优点是:一、贵,二、内容不全,有时候要学个什么东西查一本书还不行,唉,写书都是为了赚钱,哪有什么好书。

前些天我借了本《算法设计与分析》的书,读着读着发觉怎么好像读过,一想,哦,是这本书里看过,那时候觉得看英文累,印象深刻,原来写中文书还有这么一招,直接翻译来就OK。老外不喜欢东搞西搞,一本书,比砖头还厚,把想说的全说完,全部搞定,简单,容易查阅,买一本一辈子带身边。哈哈,我刚好有一本,影印版的。

这个Priority Queue原来学《数据结构》的时候老师没讲,自己也没重视,现在翻出来学,发现很有用的。
Aug 13
中文版 "C C++ 语言 函数 参考手册.chm"  122KB
内有C/C++库函数以及STL容器的使用说明。
源自 http://cppreference.com/
Aug 12
C++重载operator的示例 By Felix021

以下示例中定义了一个class test, 重载了<, +, +=, =, ==, <<, >>等符号:
Tags: , ,
分页: 66/103 第一页 上页 61 62 63 64 65 66 67 68 69 70 下页 最后页 [ 显示模式: 摘要 | 列表 ]