e81086713e446d36f62b2aa2a3502b5eb155

java杂家 -凯发k8网页登录

杂七杂八。。。一家之言

凯发k8网页登录-凯发天生赢家一触即发官网 凯发k8网页登录首页
  141 posts :: 1 stories :: 174 comments :: 0 trackbacks

置顶随笔 #

如题:求连续正整数使得其和为给定的一个正整数
下面给出我的解法,几乎可以一步到位求出来
实现代码如下:
/**
*author: koth (
http://weibo.com/yovn)
*date:  2011-12-01
*/
#include 
<stdlib.h>
#include 
<stdio.h>
#include 
<stdint.h>

int solve(int y,int& x){
    
int m=0;
    
int t=y;
    
if(y<=0){
        x
=y;
        
return 1;
    }
    
while((t&1)==0){
        m
=1;
        t
=t>>1;
    }
    
if(m==32){
        x
=y;
        
return 1;
    }
    
int lastk=32;
    
for(;lastk>m1;lastk--){
        
if(y &(1<<(lastk-1))){
            
            
break;
        }
            
    }

    
//its a number st. exp(2,k)
    if(lastk==(m1)){
        x
=y;
        
return 1;
    }
    
int k=1<<(m1);
    
int b=(y>>m)-(1<<(lastk-m-1));

    x
=(1<<(lastk-m-2))(b1-k)/2;

    
if(x<=0){
        k
=k-1-((0-x)<<1);
        x
=0-x1;
    }
    
    
return k;

}

int main(int argc,char* argv[]){
    
if(argc<=1){
        fprintf(stdout,
"usage:%s number\n",argv[0]);
        
return 0;
    }
    
int y=atoi(argv[1]);
    
int x=0;
    
int k=solve(y,x);
    fprintf(stdout,
"%d=",y);
    
for(int i=0;i<k;i){
        fprintf(stdout,
"%d",xi);
        
if(i<(k-1)){
            fprintf(stdout,
"");
        }
    }
    fprintf(stdout,
"\n");
    
return 0;
}
posted @ doubleh 阅读(1749) | |  

原题:

command network

description

after a long lasting war on words, a war on arms finally breaks out between littleken’s and knuthocean’s kingdoms. a sudden and violent assault by knuthocean’s force has rendered a total failure of littleken’s command network. a provisional network must be built immediately. littleken orders snoopy to take charge of the project.

with the situation studied to every detail, snoopy believes that the most urgent point is to enable littenken’s commands to reach every disconnected node in the destroyed network and decides on a plan to build a unidirectional communication network. the nodes are distributed on a plane. if littleken’s commands are to be able to be delivered directly from a node a to another node b, a wire will have to be built along the straight line segment connecting the two nodes. since it’s in wartime, not between all pairs of nodes can wires be built. snoopy wants the plan to require the shortest total length of wires so that the construction can be done very soon.

input

the input contains several test cases. each test case starts with a line containing two integer n (n ≤ 100), the number of nodes in the destroyed network, and m (m ≤ 104), the number of pairs of nodes between which a wire can be built. the next n lines each contain an ordered pair xi and yi, giving the cartesian coordinates of the nodes. then follow m lines each containing two integers i and j between 1 and n (inclusive) meaning a wire can be built between node i and node j for unidirectional command delivery from the former to the latter. littleken’s headquarter is always located at node 1. process to end of file.

output

for each test case, output exactly one line containing the shortest total length of wires to two digits past the decimal point. in the cases that such a network does not exist, just output ‘poor snoopy’.


一开始没仔细读题,一看以为是最小生成树呢,结果krusal算法上去wa了,prim算法也wa,修修改改一直wa,第二天发现本题是要在有向图上面构造最小树形图。

按照著名的zhu-liu算法,仔细实现了一边,终于ac了。
按照我的理解总结下该算法,该算法对每个结点,除根节点外寻找最小入边,
1)如果这些入边不构成环,那么容易证明这些边构成最小树形图。
证明:设加上根节点r一共n个点,那么一共有n-1条边,证明从r能到每个点,若存在一点v,使得从r到v没有路径,那么,假设从v反向回退必然构成环,因为每个点除了r都有入边,如果不构成环,该路径可以无穷大。
2)如果存在环,我们把环收缩成一个点,更新相应的入边和出边,得到新的图g',使得原问题在g'中等效:
怎么收缩呢?
假设我们把环收缩成环上的任意一点v,所有进环的边和出环的边自动变成v的边(如果已有,取长度最短的),其余点标记为删除,更新不在环上的所有点进入该环的长度cost为cost-cost(prev[x],x);其中点x为进入环的边在环上的端点。出边保持不变。

这里为什么这么更新?因为这样更新使得我们的算法在新图上是等效的。任何环的解决后意味着在新图里面得为改收缩后的点寻找新的入边,而实际的花费应该是新的入边减去原有的入边长度,我们的算法在找到环的时候就把环上所有的边的长度计算在花费内了.而对出边是没有影响的。


到这算法的框架基本出来了。当为某点没找到入边的时候,意味着无解。为了加快无解的侦测,我们先运行一遍dfs搜索,假如从根节点出发,可触及的节点数小于n-1(不含r)则意味着无解。反之,肯定有解。
为什么?
因为如果可触及数小于n-1,意味着某点是不可触及的,也就是原图不是弱连通。对该点来说不存在从r到它的路径。反之,从r到某点都有一条路径,沿着该路径就能找到结点的入边。


第二个问题是,如何快速侦测环呢?
我使用了一个不相交集。回忆krusal的算法实现里面也是使用不相交集来避免找产生环的最小边。

下面是我的代码:

// 3164.cpp : defines the entry point for the console application.
//

#include 
<iostream>
#include 
<cmath>


using namespace std;

typedef 
struct _point
{

    
double x;
    
double y;

    
double distanceto(const struct _point& r)
    {
          
return sqrt((x-r.x)*(x-r.x)(y-r.y)*(y-r.y));

    }

}point;



