C++博客-yesg-随笔分类-算法

扩展欧几里德和中国剩余定理

nk_ysg — Sat, 02 Jun 2012 06:31:00 GMT

/*
*/
#include
using namespace std;

/*
注意到对于gcd(a,b) = d 我们对(a, b)用欧几里德辗转相除会最终得到
(d, 0)此时对于把a =d, b = 0 带入a*x + b*y = d，显然x = 1，y可以为任意值，
这里y可以为任意值就意味着解会有无数个。我们可以用a = d, b = 0的情况逆推出来
任何gcd(a, b) = d 满足a*x + b*y = d的解。如果x0, y0是b*x + (a%b)*y = d 的解，
那么对于a*x + b*y = d的解呢？
b*x0 + (a%b)*y0 = d => b*x0 + (a - [a/b]*b)*y0 = a*y0 + b*(x0 - [a/b]*y0)，
所以a*x + b*y = d的解x1 = y0， y1 = x0 - [a/b]*y0; 这样我们可以程序迭带了。
*/
int extEuclid(int a, int b, int &x, int &y)
{
    if (b == 0)
    {
        x = 1;
        y = 0;
        return a;
    }
    int d = extEuclid(b, a % b, x, y);
    int iTemp = x;
    x = y;
    y = iTemp - (a / b)* y;
    return d;
}
//解同余方程ax = b(mod n) (返回最小的正数x)

int modularLinearEquation(int a, int b, int n)
{
    //等价于求ax + cn = b;
    //先求a*x1 + c1*n = gcd(a, n)
    int x, y, d;
    d = extEuclid(a, n, x, y);
    if (b % d != 0)
       return -1;
    x = x * (b / d);
    x = (( x % n) + n) % n;
    return x;
}

//中国剩余定理，推导都是数学
int solModularEquations(int b[], int m[], int k)
{

    int iTemp;
    int y;
    int result;

    int M = 1;
    for (int i = 0; i < k; i++)
      M *= m[i];

    result = 0;
    for (int i = 0; i < k; i++)
    {
      iTemp = M / m[i];
      y = modularLinearEquation(iTemp, 1, m[i]);
      result = (result + b[i] * iTemp * y) % M;
    }
    return result;
}

int main()
{
    int x, y , d;
    d = extEuclid(1001, 767, x, y);
    cout << x << endl;
    cout << y << endl;
    cout << d << endl;
    cout << "1001 * x + 767 * y = " << (1001 * x + 767 * y) << endl;
    cout << modularLinearEquation(3, 2, 100) << endl;
    return 0;
}

nk_ysg 2012-06-02 14:31 发表评论

数字三角形(记忆化搜索和递归)

nk_ysg — Fri, 07 Oct 2011 17:24:00 GMT

/*
给定一个由n行数字组成的数字三角形如下图所示。
试设计一个算法，计算出从三角形的顶至底的一条路径，使该路径经过的数字总和
对于给定的由n行数字组成的数字三角形，
编程计算从三角形的顶至底的路径经过的数字和的最大值

*/
#include <stdio.h>
#include <string.h>

const int MAXN = 100 + 10;

int a[MAXN][MAXN],d[MAXN][MAXN];
int n;

inline int max(int a,int b) { return a > b ? a : b;}

int dfs(int i,int j)
{
    if (d[i][j] >= 0) return d[i][j];
    return d[i][j] = a[i][j] + (i == n ? 0 : max(dfs(i+1,j),dfs(i+1,j+1)));
}

int main()
{
    int i;
    int j;
    freopen("1.txt","r",stdin);
    while (scanf("%d",&n) != EOF)
    {
        memset(d,-1,sizeof(d));
        for (i = 1; i <= n; ++i)
        {
            for (j = 1; j <= i; ++j)
                scanf("%d",&a[i][j]);
        }
        printf("%d\n",dfs(1,1));
    }
    return 0;
}

#include <stdio.h>

const int MAXN = 100 + 10;

int a[MAXN][MAXN],d[MAXN][MAXN];

int n;

inline int max(int a, int b) { return a > b ? a : b; }

