C++博客-zoyi&&菠菜-随笔分类-数学

大素数的检验

zoyi — Tue, 23 Sep 2008 05:26:00 GMT

摘要: 阅读全文

zoyi 2008-09-23 13:26 发表评论

spoj Sum of one-sequence

zoyi — Sat, 12 Jul 2008 04:08:00 GMT

在poi的列表上是easy。。。。
我倒。。这也easy。。。。。。
这道题多亏了javaman。。呵呵。。这一次。。他讲的挺清楚的。。。谢谢他一下

题目大意：
根据输入，序列的长度，序列和，求这样一串序列
序列满足以下情况：

for any k, such that 1 <= k < n : |a_k - a_k₊₁| = 1 and
a₁ = 0

如果没有如此序列则输出"No"

我们假设n-->序列长度，S-->序列和；

S=sigma a[i](1<=i<=n)
令bi=a[i+1]-a[i] (1<=i

S=sigama b[i]+sigama a[i] (1<=i S= sgama (n-i)b[i]; (1<=i
b[i]-->-1 or 1

S=sigama (n-i)(b[i]+1)-sigama(n-i)
=2*sigama ((n-i)(b[i]+1)/2) - n*(n-1)/2

d[i]=(b[i]+1)/2-->0 or 1
根据上式得：
no-->
1: S+n*(n-1)/2 为奇数
2: abs(S)>n*(n-1)/2 //important ....我就是这里一直错...绝对值...没考虑....

#include<iostream>
#include<math.h>
int N,S;
int main()
{
    int std,i,d,pre;
    int T;
    scanf("%d",&T);
    for(int j=1;j<=T;j++)
    {
        scanf("%d%d",&N,&S);
        std=(N*(N-1))/2;
        if((S+std)%2||abs(S)>std)
        {
            printf("No\n");
            if(j!=T)printf("\n");
            continue;
        }
        printf("0\n");
        for(pre=0,std+=S,std/=2,i=1;i<N;i++)
        {
            if(std>=N-i)
            {
                d=1;
                std-=N-i;
            }
            else d=0;
            printf("%d\n",pre=2*d-1+pre);
        }
        if(j!=T)printf("\n");
    }
    return 0;
}

zoyi 2008-07-12 12:08 发表评论

pku 2917 Diophantus of Alexandria

zoyi — Wed, 09 Jul 2008 02:41:00 GMT

一道数学题。弱弱不会做。。

开始一直钻牛角尖。。用暴力把数据跑了一下。然后找规律。。发现好像是找n的互质的质因数对数。。
小的数据用手算是可以算出来的。。但是不知道怎么用程序来算。。

后面参考pc的。。。。知道了是怎么回事。。
我想以后如果有需要求质因数对数的个数。。可以用这种方法反过来求。。

下面是代码：

#include<iostream>
int N;
int solve(int N)
{
    int i,ans=1,t;
    for(i=2;i*i<=N;i++)
    {
        t=0;
        while(!(N%i))
        {
            t++;
            N/=i;
        }
        ans*=2*t+1;
    }
    if(ans==1||N!=1)ans*=3;
    return ans;
}
int main()
{
    int T,cas,id;
    scanf("%d",&T);
    for(cas=1;cas<=T;cas++)
    {
        scanf("%d",&N);
        id=solve(N);
        if(N==1)id=1;
        printf("Scenario #%d:\n%d\n\n",cas,(id+1)/2);
    }
    return 0;
}

zoyi 2008-07-09 10:41 发表评论

平面中的直线问题

zoyi — Thu, 10 Apr 2008 15:36:00 GMT

2008年4月10日23:52:47
第一个问题：
由平面中的n条直线确定的最大区域数Ln
L0=1；
Ln=Ln-1+n(n>0);
Ln=1+n*(n+1)/2;

第二个问题：
是平面直线的变形问题，用弯曲的线来代替直线，每一个弯曲线含有一个“锯齿形的转角”，同样确定平面区域的最大个数Zn
（我们把一条弯曲折线抽象为两条，但是合并了某些区域）
Zn=L2n-2*n=2*n*n-n+1;(n>=0);

第三个问题：
就是一下这个问题：
count the regions

若当前有n-1条边，那么在往里面加一条边，这条新加的边最多和以前的边有9*(n-1)个交点，
那么会添加 9*（n-1）+1个面
这条规律对于上面两种也是用，加x个点，那么就会添加x+1个面

那么总结下来 Xn=Xn-1+9*(n-1)+1
即Xn=(9/2)*n*n-(7/2)*n+1;

以下是我的代码：

#include<stdio.h>
int main()
{
    long long  n;
    long long  ans;
    while(scanf("%lld",&n)!=EOF){
        ans=9*n*n-7*n+2;
        ans/=2;
        printf("%lld\n",ans);
    }
    return 0;
}

zoyi 2008-04-10 23:36 发表评论