const int max_v=100;
const int max_e=10000;
const double no_edge=1.7976931348623158e 308;

point vertexes[max_v]
={0};
int parents[max_v]={0};
int ranks[max_v]={0};
double g[max_v][max_v]={0};
bool visited[max_v]={0};
bool deleted[max_v]={0};
int prev[max_v]={0};

int nvertex=0;
int nedge=0;




int u_find(int a)
{
    
if(parents[a]==a)return a;
    parents[a]
=u_find(parents[a]);
    
return parents[a];
}
void u_union(int a,int b)
{
    
int pa=u_find(a);
    
int pb=u_find(b);
    
if(ranks[pa]==ranks[pb])
    {
        ranks[pa]
;
        parents[pb]
=pa;
    }
else if(ranks[pa]<ranks[pb])
    {
        parents[pa]
=pb;
    }
    
else
    {
        parents[pb]
=pa;
    }
}

void dfs(int v,int& c)
{

    visited[v]
=true;
    
for(int i=1;i<nvertex;i)
    {
        
if(!visited[i]&&g[v][i]<no_edge)
        {
            c
=1;

            dfs(i,c);
        }

    }

}



void docycle(int s,int t,double& cost)
{
    memset(visited,
0,sizeof(bool)*nvertex);
    
int i=s;
    
do
    {
        visited[i]
=true;
        cost
=g[prev[i]-1][i];
        
//cout<<"from "<<(prev[i]-1)<<" to "<
        i=prev[i]-1;

    }
while(i!=s);


    
do
    {
        

        
for(int k=0;k<nvertex;k)
        {
            
if(!deleted[k]&&!visited[k])
            {
            
                
if(g[k][i]<no_edge)
                {
                    
if(g[k][i]-g[prev[i]-1][i]<g[k][s])
                    {


                        g[k][s]
=g[k][i]-g[prev[i]-1][i];
                        
//cout<<"1.update ["<
                    }
                    

                }
                
if(g[i][k]<no_edge)
                {
                    
if(g[i][k]<g[s][k])
                    {

                        g[s][k]
=g[i][k];
                        
//cout<<"2.update ["<
                    }
                    
                }
            }
        }


        
if(i!=s)
        {
            deleted[i]
=true;
            
//cout<<"mark "<
        }
        i
=prev[i]-1;


    }
while(i!=s);



}




int main(void)
{



    
    
while(cin>>nvertex>>nedge)
    {

        
int s,t;

        
int nv=0;
        
bool cycle=0;
        
double cost=0;
        memset(vertexes,
0,sizeof(vertexes));
        memset(visited,
0,sizeof(visited) );

        memset(deleted,
0,sizeof(deleted));
        memset(g,
0,sizeof(g));
        memset(prev,
0,sizeof(prev));


        memset(ranks,
0,sizeof(ranks));

        memset(parents,
0,sizeof(parents));

        
for(int i=0;i<nvertex;i)
        {

            cin
>>vertexes[i].x>>vertexes[i].y;
            parents[i]
=i;
            
for(int j=0;j<nvertex;j)
            {
                g[i][j]
=no_edge;
            }

        }
        
for(int i=0;i<nedge;i)
        {

            cin
>>s>>t;
            
if(t==1||s==t)continue;
            g[s
-1][t-1]=vertexes[s-1].distanceto(vertexes[t-1]);

        }



        dfs(
0,nv);
        
if(nv<nvertex-1)
        {

            cout
<<"poor snoopy"<<endl;
            
continue;
        }



        
do {
            cycle
=false;
            
for(int i=0;i<nvertex;i){parents[i]=i;}
            memset(ranks,
0,sizeof(bool)*nvertex);
        

            
for (int i = 1; i < nvertex; i) {
                
double minimum=no_edge;
                
if(deleted[i])continue;
                
for(int k=0;k<nvertex;k)
                {
                    
if(!deleted[k]&&minimum>g[k][i])
                    {
                        prev[i]
=k1;
                        minimum
=g[k][i];
                    }
                }
                
if(minimum==no_edge)
                {
                

                    
throw 1;
                }
                
if (u_find(prev[i]-1== u_find(i)) {
                    docycle(prev[i]
-1,i, cost);
                    cycle 
= true;
                    
break;
                }
                
else
                {
                    u_union(i,prev[i]
-1);
                }

            }


        } 
while (cycle);

        
for (int i = 1; i < nvertex; i) {
            
if(!deleted[i])
            {
                cost
=g[prev[i]-1][i];
                
//cout<<"from "<<(prev[i]-1)<<" to "<
            }

        }

        printf(
"%.2f\n",cost);


    }


}



posted @ doubleh 阅读(2395) | |  

   今天买了本《算法概论》影印注释版,仔细看了第一章,果然名不虚传,很是吸引人。
第一章的习题难度适中,这里抽出第35题来,这题是证明wilson定理。
wilson定理:
   n是一个素数当且仅当: (n-1)! ≡ -1(mod n)

证明:
  首先我们证明若n是素数,那么等式成立,对于n=2,这是很明显的。以下证明n>2 的情形。
  1)若n是素数,那么关于n同余的乘法群g={1,2,3....n-1}
    群g每个元素都有逆元,映射 f:a -> a^-1 ,f(a)=a^-1 是一个一一映射。现在,任取一个元素,它的逆元要么是其它一个元素,或者是它本身。 我们假设其中元素x的逆元是它本身,那么x*x ≡1(mod n) =>(x 1)*(x-1)=k*n,而n是素数,所以要么x=n-1,要么x=1。也就是说,除了这两个元素,其它的元素的逆元都是映射到别的元素的。而n>2,是奇数,所以元素共有n-1个,也就是偶数个元素。这样,除了1和n-1外剩余的n-3个元素刚好结成两两一对(x,y)(逆元是唯一的)使得f(x)=y=x^-1,也就是xy≡1(mod n).
