2042 Lagranges Four-Square Theorem

2042 Lagranges Four-Square Theorem

問題

解答方針

問題自体は単純ですが,時間が非常に厳しいです.平方数および二つの平方数の和で表せる数をあらかじめHashSetに登録しておくことで高速化しています.

解答例

import java.util.*;

public class Main {
    public static void main(String[] args) {
        Scanner sc = new Scanner(System.in);
        Solver sol = new Solver();
        while(true){
            int n = sc.nextInt();
            if(n==0) break;
            System.out.println(sol.solve(n,4,1));
        }
    }
}

class Solver{
    HashSet<Integer> square;
    HashSet<Integer> square2;
    Solver(){
        square = new HashSet<Integer>();
        square2 = new HashSet<Integer>();
        for(int i=1;i<=32768;i++){
            if(presolve(i,1,1)>0) square.add(i);
        }
        for(int i=1;i<=32768;i++){
            if(presolve(i,2,1)>0) square2.add(i);
        }
    }
    
    public int presolve(int n,int m,int a){
        // n==0 m>=0
        if(n==0) return 1;
        // n>0 m==0
        else if(m==0) return 0;
        // n>0 m>0
        else{
            int cnt = 0;
            for(int i=a;i*i<=n;i++){
                cnt += solve(n-i*i,m-1,i);
            }
            return cnt;
        }
    }
    
    public int solve(int n,int m,int a){
        // n==0 m>=0
        if(n==0) return 1;
        // n>0 m==0
        else if(m==0) return 0;
        // n>0 m>0
        else if(m==1){
            if(!square.contains(n)) return 0;
            else if(n>=a*a) return 1;
            else return 0;
        }
        else if(m==2){
            if(!square2.contains(n)) return 0;
            
            int cnt = 0;
            for(int i=a;i*i<=n;i++){
                cnt += solve(n-i*i,m-1,i);
            }
            return cnt;
        }
        else{
            int cnt = 0;
            for(int i=a;i*i<=n;i++){
                cnt += solve(n-i*i,m-1,i);
            }
            return cnt;
        }
    }
}

タグ:

+ タグ編集
  • タグ:

このサイトはreCAPTCHAによって保護されており、Googleの プライバシーポリシー利用規約 が適用されます。

最終更新:2006年04月17日 02:05
ツールボックス

下から選んでください:

新しいページを作成する
ヘルプ / FAQ もご覧ください。