Summer day 13 & 14

这两天没做什么题u，主要在看ppt和玩，以及越狱升级。所以也就只能搞搞0725训练赛的题目了。

A. UVA 116 Unidirectional TSP

会玩的旅行商系列。
走法如图。

求最左到最右的最小花费。
DP, 字典序需注意。
模仿高玩代码。

#include <iostream>
#include <cstdlib>
#include <cstring>
#include <cstdio>

using namespace std;

int maps[11][101];
int smap[11][101];
int fath[11][101];

int main()
{
    int n,m;
    while (~scanf("%d%d",&n,&m)) {
        for (int i = 1 ; i <= n ; ++ i)
        for (int j = 1 ; j <= m ; ++ j)
            scanf("%d",&maps[i][j]);

        memset(smap, 0, sizeof(smap));

        for (int i = m ; i >= 1 ; -- i)
        for (int j = 1 ; j <= n ; ++ j) {
            smap[j][i] = smap[j][i+1]+maps[j][i];
            fath[j][i] = j;
            if (j > 1 && smap[j][i] >= smap[j-1][i+1]+maps[j][i]) {
                smap[j][i] = smap[j-1][i+1]+maps[j][i];
                fath[j][i] = j-1;
            }
            if (j == n && smap[j][i] >= smap[1][i+1]+maps[j][i]) {
                smap[j][i] = smap[1][i+1]+maps[j][i];
                fath[j][i] = 1;
            }

            if (j < n && smap[j][i] > smap[j+1][i+1]+maps[j][i]) {
                smap[j][i] = smap[j+1][i+1]+maps[j][i];
                fath[j][i] = j+1;
            }
            if (j == 1 && smap[j][i] > smap[n][i+1]+maps[j][i]) {
                smap[j][i] = smap[n][i+1]+maps[j][i];
                fath[j][i] = n;
            }
        }

        int spa = 1;
        for (int i = 2 ; i <= n ; ++ i)
            if (smap[spa][1] > smap[i][1])
                spa = i;

        int min = smap[spa][1];
        for (int i = 1 ; i <= m ; ++ i) {
            if (i < m) printf("%d ",spa);
            else printf("%d\n%d\n",spa,min);
            spa = fath[spa][i];
        }
    }
    return 0;
}

很荣幸去讲了一道题，就是B题，小数进制转换。

B. POJ 1131 Octal Fractions

可以用Java大小数类来做，我不会。
于是
强行上C。
思路：（直接复制ppt）
不能用浮点运算
对于0.d1 d2 d3 … dk，
十进制数为ans=∑▒〖8^(−n)∗����〗
计算方法：（dk/8 +dk-1）/8+……
用num表示余数，从八进制最后一位开始读取，类似于无小数的进制转换，每读取一个数更新一下要输出的数组。
有小数进制转换的模板。
注释说明了一切。

#include<cstdio>
#include<cstring>

const int maxN 100

int main()
{
    char num_8[maxN];
    while(~scanf("%s",num_8))
    {
        memset(num_10, 0, sizeof(num_10));
        int len =0;     //表示十进制长度
        int j;
        for(int i=strlen(num_8)-1;i>1;i--)  从八进制中的最后一位开始做除法
        {   int num = num_8[i]-'0';     //读取新位数
            for( j=0;(j<len||num);j++)  //如果没有到10进制的最后一位或者余数不为0就继续运行
            {
                int temp = num*10+(j<len?num_10[j]-'0':0);
                num_10[j]=temp/8+'0';
                num=temp%8;
            }
            len=j;      //更新十进制长度
            for(int k = 0; k<len; k++)
                printf("%c ",num_10[k]);
            printf("\n");
        }
        num_10[j]='\0';
        printf("%s [8] = 0.%s [10]\n",num_8,num_10);
    }
    return 0;
}

C是入门经典原题，
D是求点集最短连线。
高玩略去。

E. CodeForces 688D

题意：给定x，k,和n个ci。你可以知道x%ci,问是否能确定x%k。
解法：自豪地介绍（误）%*￥#@！——中国剩余定理。
道理大家都（几乎不）懂，经过几个懂的人抽象，题意变成了这样：
题意‘：C 数组整体的最小公倍数 lcm(c) 是否是 K 的倍数，如果是，则能确定输出 yes，否则输出 no.

gcd（）和lcm（）就是其核心了：

LL gcd(LL a, LL b)
{
    return b?gcd(b, a%b):a;
}

LL lcm(LL a, LL b)
{
    return a*b/gcd(a, b);
}

然后经验告诉我开个（够大）的数组。

//RE:arrary too big
#include<cstdio>
#define LL long long

const int maxn = 1000005; 

int main()
{
    LL n, k;
    LL c[maxn];
    scanf("%ld%ld",&n, &k);
    for(int i = 0; i<n; i++)
        scanf("%ld", &c[i]);
    LL temp = c[0];
    for(int i = 1; i<n; i++)
    {   temp = lcm(temp, c[i]);
    }
    if(temp%k == 0)
        printf("Yes\n");
    else
        printf("No\n");
    return 0;
}

= =。
所以还是一次读一个计算一个比较老实靠谱。