现在把g的元素全部乘起来,让互为逆元的元素组成一对那么(n-1)!=1*(n-1)*(x1*y1)*(x2*y2)...(xk*yk)≡1*(n-1)*1*1...1 (mod n)≡-1(mod n).
    这样,我们证明了一个方向了,下面我们证明另一个方向

  2)若(n-1)! ≡ -1(mod n)成立,则n是素数,若不然,令 n是和数。则(n-1)!肯定整除n,因为n的每个因子p满足p>1,p    现在令(n-1)!=k1*n.又(n-1)! ≡ -1(mod n) => k1*n ≡ -1(mod n),由同余性质知,存在k2使得k1*n 1=k2*n,两边同时除以n得k1 1/n=k2.显然,这是不可能的,所以若(n-1)! ≡ -1(mod n)成立,则n是素数

证明完毕!

这里用群的概念能够简化一定的描述,其实可以完全不用群的概念的,只不过这样一来,描述更长点,繁琐点!





posted @ doubleh 阅读(1990) | |  

2011年12月1日

如题:求连续正整数使得其和为给定的一个正整数
下面给出我的解法,几乎可以一步到位求出来
实现代码如下:
/**
*author: koth (
http://weibo.com/yovn)
*date:  2011-12-01
*/
#include 
<stdlib.h>
#include 
<stdio.h>
#include 
<stdint.h>

int solve(int y,int& x){
    
int m=0;
    
int t=y;
    
if(y<=0){
        x
=y;
        
return 1;
    }
    
while((t&1)==0){
        m
=1;
        t
=t>>1;
    }
    
if(m==32){
        x
=y;
        
return 1;
    }
    
int lastk=32;
    
for(;lastk>m1;lastk--){
        
if(y &(1<<(lastk-1))){
            
            
break;
        }
            
    }

    
//its a number st. exp(2,k)
    if(lastk==(m1)){
        x
=y;
        
return 1;
    }
    
int k=1<<(m1);
    
int b=(y>>m)-(1<<(lastk-m-1));

    x
=(1<<(lastk-m-2))(b1-k)/2;

    
if(x<=0){
        k
=k-1-((0-x)<<1);
        x
=0-x1;
    }
    
    
return k;

}

int main(int argc,char* argv[]){
    
if(argc<=1){
        fprintf(stdout,
"usage:%s number\n",argv[0]);
        
return 0;
    }
    
int y=atoi(argv[1]);
    
int x=0;
    
int k=solve(y,x);
    fprintf(stdout,
"%d=",y);
    
for(int i=0;i<k;i){
        fprintf(stdout,
"%d",xi);
        
if(i<(k-1)){
            fprintf(stdout,
"");
        }
    }
    fprintf(stdout,
"\n");
    
return 0;
}
posted @ doubleh 阅读(1749) | |  

2011年2月6日

     摘要: 年过的差不多了,今天偶尔兴起上hoj上翻几道dp练手的题来。。。,顺便把代码贴下留念  1.数塔 code highlighting produced by actipro codehighlighter (freeware) http://www.codehighlighter.com/ -->/**  *   */ pack...  
posted @ doubleh 阅读(1992) | |  

2009年12月4日

     摘要: 前一篇博客,我简单提了下怎么为nio2增加transmitfile支持,文件传送吞吐量是一个性能关注点,此外,并发连接数也是重要的关注点。 不过jdk7中又一次做了简单的实现,不支持同时投递多个acceptex请求,只支持一次一个,返回后再投递。这样,客户端连接的接受速度必然大打折扣。不知道为什么sun会做这样的实现,wsasend()/wsareceive()一次只允许一个还是可以理解,...  
posted @ doubleh 阅读(3864) | |  

2009年11月29日

jdk7的nio2特性或许是我最期待的,我一直想基于它写一个高性能的java http server.现在这个想法终于可以实施了。
本人基于目前最新的jdk7 b76开发了一个http server性能确实不错。
在windows平台上nio2采用accpetex来异步接受连接,并且读写全都关联到iocp完成端口。不仅如此,为了方便开发者使用,连iocp工作线程都封装好了,你只要提供线程池就ok。

但是要注意,iocp工作线程的线程池必须是 fix的,因为你发出的读写请求都关联到相应的线程上,如果线程死了,那读写完成情况是不知道的。

作为一个http server,传送文件是必不可少的功能,那一般文件的传送都是要把程序里的buffer拷贝到内核的buffer,由内核发送出去的。windows平台上为这种情况提供了很好的凯发天生赢家一触即发官网的解决方案,使用transmitfile接口

bool transmitfile(
    socket hsocket,
    handle hfile,
    dword nnumberofbytestowrite,
    dword nnumberofbytespersend,
    lpoverlapped lpoverlapped,
    lptransmit_file_buffers lptransmitbuffers,
    dword dwflags
);

你只要把文件句柄发送给内核就行了,内核帮你搞定其余的,真正做到zero-copy.
但是很不幸,nio2里asynchronoussocketchannel没有提供这样的支持。而为http server的性能考量,本人只好自己增加这个支持。

要无缝支持,这个必须得表现的跟 read /write一样,有完成的通知,通知传送多少数据,等等。

仔细读完sun的iocp实现以后发现这部分工作他们封装得很好,基本只要往他们的框架里加东西就好了。
为了能访问他们的框架代码,我定义自己的transmitfile支持类在sun.nio.ch包里,以获得最大的权限。

package sun.nio.ch;

import java.io.ioexception;
import java.lang.reflect.field;
import java.nio.channels.asynchronouscloseexception;
import java.nio.channels.asynchronoussocketchannel;
import java.nio.channels.closedchannelexception;
import java.nio.channels.completionhandler;
import java.nio.channels.notyetconnectedexception;
import java.nio.channels.writependingexception;
import java.util.concurrent.future;