int main()
{
    int i,j;
    freopen("1.txt","r",stdin);
    while (scanf("%d",&n) != EOF)
    {
        for (i = 1; i <= n; ++i)
            for (j = 1; j <= i; ++j)
                scanf("%d",&a[i][j]);

        for (j = 1; j <= n; j++) d[n][j] = a[n][j];
        for (i = n - 1; i >= 1; i--)
            for (j = 1; j <= i; j++)
                d[i][j] = a[i][j] + max(d[i+1][j],d[i+1][j+1]);
        printf("%d\n",d[1][1]);
    }

    return 0;
}

nk_ysg 2011-10-08 01:24 发表评论

分数拆分

nk_ysg — Fri, 07 Oct 2011 09:04:00 GMT

1 /*
2 输入正整数k,找到所有的正整数x>=y,使得1/k = 1/x + 1/y.
3 样例输入:
4 2
5 12
6 样例输出
7 2
8 1/2=1/6+1/3
9 1/2=1/4+1/4
10
11 8
12 1/12 = 1/156 + 1/13
13 1/12 = 1/84 + 1/14
14 1/12 = 1/60 + 1/15
15 1/12 = 1/48 + 1/16
16 1/12 = 1/36 + 1/18
17 1/12 = 1/30 + 1/20
18 1/12 = 1/28 + 1/21
19 1/12 = 1/24 + 1/24
20 */
21 #include <stdio.h>
22
23 const int MAXN = 1000 + 10;
24 int x[MAXN],y[MAXN];
25
26 int main()
27 {
28     int k;
29     int ans;
30     int i;
31     while (scanf("%d",&k) != EOF)
32     {
33         ans = 0;
34         for (i = k + 1; i <= 2*k; ++i)
35         {
36             if ((k * i) % (i-k) == 0) { x[ans] = (k*i)/(i-k),y[ans++] = i; }
37         }
38         printf("%d\n",ans);
39         for (i = 0; i < ans; ++i)
40         {
41             printf("1/%d = 1/%d + 1/%d\n",k,x[i],y[i]);
42         }
43     }
44     return 0;
45 }
46

nk_ysg 2011-10-07 17:04 发表评论

拓扑排序(感谢RoBa和NOCOW)

nk_ysg — Thu, 06 Oct 2011 11:19:00 GMT

Genealogical tree. 家谱树 Time Limit: 1 second Memory Limit: 16M

Background 大意：火星人的种族关系比较混乱。宗族管理委员会进行开会发言等活动时，又想遵守先长辈，后晚辈的顺序。于是有下面的问题：

Problem 编写一个程序对给定的人员，决定一个先后顺序，这个顺序必须遵守先长辈，后晚辈的原则。

Input 首行为N， 1 <= N <= 100 ,N为总人数。根据百年传统，给人员用自然数编号为1至N。以下的N行，第I行为第I个人的子孙列表，子孙的顺序是任意的，用空格隔开，且以0为结束。子孙列表可以是空的。

Output 在一行内输出发言的顺序。如有多个顺序满足条件，则输出任一个。至少会有一个满足的顺序的。

Sample Input

5 0 4 5 1 0 1 0 5 3 0 3 0

Sample Output

2 4 5 3 1

/* URAL 1022 Genealogical tree
* 900803 09:04:18 21 Aug 2005 RoBa 1022 C++ Accepted 0.001 190 KB
* 拓扑排序
*/

#include <stdio.h>
#include <stdlib.h>
#include <string.h>

const int MAX = 101;
int map[MAX][MAX],inde[MAX],taken[MAX];

int main()
{
    int i,n,tmp,j;
    while (scanf("%d",&n)!=EOF)
    {
        memset(map,0,sizeof(map));
        memset(inde,0,sizeof(inde));
        for (i = 1 ; i <= n ; i++)
            while (scanf("%d",&tmp),tmp)
            {
                map[i][tmp] = 1;
                ++inde[tmp];
            }
        while (1)
        {
            for (i = 1 ; i <= n ; i++)
                if (inde[i] == 0 && taken[i] == 0) break;
            if (i > n) break;
            for (j = 1 ; j <= n ; j++)
                if (map[i][j])
                {
                    map[i][j] = 0;
                    --inde[j];
                }
            taken[i] = 1;
            printf("%d ",i);
        }
        printf("\n");
    }
    return 0;
}

nk_ysg 2011-10-06 19:19 发表评论