/**

 * 
@author yvon
 * 
 
*/
public class windowstransmitfilesupport {
   
   //sun's nio2 channel  implementation class
    
private windowsasynchronoussocketchannelimpl channel;
   
    //nio2 framework core data structure
    pendingiocache iocache;

   //some field retrieve from sun channel implementation class
    
private object writelock;
    
private field writingf;
    
private field writeshutdownf;
    
private field writekilledf; // f

    windowstransmitfilesupport()
    {
        
//dummy one for jni code
    }

    
/**
     * 
     
*/
    
public windowstransmitfilesupport(
            asynchronoussocketchannel
             channel) {

        
this.channel = (windowsasynchronoussocketchannelimpl)channel;
        
try {
        // initialize the fields
            field f 
= windowsasynchronoussocketchannelimpl.class
                    .getdeclaredfield(
"iocache");
            f.setaccessible(
true);
            iocache 
= (pendingiocache) f.get(channel);
            f 
= asynchronoussocketchannelimpl.class
                    .getdeclaredfield(
"writelock");
            f.setaccessible(
true);
            writelock 
= f.get(channel);
            writingf 
= asynchronoussocketchannelimpl.class
                    .getdeclaredfield(
"writing");
            writingf.setaccessible(
true);

            writeshutdownf 
= asynchronoussocketchannelimpl.class
                    .getdeclaredfield(
"writeshutdown");
            writeshutdownf.setaccessible(
true);

            writekilledf 
= asynchronoussocketchannelimpl.class
                    .getdeclaredfield(
"writekilled");
            writekilledf.setaccessible(
true);

        } 
catch (nosuchfieldexception e) {
            
// todo auto-generated catch block
            e.printstacktrace();
        } 
catch (securityexception e) {
            
// todo auto-generated catch block
            e.printstacktrace();
        } 
catch (illegalargumentexception e) {
            
// todo auto-generated catch block
            e.printstacktrace();
        } 
catch (illegalaccessexception e) {
            
// todo auto-generated catch block
            e.printstacktrace();
        }
    }

    
    
/**
     * implements the task to initiate a write and the handler to consume the
     * result when the send file completes.
     
*/
    
private class sendfiletask<v, a> implements runnable, iocp.resulthandler {
        
private final pendingfuture<v, a> result;
        
private final long file;//file is windows file handle

        sendfiletask(
long file, pendingfuture<v, a> result) {
            
this.result = result;
            
this.file = file;
        }

    

        @override
        
// @suppresswarnings("unchecked")
        public void run() {
            
long overlapped = 0l;
            
boolean pending = false;
            
boolean shutdown = false;

            
try {
                channel.begin();

        

                
// get an overlapped structure (from the cache or allocate)
                overlapped = iocache.add(result);
                
int n = transmitfile0(channel.handle, file, overlapped);
                
if (n == iostatus.unavailable) {
                    
// i/o is pending
                    pending = true;
                    
return;
                }
                
if (n == iostatus.eof) {
                    
// special case for shutdown output
                    shutdown = true;
                    
throw new closedchannelexception();
                }
                
// write completed immediately
                throw new internalerror("write completed immediately");
            } 
catch (throwable x) {
                
// write failed. enable writing before releasing waiters.
                channel.enablewriting();
                
if (!shutdown && (x instanceof closedchannelexception))
                    x 
= new asynchronouscloseexception();
                
if (!(x instanceof ioexception))
                    x 
= new ioexception(x);
                result.setfailure(x);
            } 
finally {
                
// release resources if i/o not pending
                if (!pending) {
                    
if (overlapped != 0l)
                        iocache.remove(overlapped);
                
                }
                channel.end();
            }

            
// invoke completion handler
            invoker.invoke(result);
        }

        

        
/**
         * executed when the i/o has completed
         
*/
        @override
        @suppresswarnings(
"unchecked")
        
public void completed(int bytestransferred, boolean caninvokedirect) {
    

            
// release waiters if not already released by timeout
            synchronized (result) {
                
if (result.isdone())
                    
return;
                channel.enablewriting();

                result.setresult((v) integer.valueof(bytestransferred));

            }
            
if (caninvokedirect) {
                invoker.invokeunchecked(result);
            } 
else {
                invoker.invoke(result);
            }
        }

        @override
        
public void failed(int error, ioexception x) {
            
// return direct buffer to cache if substituted
        

            
// release waiters if not already released by timeout
            if (!channel.isopen())
                x 
= new asynchronouscloseexception();

            
synchronized (result) {
                
if (result.isdone())
                    
return;
                channel.enablewriting();
                result.setfailure(x);
            }
            invoker.invoke(result);
        }

    }

    
public <extends number, a> future<v> sendfile(long file, a att,
            completionhandler
<v, ? super a> handler) {

        
boolean closed = false;
        
if (channel.isopen()) {
            
if (channel.remoteaddress == null)
                
throw new notyetconnectedexception();

            
            
// check and update state
            synchronized (writelock) {
                
try{
                
if (writekilledf.getboolean(channel))
                    
throw new illegalstateexception(
                            
"writing not allowed due to timeout or cancellation");
                
if (writingf.getboolean(channel))
                    
throw new writependingexception();
                
if (writeshutdownf.getboolean(channel)) {
                    closed 
= true;
                } 
else {
                    writingf.setboolean(channel, 
true);
                }
                }
catch(exception e)
                {
                    illegalstateexception ise
=new illegalstateexception(" catch exception when write");
                    ise.initcause(e);
                    
throw ise;
                }
            }
        } 
else {
            closed 
= true;
        }

        
// channel is closed or shutdown for write
        if (closed) {
            throwable e 
= new closedchannelexception();
            
if (handler == null)
                
return completedfuture.withfailure(e);
            invoker.invoke(channel, handler, att, 
null, e);
            
return null;
        }



        
return implsendfile(file,att,handler);
    }


    
<extends number, a> future<v> implsendfile(long file, a attachment,
            completionhandler
<v, ? super a> handler) {
        
// setup task
        pendingfuture<v, a> result = new pendingfuture<v, a>(channel, handler,
                attachment);
        sendfiletask
<v,a> sendtask=new sendfiletask<v,a>(file,result);
        result.setcontext(sendtask);
        
// initiate i/o (can only be done from thread in thread pool)
        
// initiate i/o
        if (iocp.supportsthreadagnosticio()) {
            sendtask.run();
        } 
else {
            invoker.invokeonthreadinthreadpool(channel, sendtask);
        }
        
return result;
    }
    
    
private native int transmitfile0(long handle, long file,
            
long overlapped);
    
}

这个操作跟默认实现的里的write操作是很像的,只是最后调用的本地方法不一样。。

接下来,我们怎么使用呢,这个类是定义在sun的包里的,直接用的话,会报illegalaccesserror,因为我们的类加载器跟初始化加载器是不一样的。
解决办法一个是通过启动参数-xbootclasspath,让我们的包被初始加载器加载。我个人不喜欢这种办法,所以就采用jni来定义我们的windows transmitfile支持类。

这样我们的工作算是完成了,注意,发送文件的时候传得是文件句柄,这样做的好处是你可以更好的控制,一般是在发送前,打开文件句柄,完成后在回调通知方法里关闭文件句柄。



有兴趣的同学可以看看我的http server项目:


目前基本功能实现得差不多,做了些简单的测试,性能比较满意。这个服务器不打算支持servlet api,基本是专门给做基于长连接模式通信的定做的。






posted @ doubleh 阅读(2565) | |  

2009年5月4日

问题:
有个链表(list),有n个元素,当n很大的时候,我们通常想分批处理该链表。假如每次处理m条(0
问题很简单,我们需要次,这里我们用<>表示向上取整,[]表示向下取整,那么怎么来表示这个值呢?
我们可以证明:
=[(n-1)/m] 1    (0
不失一般性,我们设n=mk r(0<=r 1)当r>0时,

左边:=<(mk r)/m>==k =k 1
右边:[(n-1)/m] 1=[(mk r-1)/m] 1=[k (r-1)/m] 1=k 1 [(r-1)/m]=k 1
2)当r=0
左边:=k
右边:[(n-1)/m] 1=[(mk-1)/m] 1=[(m(k-1) m-1)/m] 1=[k-1 (m-1)/m] 1=k [(m-1)/m]=k

命题得证。

有了这个公式,我们在java代码里可以这样计算:
int nn=(n-1)/1
.


因为'/'是往下取整的。








posted @ doubleh 阅读(3936) | |  

2009年4月24日

原题:

command network

description

after a long lasting war on words, a war on arms finally breaks out between littleken’s and knuthocean’s kingdoms. a sudden and violent assault by knuthocean’s force has rendered a total failure of littleken’s command network. a provisional network must be built immediately. littleken orders snoopy to take charge of the project.

with the situation studied to every detail, snoopy believes that the most urgent point is to enable littenken’s commands to reach every disconnected node in the destroyed network and decides on a plan to build a unidirectional communication network. the nodes are distributed on a plane. if littleken’s commands are to be able to be delivered directly from a node a to another node b, a wire will have to be built along the straight line segment connecting the two nodes. since it’s in wartime, not between all pairs of nodes can wires be built. snoopy wants the plan to require the shortest total length of wires so that the construction can be done very soon.

input

the input contains several test cases. each test case starts with a line containing two integer n (n ≤ 100), the number of nodes in the destroyed network, and m (m ≤ 104), the number of pairs of nodes between which a wire can be built. the next n lines each contain an ordered pair xi and yi, giving the cartesian coordinates of the nodes. then follow m lines each containing two integers i and j between 1 and n (inclusive) meaning a wire can be built between node i and node j for unidirectional command delivery from the former to the latter. littleken’s headquarter is always located at node 1. process to end of file.

output

for each test case, output exactly one line containing the shortest total length of wires to two digits past the decimal point. in the cases that such a network does not exist, just output ‘poor snoopy’.


一开始没仔细读题,一看以为是最小生成树呢,结果krusal算法上去wa了,prim算法也wa,修修改改一直wa,第二天发现本题是要在有向图上面构造最小树形图。

按照著名的zhu-liu算法,仔细实现了一边,终于ac了。
按照我的理解总结下该算法,该算法对每个结点,除根节点外寻找最小入边,
1)如果这些入边不构成环,那么容易证明这些边构成最小树形图。
证明:设加上根节点r一共n个点,那么一共有n-1条边,证明从r能到每个点,若存在一点v,使得从r到v没有路径,那么,假设从v反向回退必然构成环,因为每个点除了r都有入边,如果不构成环,该路径可以无穷大。
2)如果存在环,我们把环收缩成一个点,更新相应的入边和出边,得到新的图g',使得原问题在g'中等效:
怎么收缩呢?
假设我们把环收缩成环上的任意一点v,所有进环的边和出环的边自动变成v的边(如果已有,取长度最短的),其余点标记为删除,更新不在环上的所有点进入该环的长度cost为cost-cost(prev[x],x);其中点x为进入环的边在环上的端点。出边保持不变。

这里为什么这么更新?因为这样更新使得我们的算法在新图上是等效的。任何环的解决后意味着在新图里面得为改收缩后的点寻找新的入边,而实际的花费应该是新的入边减去原有的入边长度,我们的算法在找到环的时候就把环上所有的边的长度计算在花费内了.而对出边是没有影响的。


到这算法的框架基本出来了。当为某点没找到入边的时候,意味着无解。为了加快无解的侦测,我们先运行一遍dfs搜索,假如从根节点出发,可触及的节点数小于n-1(不含r)则意味着无解。反之,肯定有解。
为什么?
因为如果可触及数小于n-1,意味着某点是不可触及的,也就是原图不是弱连通。对该点来说不存在从r到它的路径。反之,从r到某点都有一条路径,沿着该路径就能找到结点的入边。


第二个问题是,如何快速侦测环呢?
我使用了一个不相交集。回忆krusal的算法实现里面也是使用不相交集来避免找产生环的最小边。

下面是我的代码:

// 3164.cpp : defines the entry point for the console application.
//

#include 
<iostream>
#include 
<cmath>


using namespace std;

typedef 
struct _point
{

    
double x;
    
double y;

    
double distanceto(const struct _point& r)
    {
          
return sqrt((x-r.x)*(x-r.x)(y-r.y)*(y-r.y));

    }

}point;



const int max_v=100;
const int max_e=10000;
const double no_edge=1.7976931348623158e 308;

point vertexes[max_v]
={0};
int parents[max_v]={0};
int ranks[max_v]={0};
double g[max_v][max_v]={0};
bool visited[max_v]={0};
bool deleted[max_v]={0};
int prev[max_v]={0};

int nvertex=0;
int nedge=0;




int u_find(int a)
{
    
if(parents[a]==a)return a;
    parents[a]
=u_find(parents[a]);
    
return parents[a];
}
void u_union(int a,int b)
{
    
int pa=u_find(a);
    
int pb=u_find(b);
    
if(ranks[pa]==ranks[pb])
    {
        ranks[pa]
;
        parents[pb]
=pa;
    }
else if(ranks[pa]<ranks[pb])
    {
        parents[pa]
=pb;
    }
    
else
    {
        parents[pb]
=pa;
    }
}

void dfs(int v,int& c)
{

    visited[v]
=true;
    
for(int i=1;i<nvertex;i)
    {
        
if(!visited[i]&&g[v][i]<no_edge)
        {
            c
=1;

            dfs(i,c);
        }

    }

}



void docycle(int s,int t,double& cost)
{
    memset(visited,
0,sizeof(bool)*nvertex);
    
int i=s;
    
do
    {
        visited[i]
=true;
        cost
=g[prev[i]-1][i];
        
//cout<<"from "<<(prev[i]-1)<<" to "<
        i=prev[i]-1;

    }
while(i!=s);


    
do
    {
        

        
for(int k=0;k<nvertex;k)
        {
            
if(!deleted[k]&&!visited[k])
            {
            
                
if(g[k][i]<no_edge)
                {
                    
if(g[k][i]-g[prev[i]-1][i]<g[k][s])
                    {


                        g[k][s]
=g[k][i]-g[prev[i]-1][i];
                        
//cout<<"1.update ["<
                    }
                    

                }
                
if(g[i][k]<no_edge)
                {
                    
if(g[i][k]<g[s][k])
                    {

                        g[s][k]
=g[i][k];
                        
//cout<<"2.update ["<
                    }
                    
                }
            }
        }


        
if(i!=s)
        {
            deleted[i]
=true;
            
//cout<<"mark "<
        }
        i
=prev[i]-1;


    }
while(i!=s);



}




int main(void)
{



    
    
while(cin>>nvertex>>nedge)
    {

        
int s,t;

        
int nv=0;
        
bool cycle=0;
        
double cost=0;
        memset(vertexes,
0,sizeof(vertexes));
        memset(visited,
0,sizeof(visited) );

        memset(deleted,
0,sizeof(deleted));
        memset(g,
0,sizeof(g));
        memset(prev,
0,sizeof(prev));


        memset(ranks,
0,sizeof(ranks));

        memset(parents,
0,sizeof(parents));

        
for(int i=0;i<nvertex;i)
        {

            cin
>>vertexes[i].x>>vertexes[i].y;
            parents[i]
=i;
            
for(int j=0;j<nvertex;j)
            {
                g[i][j]
=no_edge;
            }

        }
        
for(int i=0;i<nedge;i)
        {

            cin
>>s>>t;
            
if(t==1||s==t)continue;
            g[s
-1][t-1]=vertexes[s-1].distanceto(vertexes[t-1]);

        }



        dfs(
0,nv);
        
if(nv<nvertex-1)
        {

            cout
<<"poor snoopy"<<endl;
            
continue;
        }



        
do {
            cycle
=false;
            
for(int i=0;i<nvertex;i){parents[i]=i;}
            memset(ranks,
0,sizeof(bool)*nvertex);
        

            
for (int i = 1; i < nvertex; i) {
                
double minimum=no_edge;
                
if(deleted[i])continue;
                
for(int k=0;k<nvertex;k)
                {
                    
if(!deleted[k]&&minimum>g[k][i])
                    {
                        prev[i]
=k1;
                        minimum
=g[k][i];
                    }
                }
                
if(minimum==no_edge)
                {
                

                    
throw 1;
                }
                
if (u_find(prev[i]-1== u_find(i)) {
                    docycle(prev[i]
-1,i, cost);
                    cycle 
= true;
                    
break;
                }
                
else
                {
                    u_union(i,prev[i]
-1);
                }

            }


        } 
while (cycle);

        
for (int i = 1; i < nvertex; i) {
            
if(!deleted[i])
            {
                cost
=g[prev[i]-1][i];
                
//cout<<"from "<<(prev[i]-1)<<" to "<
            }

        }

        printf(
"%.2f\n",cost);


    }


}



posted @ doubleh 阅读(2395) | |  

2009年4月10日

   今天买了本《算法概论》影印注释版,仔细看了第一章,果然名不虚传,很是吸引人。
第一章的习题难度适中,这里抽出第35题来,这题是证明wilson定理。
wilson定理:
   n是一个素数当且仅当: (n-1)! ≡ -1(mod n)

证明:
  首先我们证明若n是素数,那么等式成立,对于n=2,这是很明显的。以下证明n>2 的情形。
  1)若n是素数,那么关于n同余的乘法群g={1,2,3....n-1}
    群g每个元素都有逆元,映射 f:a -> a^-1 ,f(a)=a^-1 是一个一一映射。现在,任取一个元素,它的逆元要么是其它一个元素,或者是它本身。 我们假设其中元素x的逆元是它本身,那么x*x ≡1(mod n) =>(x 1)*(x-1)=k*n,而n是素数,所以要么x=n-1,要么x=1。也就是说,除了这两个元素,其它的元素的逆元都是映射到别的元素的。而n>2,是奇数,所以元素共有n-1个,也就是偶数个元素。这样,除了1和n-1外剩余的n-3个元素刚好结成两两一对(x,y)(逆元是唯一的)使得f(x)=y=x^-1,也就是xy≡1(mod n).
现在把g的元素全部乘起来,让互为逆元的元素组成一对那么(n-1)!=1*(n-1)*(x1*y1)*(x2*y2)...(xk*yk)≡1*(n-1)*1*1...1 (mod n)≡-1(mod n).
    这样,我们证明了一个方向了,下面我们证明另一个方向

  2)若(n-1)! ≡ -1(mod n)成立,则n是素数,若不然,令 n是和数。则(n-1)!肯定整除n,因为n的每个因子p满足p>1,p    现在令(n-1)!=k1*n.又(n-1)! ≡ -1(mod n) => k1*n ≡ -1(mod n),由同余性质知,存在k2使得k1*n 1=k2*n,两边同时除以n得k1 1/n=k2.显然,这是不可能的,所以若(n-1)! ≡ -1(mod n)成立,则n是素数

证明完毕!

这里用群的概念能够简化一定的描述,其实可以完全不用群的概念的,只不过这样一来,描述更长点,繁琐点!





posted @ doubleh 阅读(1990) | |  

2009年1月17日





今天去网上看了一下09年的考研试题,看见该题目(图片):



先来定义结点(为了简便,省略set/get):
public class node
{
 
public int data;
 
public node link;
}
我能想到的两种解法,一个基于递归:

递归版的思路就是,基于当前结点,如果后一个是倒数第k-1,那么当前结点是所求,若不然,返回当前是倒数第几个。
public int printrkwithrecur(node head,int k)
    {
        
if(k==0||head==null||head.link==null)return 0;
        
if(_recurfind(head.link,k)>=k)return 1;
        
return 0;
    }
    
private final int _recurfind(node node, int k) {
        
if(node.link==null)
        {
            
return 1;
        }
        
int sret=_recurfind(node.link,k);
        
if(sret==k-1)
        {
            system.out.println(
"got:"node.data);
            
return k;
        }
        
return sret1;
    }


对每个结点,该算法都只访问一次,因此复杂度o(n).

第二解法,相对递归来说,这种方法可以算是消除递归版,而且从某种意义上来说比递归更高效,跟省空间,递归版实际上是把回溯的数据存在栈上,而版方法是自己存储,且利用数组实现一个循环队列,只存储k个元素。

public static class cycleintqueue
    {
        
int[] datas;
        
int top=0;
        
int num=0;
        
public cycleintqueue(int n)
        {
            datas
=new int[n];
        }
        
        
public void push(int i)
        {
            datas[(top
)%datas.length]=i;
            num
;
            
        }
        
public int numpushed()
        {
            
return num;
        }
        
        
        
public int getbuttom()
        {
            
return datas[top%datas.length];
        }
    }
    
public int printrkwithcyclequeue(node head,int k)
    {
        
if(k==0||head==null)return 0;
        cycleintqueue queue
=new cycleintqueue(k);
        node cur
=head.link;
        
while(cur!=null)
        {
            queue.push(cur.data);
            cur
=cur.link;
        }
        
if(queue.numpushed()<k)return 0;
        
        system.out.println(
"got:"queue.getbuttom());
        
return 1;
    }

本算法,都每个结点也只放一次,另外进行一次入队操作,该操作复杂度o(1),从而,整个算法复杂度仍是o(n).


posted @ doubleh 阅读(2252) | |  

2009年1月6日


近日在csdn上看到中软一道面试题,挺有意思的。
题目:一条小溪上7块石头,如图所示:

分别有六只青蛙:a,b,c,d,e,f。a,b,c三只蛙想去右岸,它们只会从左向右跳;d,e,f三只蛙想去左岸,它们只会从右向左跳。青蛙每次最多跳到自己前方第2块石头上。请问最少要跳几次所有青蛙上岸。写出步骤。

这个题是个路径搜索的问题,在解空间搜索所有的解,并找出最优的解法(即步骤最少的)。
那么怎么算是一个解呢?具体而言就是最后石头上没有青蛙了。



我们先给题目建模,7块石头,其上可以是没青蛙,可以有一只往左跳的青蛙,也可以有一只往右跳的青蛙。可以把这7块石头看成一个整体,来表示一个状态。这里我们把这7块石头看成一个数组,里面只能有0,1,2三种值,这样表示,那么初始时为:
1,1,1,0,2,2,2
我们把它再表示成一个数字,来表示状态值,这个值把这个数组按三进制拼成一个数字,我们用一个辅助函数来做这件事情:
private final int makes() {
        
        
int r=0;
        
int p=1;
        
for(int i=0;i<7;i)
        {
            r
=p*states[i];
            p
*=3;
        }
        
return r;
    }

那么题目现在变成从状态111022转换成状态0000000,所需最少的步骤.

那么状态是怎样转换的呢?
很显然。,每次青蛙跳都会触发状态的转换,我们在每个状态时搜索每种可能的转换,我们记初始状态为s(s等于三进制111022)记要求解的值为opt(s),假如可以转换到t1,t2,...tk.
那么,显然
opt(s)=min(1opt(t1),1opt(t2),.,1opt(tk));

另外,由于最终状态为0,所以opt(0)=0,就是说已经在最终状态了,就不需要一步就可以了。
有了上面这个等式,我们可以递归求解了,但是如果单纯的递归,会导致大量的重复计算,所以这里我们用备忘录的方法,记下已经求解出来的opt(x),放在一个数组里,由于只有7块石头,所以最多我们需要3^7=2187个状态。我们用一个2187个元素的数组,  其中第i个元素表示opt(i),初始化每个元素用-1表示还未求解。opt(0) 可直接初始化为0.

到此我们还有一个问题,怎么能够在算法结束的时候打印出最优的步骤呢?按照这个步骤,我们可以重建出青蛙是如何在最优的情况下过河的。为此,我们可以再用一个步骤数组,每次在采取最优步骤的时候记录下来。

整个算法如下:
package test;

import java.util.arrays;
/**
 *
 * @author yovn
 *
 */
public class frogjump {
   
    private int steps[];
    private int states[];
   
   
    private static class step
    {
        int offset=-1;
        int jump;
        int jumpto;
    }
   
   
    private step jumps[];
    private int inits;
    public frogjump()
    {
        steps=new int[81*27];
        states=new int[7];
        for(int i=0;i<3;i )states[i]=1;
        for(int i=4;i<7;i )states[i]=2;
        arrays.fill(steps, -1);
        steps[0]=0;
        jumps=new step[81*27];
        inits=makes();
    }
   
    public int shorteststeps(int s)
    {
        if(steps[s]==-1)
        {

            int minstep=integer.max_value;
            step onestep=new step();
            for(int i=0;i<7;i )
            {
                if(states[i]==1)
                {
                    if(i>4)
                    {
                        states[i]=0;
                        minstep = recurfind(minstep,onestep,i,7-i);
                        states[i]=1;
                    }
                    else
                    {
                        if(states[i 1]==0)
                        {
                            states[i]=0;
                            states[i 1]=1;
                            minstep = recurfind(minstep,onestep,i,1);
                            states[i]=1;
                            states[i 1]=0;
                           
                        }
                        if(states[i 2]==0)
                        {
                            states[i]=0;
                            states[i 2]=1;
                            minstep = recurfind(minstep,onestep,i,2);
                            states[i]=1;
                            states[i 2]=0;
                           
                        }
                    }
                }
                else if(states[i]==2)
                {
                    if(i<2)
                    {
                        states[i]=0;
                       
                        minstep = recurfind(minstep,onestep,i,-1-i);
                        states[i]=2;
                    }
                    else
                    {
                        if(states[i-1]==0)
                        {
                            states[i]=0;
                            states[i-1]=2;
                            minstep = recurfind(minstep,onestep,i,-1);
                            states[i]=2;
                            states[i-1]=0;
                           
                        }
                        if(states[i-2]==0)
                        {
                            states[i]=0;
                            states[i-2]=2;
                            minstep = recurfind(minstep,onestep,i,-2);
                            states[i]=2;
                            states[i-2]=0;
                           
                        }
                    }
                }
               
            }
            steps[s]=minstep;
            jumps[s]=onestep;
           
           
        }
        return steps[s];

    }

    private final int recurfind(int minstep, step onestep, int pos, int jump) {
        int tos=makes();
        int r=shorteststeps(tos);
        if(r         {
            onestep.jump=jump;
            onestep.offset=pos;
            onestep.jumpto=tos;
            minstep=r 1;
        }
        return minstep;
    }

   
   
    public void printpath()
    {
        int s=inits;
        int i=1;
       
        while(s!=0)
        {
           
           
            system.out.println("[" (i ) "] frog at #" jumps[s].offset " jumps #" jumps[s].jump);
            s=jumps[s].jumpto;
           
        }
    }
    private final int makes() {
       
        int r=0;
        int p=1;
        for(int i=0;i<7;i )
        {
            r =p*states[i];
            p*=3;
        }
        return r;
    }

    /**
     * @param args
     */
    public static void main(string[] args) {
        frogjump fj=new frogjump();
        int steps=fj.shorteststeps(fj.inits);
       
        system.out.println("use " steps " steps!");
        fj.printpath();

    }

}

运行结果:

use 21 steps!
[
1] frog at #2 jumps #1
[
2] frog at #4 jumps #-2
[
3] frog at #5 jumps #-1
[
4] frog at #3 jumps #2
[
5] frog at #1 jumps #2
[
6] frog at #0 jumps #1
[
7] frog at #2 jumps #-2
[
8] frog at #0 jumps #-1
[
9] frog at #4 jumps #-2
[
10] frog at #2 jumps #-2
[
11] frog at #0 jumps #-1
[
12] frog at #5 jumps #2
[
13] frog at #3 jumps #2
[
14] frog at #1 jumps #2
[
15] frog at #5 jumps #2
[
16] frog at #3 jumps #2
[
17] frog at #5 jumps #2
[
18] frog at #6 jumps #-1
[
19] frog at #5 jumps #-2
[
20] frog at #3 jumps #-2
[
21] frog at #1 jumps #-2








posted @ doubleh 阅读(4004) | |  

2008年12月10日

写一个函数,输出前n个数(从7开始),这n个数满足如下3个条件中的任意一个
1.整出7
2.各位上的数字之和整除7,(比如34)
3.任意位上包含数字7


附我的代码:
void printn(int n)
{

    
    
int c=0;
    
int i=7;
    
do 
    {
        
if(i%7 ==0)
        {
            printf(
"%d\n",i);
            c
;
        }
        
else
        {
            
int j=i%10;
            
int k=j;
            
int s=k;
            
int p=10;
            
while(k             {

                
if(j==7)
                {
                    printf(
"%d\n",i);
                    s
=0;
                    c
;
                    
break;

                }
                
else
                {
                    j
=((i-k)/p)%10;
                    s
=j;
                    k
=j*pk;
                    p
*=10;


                }
            }
            
if(s&&s%7==0)
            {


                printf(
"%d\n",i);
                c
;
            }
            

        }
        i
;
    } 
while (c<n);
}


posted @ doubleh 阅读(3190) | |  

仅列出标题  下一页
网站地图